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FOREWORD 


This  final  report  under  Contract  No.  DAADO7-79-C-0122  with 
Atmospheric  Sciences  Laboratory  (ASL)  at  the  White  Sands  Missile  Range 
describes  computerized  techniques  for  the  automatic  calculation  of  line- 
of-sight  integrated  concentrations,  transmittance  and  probabi.ity  of 
detection  in  a  polluted  battlefield  environment.  Specific  provision  has 
been  made  for  automated  treatment  of  atmospheric  contaminant  sources 
typically  present  in  the  battlefield.  The  program  incorporates  many 
features  of  the  Smoke  Obscuration  Routine  and  EPAMS  programs  developed 
by  the  H.  E.  Crainer  Company  under  previous  contracts  with  ASL  and  relies 
on  source  description  material  and  source  models  contained  in  the  E-0 
SAEL  programs  compiled  by  ASL.  This  report  consists  of  three  volumes. 
Volume  I  describes  the  various  model  components,  the  mathematical  rela¬ 
tionships  used  in  the  program,  and  the  procedures  used  to  relate  the 
various  program  components.  Volume  II  contains  users'  instructions  for 
operating  the  program.  Volume  III  contains  a  listing  of  the  program. 


Mr.  Robert  Umstead,  who  served  as  COR  for  ASL  under  the 
contract  until  the  spring  of  1980  and  his  replacement,  Mr.  John  Marrs, 
provided  valuable  technical  guidance  in  specifying  Army  requirements  and 
valuable  assistance  in  assembling  the  requisite  source  description  data 
and  source  models  used  in  the  program. 
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SECTION  I 
INTRODUCTION 


1.1  BACKGROUND 

The  Atmospheric  Sciences  Laboratory  (ASL) ,  White  Sands  Missile 
Range  (WSMR) ,  Nsw  Mexico,  has  been  directed  to  develop  plans  and  opera¬ 
tional  techniques  for  determining  atmospheric  effects  on  electro-optical/ 
millimeter  systems  operating  in  a  battlefield  environment.  The  atmosphere 
of  the  battlefield  is  polluted  by  smoke,  gasses  and  particles  from 
artillery  shell  and  bomb  bursts,  smoke  munitions,  vehicle  traffic  and 
other  battlefield  activities.  The  effectiveness  of  electro-optical  and 
near-millimeter  wavelength  devices  can  be  severely  affected  by  the  pre¬ 
sence  of  these  pollutants.  The  effects  of  obscurants  on  sensing  devices 
are  related  to  the  type  of  obscurant,  the  integrated  line-of-sight 
concentration  (CL),  the  response  of  the  sensor  to  attenuation  in  the 
sensor  wavelength  and  other  o^*'  :al  effects.  In  turn,  CL  is  directly 
related  to  the  amount  of  obscurant  released  to  the  atmosphere  and  the 
dispersion  and  transport  of  the  obscurant  in  the  atmosphere. 

ASL  has  been  conducting  and  sponsoring  major  efforts  in  the 
development  of  appropriate  meteorological  inputs  to  dispersion  and 
transport  models  and  improved  methods  for  accurately  and  realistically 
describing  the  transport  and  dispersion  of  natural  and  battle-induced 
contaminants.  These  efforts  include  the  implementation  of  automated 
procedures  for  estimating  cloud  trajectories  and  dispersion  characteristics 
downwind  from  sources  of  contaminate  material  in  regions  of  complex 
terrain.  The  EPAMS  program  (Durabauld  and  Bjorklund,  1977)  previously 
developed  for  ASL  contains  the  following  features  chat  are  retained  in 
the  program  described  in  this  report: 

•  Diagnostic  analysis  routines  for  calculating  the 
depth  of  the  atmospheric  surface  mixing  layer  and 
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Che  wind  velocity  within  this  layer  from  routine 
meteorological  measurements 

•  A  diagnostic  Mesoscale  Wind-Field  Model  that  uses  this 
information  to  calculate  the  wind  field  and  mixing- 
layer  depth  contours  above  regions  of  complex  terrain 

•  Diagnostic  application  routines  which  use  the  output 
from  the  mesoscale  model  routines  with  input  source 
characteristics  end  dispersion  models  to  calculate 
dosages  and  concentrations  of  contaminate  material 

•  Graphics  routines  to  produce  displays  of  the  wind- 
field  solutions,  mixing-layer  depth  contours  and 
dosage  and  concentration  isopleths  as  overlays  on 
the  complex  terrain  contours 

Under  Contract  No.  DAEA-18-77-C-0060,  the  H.  E.  Craraer  Company  added  the 
MS3  smoke  obscuration  routine  to  the  EPAMS  program  (Dumbauld,  Saterlie 
and  Cheney,  1979).  The  MS3  Routine  allows  the  program  user  to  calculate 
line-of-sighc  concentrations,  transmittance  and  probability  of  detection 
downwind  from  standard  chemical  smoke  sources  employed  by  the  U.  S.  Anuy. 
it  is  coupled  with  the  EPAMS  program  in  that  the  diagnostic  analysis 
routines  and  diagnostic  Mesoscale  Wind-Field  Model  are  used  to  calculate 
meteorological  inputs  representative  of  the  specific  battlefield  area  of 

interest  for  use  in  the  smoke  dispersion  models.  Many  of  the  routines 
used  in  the  MS 3  are  also  used  in  the  new  program.  As  explained  below,  the 
technical  objectives  of  the  work  described  in  this  report  are  a  logical 
continuation  of  the  overall  ASL  efforts  to  describe  atmospheric  effects 
on  EG  systems. 
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1.2 


PURPOSE 


The  primary  purpose  of  the  work  described  in  this  report  is  to 
provide  a  working  prototype  computerized  model  for  estimating  the  dis¬ 
persion  and  downwind  transport  of  contaminants  introduced  into  the  lower 
atmosphere  by  the  multiplicity  of  sources  typically  present  in  the 
battlefield  ares  (dust  and  smoke  from  exploding,  munitions,  moving  vehicles, 
muzzle  blasts;  smoke  from  burning  vegetation).  Further,  the  computer 
program  is  to  provide  a  means  of  estimating  the  line-of-sight  (LOS) 
concentration,  transmittance  and  probability  of  detection  as  affected  by 
the  transport  and  dispersion  of  the  contaminants.  The  program  is  to 
perform  these  functions  with  limited  meteorological  input  data  and  a 
minimum  of  user  specification  of  the  source  characteristics  required  by 
the  computer  program. 


1.3  THE  AMSORB  PROGRAM 

The  AMSORB  (Analysis  of  Multiple  Source  Obscurants  on  the 
Realistic  Battlefield)  prototype  computer  program  has  been  developed 
in  response  to  the  purposes  stated  in  Section  1.2  above.  Figure  1-1  is  a 
simplified  diagram  showing  the  functional  elements  of  the  AMSORB  program, 
which  is  based  on  the  previously-developed  EPAMS  program  functional 
concepts.  The  Analysis  Phase  is  called  by  the  AMSORB  Executive  as  a 
result  of  A  standing  request  for  a  routine  analysis  of  meteorological 
conditions  in  the  battlefield  surface  mixing  layer.  The  Analysis  Phase 
extracts  available  upper-air  and  surface  meteorological  data  f ro«a  the 
Data  Base.  As  shown  in  Figure  1-1,  the  Analysis  Phase  is  comprised  of 
two  major  elements.  The  Mixing-Layer  Analysis  Routine  calculates  initial¬ 
ization  parameters  for  the  Mesoscale  Wind-Field  Model,  as  well  as  some 
meteorological  inputs  for  direct  use  in  the  ApplieacU  Phase,  using 
Data  Base  Inputs.  The  Mesoscale  Wind-Field  Model  calculates  wind-velocity 
patterns  and  mixing  heights  for  use  in  the  crajeetory/cranaport  subeleaent 
of  the  Applications  Phase.  The  Analysis  Phase  does  not  eoeaunieace 
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ANALYSIS 

PHASE 


directly  with  the  Applications  Phase,  but  supplies  solutions  to  the  Data 
Base  for  subsequent  use.  The  Applications  Phase  is  called  by  the  AMSORB 
Executive  as  a  result  of  a  request  for  information  routed  through  the 
ABSORB  Executive.  The  Trajectory/Transport  Routine,  Dispersion  Model 
Routine  and.  the  MS3  Routine  in  the  Application  Phase  are  identict  .  to 
the  corresponding  routines  previously  included  in  the  EPAMS  prog  -ua  and 
can  be  called  and  used  in  the  same  maimer.  The  Battlefield  Kite'  onment 
Routine  is  the  new  addition  to  the  program  for  accomplishing  .e  calcu¬ 
lation  of  LOS  concentrations,  transmittance  and  probabil*ry  of  detection 
when  all  types  of  sources  are  present  on  the  battlefield.  The  Battle¬ 
field  Environment  Routine  cannot  be  called  by  the  user  mless  the  user 
has  first  requested,  through  the  AMSORB  Executive,  source  character¬ 
istics  for  the  battlefield  scenario  from  the  Data  Processor  Phase.  The 
Battlefield  Source  Characteristics  Routine  automatically  passes  the 
requisite  source  characteristics  to  the  Data  Base  depending  on  the  types 
of  sources  the  user  specifies.  The  source  characteristics  routine  has 
been  placed  as  a  separate  routine  in  the  AMSORB  program  to  facilitate 
updating  as  new  information  on  source  characteristics  becomes  available. 

1.4  ORGANIZATION  Of  THE  REPORT 

Many  of  the  automated  operations  of  the  AMSORB  Program  are 
based  on  elements  of  the  EPAMS  program  described  in  the  report  by  Du»- 
bauld  and  Bjorklund  (1977)  and  the  Smoke  Obscuration  Routine  fot  EPAMS 
described  by  Ourabauld,  Saterlie  and  Cheney  (1979).  To  provide  the  back¬ 
ground  information  required  to  understand  the  concepts  of  the  AMSORB 
program,  we  have  abstracted  or  reproduced  in  this  report  some  of  the 
original  reports.  Scetion  2  of  Volume  i  below  contains  a  description  of 
the  Meseseaie  Wind-Field  Model,  the  transport  and  dispersion  models,  and 
the  obscuration  (visibility  and  probability  of  detection)  models  used  in 
the  AMSORB  Program.  The  models  and  procedures  used  to  develop  source 
characteristics  for  the  Battlefield  Environment  Routine  are  described 


in  Section  3.  Section  4  describes  the  various  model  routines  in  the 
AMSQRB  Program  and  the  analysis  routines  and  other  procedures  used  to 
develop  model  input  parameters  for  these  routines.  The  computer  program 
structure  of  AMSORB  is  outlined  in  Section  5.  Section  6  describes 
the  results  of  example  calculations  made  using  the  AMSORB  program. 
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SECTION  2 


MESOSCAI.E  WIND-FIELD  MODEL,  TRAJECTORY/TRANSPQRT  ROUTINE, 
AND  DISPERSION  MODELS 


2.1  MESCSCALE  WIND-FIELD  MODEL 


The  Analysis  Phase  of  AMSORB  uses  a  numerical  modeling  tech¬ 
nique  developed  for  ASL  by  Tingle  and  Bjorklund  (1973)  to  calculate 
details  of  the  mesoscale  wind-field  and  depth  of  the  surface-mixing 
layer  above  complex  terrain.  This  modeling  technique  is  used  in  AMSORB 
for  the  same  purpose.  A  brief  description  of  the  technique,  abstracted 
from  the  report  by  Dumbauld  and  Bjorklund  (1977)  ,  follows. 


The  computer  algorithm,  based  on  the  shallow  -  fluid  equations 
of  oceanography,  describes  the  motion  of  fluid  in  two  dimensions  by  the 
expressions : 
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where 

u,v  a  x, y  components  of  the  iluiJ  velocity  in  the  layer  next  to 
the  terrain 


H  ■  depth  of  the  surface  mixing  layer 
fa 

2^  =  height  of  the  terrain 

g'  »  reduced  acceleration  of  gravity 
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acceleration  of  gravity 


P  =■  fluid  density  in  the  surface  mixing  layer 
0 

^  fluid  density  in  the  passive  layer  above  the  surface 
mixing  layer 

These  equations  describe  mesosca.le  wind-patterns  above  complex  terrain 
when  the  atmosphere  above  the  terrain  can  to  distinctively  divided  into 
two  layers  of  different  densities.  This  condition  typically  exists  when 
the  mixing  layer  above  the  terrain  is  capped  by  an  elevated  inversion. 
The  interaction  of  the  upper  layer  with  the  lower  layer  is  limited  to 
reducing  the  speed  of  gravity  wives  in  the  lower  layer.  Equations  (2-11 
(2-2)  and  (2-3)  are  written  in  their  momentum  form.  Mesoscale  wind 
patterns  are  generated  by  impulsively  accelerating  the  momentum  in  the 
lower  layer  to  a  constant  value  and,  after  an  empirically-determined 
model  time  (see  Section  5.1/  in  which  the  flow  in  the  region  of  interest 
is  assumed  to  become  stabilised,  using  the  steady-state  solution  tnus 
obtained . 


The  computer  algorithm  contains  an  explicit  Lax-Wendroff 
finite-differencing  procedure  using  nine  grid  points  and  one  time  level 
in  combination  with  a  nine-point  low-pass  filter.  The  finite  differ¬ 
ences  are  centered  to  preserve  high-order  accuracy.  The  grid  arrange¬ 
ment  and  finite-differencing  form  of  the  equations  art  described  in  the 
report  by  tingle  and  Bjorklund  (1973).  High-frequency  oscillations 
generated  In  critical  flow  situations  near  regions  where  hydraulic  jumps 
occur  arc  removed  by  the  filtering  technique,  vhieh  is  based  on  a  sug¬ 
gestion  by  Shapiro  (19/0)  and  has  the  advantage  chat  its  smooching 
effect  is  easily  determined  from  its  response  function.  In  the  Analysis 
Phase,  the  filter  is  applied  every  seventh  time  step.  It  is  well  known 
that  the  length  of  the  time  step  Ac  used  In  obtaining  the  solution 
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also  affects  the  stability  of  the  finite  difference  solution.  The  time 
step  is  altered  automatically  in  this  algorithm  in  accordance  with  the 
expression 


X 

s 

Ax  . 
min 

n 

jVu2  +  v2  +  /g  »m 

max 

where  the  stability  parameter  Xg  is  dependent  on  the  grid  spacing  Ax 
and  varies  between  zero  and  one.  for  the  present  application  which  uses 
a  minimum  spacing  (Axmin)  ^  kilometers,  the  stability  parameter  has 
been  set  equal  to  0.4.  Because  the  entire  solution  grid  must  be  searched 
to  determine  the  maximum  wave  speed  in  the  denominator  of  Equation  (2-4), 
die  value  of  At  is  calculated  every  second  time  step  rather  than  every 
time  step  to  save  computation  time. 

A  sensitivity  analysis  of  the  two-layer  Mesoscale  Wind-field 

Model  indicated  that  model  instability  can  frequently  occur  in  localized 

terrain  regions  where  the  inversion  height  H  is  initialized,  or  is 

m 

adjusted  by  die  model  to  a  few  tens  of  meters  above  the  higher  terrain 
elevations.  Vue  momentum  conservation  properties  of  the  computational 
scheme  arc  such  chat  wind  velocities  in  these  localized  regions  may  be¬ 
come  unreasonably  1  :ge  as  the  depth  of  the  lower  layer  approaches  a 
minimum  vstue.  In  order  to  prevent  this  type  of  solution  instability, 

.ue  u-  and  v-coapouents  of  the  solution  at  any  given  time  step  are  auto¬ 
matically  set  to  zero  and  du  layer  depth  set  equal  to  a  minimum  depth 
of  10  merers  if  t:..  calculated  depth  for  a  grid  point  is  less  than  30 
meters.  Limited  experience  in  the  application  of  this  feature  of  the 
inesoscai*-  model  indicates  that  computational  stability  is  preserved 
without  significant  loss  of  momentum  over  the  entire  grid  and  with  the 
result  dial  wind  fields  in  the  immediate  vicinity  of  higher  elevations 
under  these  conditions  appear  to  be  reason -.ble. 


The  computerized  algorithm  uses  an  expanding  grid  in  the  boundary 
area  beyond  the  terrain  region  of  interest  to  damp  outgoing  waves  and 
eliminate  incoming  waves  based  on  the  technique  suggested  by  Lavoie  (1972). 
Grid  points  in  the  boundary  region  are  placed  at  distances  of  10,  20,  40, 

80  and  160  kilometers  from  all  four  sides  of  the  interior  grid.  Hie  ter¬ 
rain  heights  at  all  points  in  the  boundary  region  are  set  equal  to  the 
minimum  elevation  within  the  interior  grid.  The  initial  momentum  is  set 
equal  to  the  momentum  used  in  initializing  parameters  in  the  interior 
grid. 


The  form  of  the  f inite-dif f erenced  equations  and  other  details 
on  which  the  mesoscale  wind-field  computer  algorithm  is  based  are  con¬ 
tained  In  the  paper  by  Tingle  and  Bjorklund  (1973). 

2.2  TRAJECTORY/TRANSPORT  ROUTINE 

The  Tra jectory/Transport  Routine  in  the  Applications  Phase  of 
AMS0RB  shown  in  Figure  1-1  has  been  abstracted  from  the  EPAMS  program 
(Dumbauld  and  Bjorklund,  1977).  This  routine  is  designed  to  calculate: 

(1)  Cloud  trajectories  and  mean-layer  winds  and  mixing-layer 
depths  at  fixed  intervals  along  the  trajectories  for  use 
in  the  intermediate-  and  long-range  dispersion  models 
described  in  Section  2.3.1  below 

(2)  The  mean  wind  direction,  mean  wind  speed  and  mixing- 
layer  depth  at  the  location  of  sources  located  near  the 

1 ines-of-sight  for  use  in  the  obscuration-dispersion  models 
described  in  Sections  2.3.2  and  2.3.3  below 

In  both  calculations,  the  Tra jectory/Transport  Routine  uses  the  wind 
vector  and  mixing-layer  depth  values  given  by  the  Hesoscale  Wind-Field 
Model  at  the  solution  grid  points. 
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The  cloud  trajectory  from  a  release  at  any  point  within  the 
mesoscale  model  solution  grid  is  calculated  using  a  simple  Euler  pre¬ 
dictor-corrector  scheme  and  nine-point  interpolation  procedure.  The 
computer  program  uses  wind  vector  data  from  the  Mesoscale  Wind-Field 
Model  solution  at  nine  grid  points  surrounding  the  point  of  release  and, 
through  the  interpolation  procedure,  calculates  the  u-  and  v-components 
of  the  wind  at  the  point  of  release.  If  the  smoke-dispersion  models 
described  in  Section  2.3.2  are  to  be  used  in  the  ensuing  calculations, 
these  components  are  resolved  and  the  mean  wind  direction  and  speed  are 
made  available  to  the  smoke  dispersion  model.  In  calculating  cloud 
trajectories,  the  first  point  on  the  trajectory  is  determined  by  resolving 
the  components  and  placing  a  point  at  a  fixed  distance  interval  along 
the  resultant  vector.  The  fixed  distance  interval  is  set  at  one-tenth 
the  grid  spacing  used  in  the  mesoscale  wind- field  solution  matrix  (in 
the  present  program  a  fixed  interval  of  500  meters  is  used).  Wind- 
velocity  components  at  this  point  are  also  estimated  from  the  nine-point 
interpolation  procedure.  The  new  u-  and  v-components  are  then  averaged 
with  the  previously-determined  components  at  the  starting  point.  The 
final  estimate  of  the  first  point  on  the  trajectory  is  then  calculated 
by  resolving  the  averaged  components  and  moving  a  fixed  interval  along 
the  resultant  vector.  This  final  estimate  of  the  first  point  on  the 
trajectory  is  then  treated  in  the  same  manner  as  the  initial  starting 
point  (source)  and  the  interpolation  and  Iteration  scheme  repeated  to 
find  the  second  point  on  the  trajectory.  This  procedure  is  repeated 
until  the  trajectory  passes  outside  the  grid  area  containing  wind-vector 
data  or  the  number  of  iterations  exceeds  a  predetermined  limit. 

The  Interpolation  method  is  based  on  an  evaluation  of  a  trun¬ 
cated  two-dimensional  Taylor's  expansion  formula  (Carnahan,  Luther  and 
Wilkes,  1969,  p.  430)  given  by 


f (x, y) 


+  hf  +  kf 
*  y 


f 

XX 


+  hkf 

xy 


(2-5) 
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where 


x  -  x, 


y  -  y< 


and  f  is  either  the  u-or  v-component  at  x  ,  y  ,  the  nearest  grid 

*  »  J  ^  J 

point  to  the  point  of  interpolation  (x,y).  The  partial  derivatives 

f  ,  f  ,  f  ,  f  and  f  are  estimated  by  divided  differences, 
x  y  xx  xy  yy 

The  depth  of  the  mixing  layer  at  fixed  intervals  along  the  tra¬ 
jectory  and  at  the  release  point  or  smoke  source  location  is  also  inter¬ 
polated  from  data  provided  by  the  Mesoscale  Wind-Field  Model  using  a 
four-point  bivariate  interpolation  formula  given  by 


U  (x,y)  =  (1  -  p)  (1  -  q)  H  +  p  (1-q)  H 

i,j  i+l,j 


+  q  (1  -  p)  H  +  pq  H 

i,j+l  i+j,j+l 


(2-6) 


where 


(x  -  x ^ )  /  (xi+1  -  xt) 


(y  -  yj)  /  (yj+1  - 


and  H  ,  H 

mi,j  mi+l,J  mi,j+l 


,  H  ,  and  H  are  the  mixing  layer 


i+l.J+l 

depths  at  the  four  surrounding  grid  points. 


2.3 


DISPERSION  MODELS 


The  AMSORB  program  retains  the  Dispersion  Models  Routine  and  the 
MS3  Routine  previously  used  in  the  Application  Phase  of  the  EPAMS  program 
described  by  Dumbauld,  Saterlie  and  Cheney  (1979).  The  volume-source 
dispersion  models  for  instantaneous  and  continuous  sources  and  models  for 
predicting  the  buoyant  rise  of  continuous  stack  emissions  contained 
in  the  Dispersion  Models  Routine  are  used  to  calculate  isopleths  of 
dosage  and  concentration,  with  or  without  cloud  depletion  due  to  decay, 
gravitational  settling  and  precipitation  scavenging,  along  the  cloud 
trajectory  predicted  by  the  Trajectory/Transport  Routine.  The  models 
are  specifically  formulated  to  account  for  changes  in  mean  wind-speed 
and  depth  of  the  surface  mixing-layer  along  the  predicted  trajectory 
while  conserving  mass  continuity.  The  Dispersion  Models  Routine  in 
AMSORB  can  be  selected  by  the  program  user  as  a  "stand-alone"  dispersion 
routine  for  predicting  dispersion  at  intermcdiate-and  long-travel  dis¬ 
tances  from  instantaneous  and  continuous  volume  sources. 

The  MS3  Routine,  added  to  EPAMS  under  Contract  No.  DAEA18-77-C-0060 
with  ASL,  is  primarily  for  use  in  predicting  line-of-sight  integrated  con¬ 
centration  (CL)  and  dosage  (CLID)  downwind  from  quasi-continuous  and 
instantaneous  obscurant  sources.  The  model  has  no  specific  means  of 
accounting  for  changes  in  Che  mean  wind-speed  and  depth  of  the  surface 
mixing  layer  along  a  trajectory,  and  is  thus  intended  for  use  over  shorter 
distances  where  these  meteorological  parameters  can  be  considered  constant 
in  time  and  space. 

The  Battlefield  Environment  Routine  in  the  Applications  Phase 

of  AMSORB  incorporates  features  of  both  the  Dispers'd!  Models  and  MS3 

Routines  to  calculate  CL  and  CLID  along  user  specified  lines  of  sight. 

The  techniques  used  to  calculate  the  cloud  dimensions  0  ,  0  and  o  in 

x  y  2 

the  Dispersion  Models  Routine  are  used  in  the  Battlefield  Environment 
Routine  to  define  cloud  dimensions  along  the  trajectory  from  the  source 
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to  a  distance  of  about  one  mesoscale  model  grid  spacing  from  the  lines-of- 
sight.  The  dimensions  of  the  cloud  at  this  point  become  the  initial 
source  dimensions  for  use  in,  for  example,  the  quasi-continuous  obscurant 
source  model  to  calculate  CL  and  CLID  along  the  lines-of-sight.  Dis¬ 
persion  models  have  been  added  to  the  Battlefield  Environment  Routine  for 
an  obscurant  released  by  traveling  vehicles.  The  dispersion  models  used 
in  the  routines  are  described  in  more  detail  in  the  following  paragraphs. 


2.3.1  Dispersion  Models  Routine 

The  volume-  and  continuous-source  models  in  the  Dispersion 
Models  Routine  of  AMSORB  are  similar  to  models  described  by  Cramer,  et  al. , 
(1972). 


Instantaneous  Volume-Source  Model 


The  concentration  for  an  instantaneous  volume-source  can  be 
expressed  as  the  product  of  four  terms: 


=  {Peak  Term}  {Lateral  Term}  {Vertical  Term}  {Depletion  Term}  (2-7) 


The  Peak  Term  describes  the  concentration  at  the  cloud  centerline  and  is 
given  by  the  expression 


Peak  Term 


(2tt) 3/2  a  o  o 
xii  yn  zn 


(2-8) 


where 

Q  ®  total  source  strength 


0  -  standard  deviation  of  che  concentration  distribution  in 

Xf\  i 

the  alongviud  direction  at  the  ntl>  distance  along  the 
cloud  trajectory 
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a  =  standard  deviation  of  the  concentration  distribution 
jTl  in  the  crosswind  direction  at  the  nth  distance  along 

the  cloud  trajectory 


o 


zn 


standard  deviation  of  the  concentration  distribution 
in  the  vertical  at  the  nch  distance  along  the  cloud 
trajectory 


The  Vertical  Term  describes  the  vertical  distribuiton  of  mate¬ 
rial  at  the  n^  distance  along  the  cloud  trajectory  as  modified  by  re¬ 
flections  of  material  at  the  top  of  the  surface  mixing  layer  and  at  the 
earth's  surface  and  by  gravitational  settling  of  particulates.  The  expres¬ 
sion  for  the.  Vertical  Term  is 


/ 

r—  — 

/  w  Ax  \  ^ 

/2a  H  -  H'+  z  +  V  >  '  — M 

Vertical  _  ^ 

a 

1  [  m>n  S  °n  1 

Term  /  j 

ir  exp 

1 

"  2  \  / 

To{ 

I 

c 

M 

o 

where,  for  convenience  in  writing  the  Vertical  Term,  0^  (zero  to  the 
zeroeth  power)  is  set  equal  to  unity  and 


H1  =  effective  release  height  (see  Equation  (2-21)) 
z  =  calculation  height 

H  =  mixing  layer  depth  from  mesoscale  model  for  n^  dis- 
m.n  ,  , 

tance  along  the  trajectory 

r  =  fraction  of  material  reflected  at  the  surface  (1  for 
complete  reflection  and  0  for  no  reflection) 

V  =  settling  velocity  for  material  in  a  given  size  cate¬ 
gory 

Ax^  =  incremental  distance  along  the  cloud  trajectory 

u  =  mean  wind  speed  in  the  layer  containing  the  cloud  at 
the  nc  '  point  along  the  trajectory  (see  Equation 
(2-18)  below) 


where 


y  -  crosswind  distance  from  the  centerline  of  the  cloud  tra¬ 
jectory 
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The  Depletion  Terra  refers  to  the  loss  of  mate-ial  by  simple 
decay  processes  and  by  precipitation  scavenging.  The  form  of  the  Deple¬ 
tion  Terra  for  each  of  these  processes  is: 


(Decay) 


exp 


Ax 

n 

u 

n 


(2-11) 


(Precipitation 

Scavenging) 


(2-12) 


where 

k  =  decay  coefficient  or  fraction  of  material  lost  per  unit 
time 

A  =  washout  coefficient  or  fraction  of  material  removed  by 
scavenging  per  unit  time 


The  subset  of  equations  describing  the  standard  deviations  of 

Che  alongwind  a  ,  crossvind  o  and  vertical  c  concentration 
xn  yn  2n 

distributions  is: 


o 

xn 


0.6  Au 

n 


4.3  a 

n 


(°x(n-l) 


(2-13) 


where 


i. 


Au 

n 


alongwind  dimension  of  the  cloud 


wind  speed  shear  at  the  n  point  along  the  tra¬ 
jectory 


1? 


(2-14) 


/_ pn 

5Pn  \T*" 


=  wind  speed  at  a  height  of  five  meters  at  the  n 
*  point  along  the  trajectory 


p  •»  wind  power-law  coefficient  at  the  n  point  along 
the  trajectory 


z  =  top  of  the  cloud  at  the  n  point 

*  i  ^ 


U'  +  2.15o  ;  zT  <  H  N 

zn  T ,  n  m ,  n  j 


;  Zn-  >  1' 

i  ,n  —  ra ,  n 


z_  =  base  of  the  cloud  at  the  n  point 
B,n  r 


H*  -  2.15  a  ;  z.,  >  0 

zn  B,n 


;  ZB,n  i  °' 


(2-15) 


(2-16) 


u  (H  -  5)  5P»  (l  ♦  p  ) 
n  \  m,u  /  V  nJ 


H(l+Pn) 

m,n 


,n+Pn) 


(2-17) 


u  »  mean  layer  wind  speed  from  mesoscale  model  at  n 
n  . 

point 


«  mean  wind  speed  in  the  layer  containing  the  cloud 


> ,  n  \  T ,  n _ _ B ,  n  / 

(HO  (S)*’ 


(2- IB) 


f. 


\.14*  - 

:  I 

f"  i 


Ax 


incremental  distance  between  Che  n  points  (e<v>al  to 
500  meters) 


yn 


A, ft 


1/a 


i\%  + 


JLiSzl 
°A,n  / 


a 


(2-19) 


whe  re 


0.  □  standard  deviation  of  the  wind  azimuth  angle  at  the 

A.n  nth  ^>0int 


°A.n  Uo)  (f“)  Uo) 


(2-20) 


m  =  power-law  exponent  for 
ubwve  ground 


variation  ot  d  with  height 


t:  ^ 


H* 


v 


H  >  5| 


(2-21) 


rl„.  rimc!  5  seconds  for  instantaneous 
t  »  source  tunc t ton  time 


sources 


\o 

!• 


0'  ~  ten-minute  s 

A,n 


tandard  deviation  of  the  wind  azimuth  angle 


0  "  °E  n 

cti  *•»“ 


Ax 

U 


♦  (1 Utt 
n  \  Ov 


vl/o 

A 

\  Id  /  J 


(2-22) 


where 


.  *  a  (s/160) 

aE.n  h.n 


(2-25) 
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q 


power-law  exponent  for  variation  of  o  with 
height  above  ground  E‘ 

standard  deviation  of  the  wind  elevation  angle 


The  cloud  expansion  coefficients  a  (lateral)  and  B  (vertical)  in  the 
above  equations  are  set  to  unity  for  instantaneous  sources. 

The  concentration  y  is  calculated  under  the  additional 
n 

restriction  that  y  <  X  , •  If  X  >  X  , »  the  values  of  o  and 
An  —  An-1  An  An-1  yn 

om  are  increased  proportionately  so  that  Xn  =  Xn  ^  and  the  calcu¬ 
lations  continued. 

Dosage  for  an  instantaneous  source  is  obtained  from  the  expres¬ 
sion 


D 

n 


u 

z  ,n 


(2-24) 


where 


u 

z,n 


h  Vi 

mean  wind  speed  at  the  calculation  height  z  and  n 
point  on  the  trajectory 


v,,(!)P”;  *>■) 


f-  /  1  Vu  ,  , 

n,  (  r  i  ;  z  <  1 
^  5 ,  n  '  > '  — 


Pn 


(2-25) 
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The  Continuous-Source  Model  used  in  AMSORB  is  a  modified  version 
of  the  Gaussian  model  for  continuous  sources  described  by  Pasquill 
(1962)  in  which  the  concentration  is  also  expressed  as  the  product  of 
four  terms: 

X  =  {Peak  TermH Lateral  Term} {Vertical  Term} {Depletion  Term}  (2-26) 
c,  n 


The  Peak  Term  of  the  continuous-source  concentration  model  is  given  by  the 
expression 


Peak  Term 


2tt  0  o'  u 
zn  yn  n 


(2-27) 


where 


source  strength  expressed  as  a  rate 


standard  deviation  of  the  concentration  distribution  in 
the  crosswind  direction  for  a  continuous  source  at  the 
nc*1  distance  along  the  plume  trajectory 


and  u^  and  o  are  respectively  defined  by  Equations  (2-18)  and  (2-22). 
The  Vertical,  Depletion  and  Lateial  Terms  (with  o  replaced  by  0^ 
where  appropriate)  are  identical  to  che  similar  terms  described  for  the 
instantaneous  volume  source.  The  subset  oi  equations  describing  0^ 
is: 


o!  Ax 

AC  ,  tt  n 


o'  ' 

y.n-l 

o' 

,  Ac,»  - 


(2-28) 
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where 


a !  =>  a . 

Ac,n  A,n 


IT 

180 


(2-29) 


and  Che  remaining  parameters  have  been  defined  previously.  The  default 
values  o f  u  and  fj  used  in  Che  program  for  continuous  sources  are  a  =  0.9 
and  i;,  =  1 . 


Plume-Rise  Models  for  Continuous 
Stack  Emissions 

The  effective  height  H  of  a  buoyant  plume  is  given  by  the  sum  of 
the  physical  stack  height  h  and  the  buoyant  rise  Ah.  For  a  neutral  or 
unstable  atmosphere,  the  effective  stack  height  is  given  by  (after  Briggs, 
1971;  1972) 


II 


h  +  t 


)  l_r> .  98 
')  ufh'j 


(2-30) 


where 


r  a  the  inner  radius  of  the  stack  at  the  exit  (m) 
s 

u  -  the  stack  exit  velocity  (m  s~^) 

T  =  the  stack  exit  temperature  (°K)  f  >r  the  continuous  source 

T  «  ambient  air  temperature  (°K) 
y 

f  »  empirical  correction  factor 
h  -  actual  stack  height  (a) 
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mean  wind  speed  at  stack  height  (m  s 


u{h}  “ 


u 


5  ,n=  1 


(2-31) 


Vhe  factor  f,  which  limits  the  plume  rise  as  u(h}  approaches  the  stack 
exit  velocity  w,  is  defined  by 


f  •“  S 


1  ;  u{h)  <  w/1.5 


^3w  -  3u{h}^ 


w/1.5  <-  u{ h}  <  w  V 


;  u{h}  >  w 


(2-32) 


Equation  (2-30)  is  used  when  the  net  radiation  index  equals 
or  exceeds  2,  uth}  equals  or  exceeds  5  meters  per  second,  or  the 
potential  temperature  lapse  rate  $  is  equal  to  or  less  than  0.  When 
these  conditions  are  not  met,  the  program  uses  the  corresponding  Briggs 
(1972)  rise  formula  for  a  stable  atmosphere  given  by 


H  -  < 


vr2X 

h  +  l-905 


1  -  cos 


10  Sl/2  h 


>  (2-33) 


(1/3 


nJifei  .  . 

- lTf  -  h 

10  s' 


uhere 


S  =  stability  parameter 
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(2-34) 


9.8  ; 

—  V 

s 


(*.  =  height-weighted  mean  potential  temperature  lapse 

rate  over  an  approximately  200-meter  interval  above 
h  ’ 


z  =  h  +  200 

V  4-  AZ 

Z—J .  z  z 

v.  h _ 

z  =  h  +  200 


(2-35) 
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T  -  T 

A,z-H  A,z 

AZ 


(2-36) 


AZ 


z 


'z+1 


Z 

z 


(2-37) 


A,  st 


T. 


Z.i+l 


potential  temperature  at  height  Z^ 
potential  temperature  at  height 


2.1.2  MS 3  Rout ine 

The  obscurant  dispersion  models  in  the  MS3  Routine  calculate 
CL  and  CL1D  downwind  from  quasi-eonc iuuous  and  instantaneous  sources. 

tjuas  1-Cantinuous  Qba curant  Source  Model 

According  to  experiments  conducted  at  Dugway  Proving  Ground  (Salomon, 
ct  al.  1978  and  Petersen,  1978),  the  amount  o£  aerosolized  smoke,  H^it). 
at  any  given  time  Jrorn  quasi-eont inuous  stmjke  tau‘<itions  is  given  by  the 
expression 
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M  ♦  MYF  •  YF 
o 


(2-38) 


w 


M  {t } 
x 


whore 

M 

o 

MYF 


YF 


t 


B 


t 

A.B.C.D 


total  mass  in  grams  before  ignition 
munition  yield  fraction 

yield  factor  accounting  for  the  effects  of  moisture 
in  the  air 


total  burn  time  in  minutes 


time  in  minutes 


experimentally-derived  coef f lcients 


The  amount  of  smoke  produced  by  the  munition  as  a  function  of  time  is, 
therefore , 

QUI  ■  «„  (c| 


or 

3t.d  40c1  \ 

»  V  s’  <) 

This  expression  is  used  to  character ice  the  smoke  munition  source  strength 
in  the  dispersion  model.  It  should  he  noted  that  Equation  f2-J9)  is  an 
extremely  important  factor  in  determining  the  shape  and  magnitude  of  the 
CL  time  profile  at  short  distances  downwind  from  a  quasi-continuous 


q{t)  “  M  •  MYF  •  YF 


source.  The  smoke  dispersion  model  performance  has  been  compared  with 
smoke  munition  trial  data  by  Carter,  Dumbauld  and  Rafferty  (1979). 


i  . 


t 

i 


Quasi-Continuous  Source  Obscui <nt  Dispersion  Model 

The  concentration  at  time  t  in  seconds  at  a  point  x,y,z  down¬ 
wind  from  a  quasi-continuous  volume  source  emitting  smoke  over  a  time 
t  is  given  by  the  expression 
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where 


Q { c }  3  source  emission  rate  (see  Equation  (2-39)) 

0  =  standard  deviation  of  the  alongvind  concentration  dis¬ 

tribution  (n) 

o  =  standard  deviation  of  the  crosswind  concentration 
^  distribution  (m) 


0  »  standard  deviation  of  the  vertical  concentration 

2 

distribution  (tn) 


u  =  mean  cloud  transport  speed  (ms  ) 

H{t}  =  effective  height  of  the  cloud  centroid  at  time  t  (m) 

H  =  depth  of  the  surface  mixing  layer  (in) 

If  we  perforin  the  Integration  over  the  variable  r,  the  solution  for  time 
t  is 
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where 


=  the  lateral  and  vertical  terns  following  the  first  set 
of  closed  brackets 


and 
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The  standard  deviation  of  the  alongwind  concentration  distribution 
used  in  the  smoke  dispersion  model  is 
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standard  deviation  of  tlto  alongwind  concentration 
distribution  at  the  source 


vertical  wind-speed  shear  in  the  layer  containing 
the  cloud 


(2-53) 


Mean  wind  speed  at  the  reference  height  of  5  meters 
(see  Equation  (2-17)) 
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effective  lower  bound  of  the  cloud 
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average  standard  deviation  of  the  vertical  concentration 
distribution  along  the  1 ine-o£-sight 
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x^,X2  11  downwind  distance  coordinates  of  the  line-o£-sight 
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The  standard  deviation  of  the  crossulud  concent  rat  ion  distribution  is 
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where 
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»  standard  deviation  of  the  azimuth  wind  angle  in 
radians  measured  over  the  time  t 
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standard  deviation  of  the  azimuth  wind  angle  in 
radians  measured  over  the  reference  time  t 
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expansion  occurs  downwind  from  a  virtual  point  source 
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change  of  wind  direction  in  degrees  with  height  in 
the  surface  mixing  layer 


The  standard  deviation  of  the  vertical  concentration  distribution  used 
in  the  smoke  dispersion  model  is 
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where 


o',  ■  standard  deviation  o£  the  elevation  wind  angle  in  radians 
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x  »  distance  over  which  rectilinear  vertical  cloud  expansion 
r  2 

occurs  downwind  from  a  virtual  point  source 
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sjatidafij  deviation  of  the  Vertical  concent  tat  iott  dlsttl- 

but  ion  at  a  reference  distance  x„  downwind  from  the 

Hi 

source 


K4u.it ion  yields  concent rav ion  as  a  function  ul  time  at  3 

specific  (joint  (x ,  v,  z ) .  in  the  computer  program,  the  1  inc-of -sight 
integrated  concentration  l*  obtained  by  numerically  integrating  Equation 


(2-41^  over  the  line-of-sight  using  a  Gauss-Legendre  technique  (Abramowit: 
and  Stegun,  1964,  Equation  25.4.30,  page  887)  with  24  arbitrary  points 
in  the  interval  containing  the  cloud. 

Tnstantaneous  source  Obscurant  Model 


The  line-of-sight  integrated  concentration  for  r.n  instantaneous 
source  can  be  obtained  from  Equation  (2-40)  by  setting  Q { X }  equal  to  the 
total  source  strength  Q  and  integrating  over  the  line-of-sight  rather 
than  time.  The  resuJ  t  of  integrating  between  the  end-points  (x-^ ,  y^ ,  z^) 
and  (x^,  y^,  of  the  line-of-sight  is  given  by  the  expression 
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(Equation  (2-64)  is  continued) 
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(Equation  (2-74)  is  continued) 
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The  remaining  parameters  in  liquations  (2-6A)  through  (2-82)  have  been 
defined  under  the  discussion  of  quasi-continuous  sources  above,  except 
that  x  for  ail  instantaneous  sources  is  defined  by 

x  =  u  t  (2-83) 


Cloud -Kise_  Models  For  Buoyant  Instantaneous 
and  Qua?;  i  -Cont  inuous  Obscurant  Sources 

The  MS3  Routine  incorporates  cloud-rise  models  for  instan¬ 
taneous  and  quasi-continuous  smoke  sources  which  release  appreciable 
heat  during  the  smoke  emission  phase.  The  models  given  below  are 


primarily  based  on  material  contained  in  a  preprint  of  a  paper  by  G.  A. 
Briggs  (1971)  presented  at  the  Second  International  Clean  Air  Congress. 

The  rise  models  given  by  Briggs  for  stable  atmospheres  are  specified. 
Experience  in  using  these  models  for  predicting  the  rise  of  ground 
clouds  from  rocket  launches  and  of  plumes  from  stacks  indicates  the 
stable  formulas  predict  nearly  the  same  rise  as  that  given  by  the 
Briggs  models  for  adiabatic  atmospheres,  if  the  vertical  gradient  of  po¬ 
tential  temperature  is  arbitrarily  set  to  a  small  value  (3.322  x  10  ^ 
degrees  meter  S  when  the  atmospheric  lapse  rate  is  adiabatic  or  unstable. 

The  effective  height  of  the  cloud  centroid  downwind  from  an 
instantaneous  source  is  given  by 
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where 


It  -  source  height  (m) 
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g  »  acceleration  due  to  gravity  (9.8  m  s”“) 
Qj  »  heat  released  (cal  ) 


;  .  'U  *  ■  - 

-i — 


39 


specific  heat  of  air  at  constant  pressure  (0.24  cal  g 

Vl) 

-3 

density  of  the  ambient  air  (g  m  ) 
ambient  air  temperature  (°K) 

entrainment  coefficient  for  an  instantaneous  source  ^~0.6) 
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height  weighted  mean  potential  temperature  lapse  rate  over 
a  height  interval  of  about  50  meters  above  the  surface 
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potential  temperature  at  height  Zr 


potential  temperature  at  height  Z 


The  effective  height  of  the  cloud  centroid  downwind  from  a  quasi- 
continuous  source  is  given  by  the  expression 
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where 
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=  effective  rate  of  heat  release  in  calories  per  second 


l  = 


entrainment  coefficient  for  a  continuous  source  ( — 0 . 66 ) 
time  after  munition  is  ignited 


Equation  (2-84)  for  instantaneous  sources  and  Equation  (2-90)  for  quasi- 
continuous  sources  are  automatically  used  to  calculate  the  effective 
height  H{ t }  when  values  of  either  Qj  or  are  specified  in  the 

lookup  dar«  tables  for  specific  obscurant  sources. 
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2.3.3 


Battlefield  Environment  Routine  Dispersion  Models 


As  noted  in  the  Introduction  to  Section  2.3,  the  Battlefield 
Environment  Routine  incorporates  the  same  models  described  in  Sections 
2.3.1  and  2.3.2  above.  When  the  Battlefield  Environment  Routine  is 
selected  by  the  user,  the  program  automatically  selects  the  proper  model 
for  use  in  calculating  CL  and  CLID  along  the  lines-of-sight  specified  by 
the  user  depending  on  the  source  category  and  the  distance  of  the  source 
from  the  lines-of-sight.  For  example,  suppose  the  user  has  identified  an 
explosive  munition  source  (see  Section  3.2  below  for  a  description  of 
explosive  munition  sources  handled  by  the  program)  functioned  at  7  kilo¬ 
meters  from  the  line-of-sight .  Further  assume  that  the  program  has  used 
the  Mesoscale  Wind-Field  Model  and  Trajectory/Transport  Routine  to  deter¬ 
mine  that  the  source  contributes  to  the  obscurant  concentration  along  the 
line-of-sight  at  the  time  of  interest  and  the  mesoscale  grid  solution  is 
based  on  a  2.5-kilometer  spacing.  The  program  would  then  use  the  tech¬ 
niques  described  in  Section  2.3.1  for  calculating  the  concentration  and 
cloud  dimensions  along  the  cloud  trajectory  until  the  cloud  trajectory 
intersects  the  boundaries  of  the  grid  square  containing  the  intersection 
of  the  cloud  trajectory  and  the  line-of-sight.  The  cloud  dimensions  0 

xn 

O  and  J  determined  to  exist  at  this  boundary  from  application  of  the 
yn  zn  '  r 

Instantaneous  Volume  Source  Model  described  in  Section  2.3.1  and  the 
initial  source  strength  Q  are  then  used  to  define  the  initial  conditions 
for  the  application  of  the  Instantaneous  Source  Obscurant  Model  described 
in  Section  2.3.2  above.  The  Instantaneous  Source  Obscurant  Model  uses  the 
meteorological  parameters  for  the  grid  square  containing  the  line-of-sight 
in  the  calculation  of  CL  and  CLID  along  the  line-of-sight.  The  probability 
of  detecting  a  target  is  then  calculated,  if  the  user  desires,  using  the 
Obscuration  Model  Routine  described  in  Section  2. A  below.  If  the  source 
had  been  determined  by  the  program  to  be  within  the  same  grid  square  as 
the  line  of  sight,  then  the  program  would  have  used  the  source  param¬ 
eters  from  the  explosive  munition  source  character istics  model  and  the 
meteorological  inputs  for  the  grid  square  directly  in  the  Instantaneous 
Source  Obscurant  Model. 
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A  similar  procedure  is  followed  when  the  user  defined  source 
category  is  for  a  quasi-continuous  obscurant  source,  except  the  Battlefield 
Environment  Routine  uses  the  Continuous-Source  Model  described  in  Section 
2.3.1  and  the  Quasi-Continuous  Source  Obscurant  Model  described  in 
Section  2.3.2.  It  should  be  noted  that  the  primary  purpose  for  using 
the  Dispersion  Models  Routine  described  in  Section  2.3.1  in  the  Battlefield 
Environment  Routine  is  to  simply  account  for  the  effects  of  changes  in 
wind  speed  and  depth  of  the  surface  mixing  layer  along  cloud  trajectories 
that  exceed  the  mesoscale  grid  spacing  while  preserving  mass  continuity. 

When  the  trajectory  between  the  source  and  line  of  sight  is  contained 
within  a  grid  square  or  enters  the  grid  square  containing  the  line  of 
sight,  the  resolution  of  the  mesoscale  model  is  such  that  wind  speed  and 
depth  of  the  surface  mixing  layer  are  assumed  constant  along  the  trajectory. 


Mobile  Line  Source  Dispersion  Model 


A  Mobile  Line  Source  Model  routine  has  been  added  to  the  battle¬ 
field  environment  routine,  which  is  not  available  for  use  in  the  MS3  Routine. 
The  line  source  geometry  is  shown  in  Figure  2-1.  The  vehicle  is  assumed  to 
move  from  0  to  point  A  at  a  constant  speed  V^,.  The  concentration  at  point 
r{a,3,z}  from  an  instantaneous  source  generated  at  point  S  by  the  moving 
vehicle  Is  given  by 
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(Equation  (2-92)  is  continued) 
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(Continuation) 


i  exp  - 
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where 


Q  =  constant  emission  rate  (g  m  ) 


V  =  vehicle  speed  (m  s  ) 


dt'  =  instant  the  source  is  generated  at  point  S 


t  *  time  concentration  is  calculated 


t'  =  travel  time  of  the  vehicle 


and  the  parameters  H  ,  H,  2,  V  ,  o  ,  o  and  0  have  their  usual  meaning.  If 

w  a  x  y  z 

we  start  the  vehicle  at  times  t'  equal  0  and  t  equal  0,  the  concentration 
at  point  r{u.h,zi  from  all  the  instantaneous  point  sources  comprising 
Che  line  source  generated  by  the  vehicle  traveling  over  a  time  t^  can  be 
calculated  by  integrating  Equation  (2-92)  over  the  interval  t^.  Note  that 
the  limit  of  integration  t^  must  he  equal  to  t  for  times  t  less  than  or 
equal  to  t ^  and  equal  to  for  times  t  greater  chan  t^.  The  integration 
can  be  analytically  performed  only  if  the  standard  deviations  0^  and  oa 
have  a  similar  luttciional  dependence  on  the  time  t'.  We  have  performed  the 
integration  under  the  assumption  chat: 
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x'  =  acosO  +  $sin9 


(2-98) 


and  where  T"  is  determined  such  that  0^  given  by  Equation  (2-95)  will  approxi 
mate  the  value  that  would  have  been  obtained  if  Equation  (2-50)  had  been  used 
to  calculate  o^.  The  results  of  integrating  Equation  (2-92)  are  given  by 
the  expression 
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(Equation  (2-99)  is  continued) 
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(VT  rfinQ  -  u)/(/2  T"  u) 
Q  V  k  T’7(tt3/2*^  u) 


(2-119) 

(2-120) 


The  purpose  of  the  parameter  T",  as  noted  above,  is  to  approximate 
Equation  (2-50)  using  Equation  (2-95)  in  a  manner  that,  after  the  integration 
of  Equation  (2-92)  over  the  time  t'  is  performed,  yields  the  average  value 
°f  °x  from  all  the  instantaneous  sources  affecting  the  concentration  at  the 
point  r{a,B,z).  To  obtain  the  appropriate  value  of  T",  we  first  determine 
the  portion  of  the  line  source  contributing  to  the  concentration  at  the  point 
r  at  time  t.  This  is  accomplished  by  calculating  the  origin  of  clouds  from 
instantaneous  sources  comprising  the  line  source  which  lie  within  approximately 
5  in  any  direction  from  the  point  r  at  time  t.  We  then  calculate  a  value 
of  T"  by  setting  Equations  (2-50)  and  (2-95)  equal  to  one  another  and  solving 
for  T",  resulting  in  the  following  expression  for  T" : 


b1  =  Jl  T"  o^h»2.5}[ut  +  xv] 

/2  T"  o^{t=2.  5} j^u(t-t^)  +  xv J ;  t  >  j 

b2  =  ( 

Jl  T"  O' {t=2. 5}  x  :  t  <  t,  I 

A  v  ~  1  J 


T"  =  k 


11/2 

x  x2+(o  /k')2  +  (o  /k')2  £n  x  +lx2+(o  /k')2] 

l  o _ J _ o _ L  o  _ j. 


2  2 

X2  "  X1  +  2xv  ^x2_xl^ 


~1  }  (2-121) 


x=X2 

where  the  symbolic  notation  means  that  the  numerator  is  evaluated  at 

x  equal  x^  and  X2.  x~xi 
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The  values  of  x^  and  are  determined  from  the  coordinate  trans 


formation  matrix 


x  j  \  /cos (£-0)  “  sin(tj-O) 
0  /  \  sin(£-0)  cos(^-0) 


\  /  cos(£-0)  -  sin(£-0) 

0  /  \sin(C-6)  cos(^-0) 


where 


X1  =  ^ai^  ~  R)/cos^(tan  ^(5a^{l=2.5}')) 


x2  =  ^ai2  +  R)/cos2(tan  1(5o^{t=2.5})) 

R  =  a^2  -  (a!2  +  B^2)  cos2(tan  1(5o^{t=2.5})) 


f'a.\  /cos(0-£)  -  sin(0-£) 


Bi/  \sin(0-r)  cos  (0-£) 


-  a  -  VT  t  cos  (90  -  0) 


B.  =  0  -  V  t  sin(90  -  6) 


£  =  tan  1  (u  sin0)/(u  cos0  -  VT) 


0. 6  Au 


4.3  u 


The  angle  £  in  Equation  (2-129)  is  the  angle  between  the  mean  wind  direction 
and  the  locus  of  points  of  the  instantaneous  clouds  generated  by  the  moving 
vehicle  prior  to  the  calculation  time  t.  Inspection  of  Equation  (2-129) 
shows  the  angle  is  dependent  on  the  wind  speed  u  and,  as  shown  by  Equation 
(2-57),  the  wind  speed  in  turn  is  dependent  on  cloud  dimensions.  For  this 
reason,  the  program  performs  an  iteration  to  define  the  appropriate  value 
of  u  and  the  angle  £  for  the  portion  of  the  line  source  generated  cloud 
which  may  be  in  the  vicinity  of  the  calculation  point  at  time  t.  This 
value  of  u  determined  in  the  iteration  is  then  used  in  the  dispersion  model 
(Equation  (2-99)  and  following  equations).  The  parameters  u  and  Au  are, 
respectively,  the  values  of  the  mean  wind  speed  and  wind  speed  sheer  re¬ 
quired  to  equate  Equations  (2-50)  and  (2-95).  When  the  cloud  is  fully 
mixed  in  the  surface  mixing  layer, 


H 

x,  >  x*  =  - 55 -  ~  x  ,  (2-131) 

2.15  a-E  v 


the  program  calculates  u  and  Au  from  the  expressions: 


u 


(l+p)  Z  p  (H  -2) 
R  m 


(h^  -  21+P) 


(2-132) 


Au  = 


(2-133) 
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However,  when  x 


1 


< 


x*, 


they  are  evaluated  from  che  following  expressions 


UR 


(2.15  o')p 


(x2-x1)(l+p)  z  P  )  •  (1+p) 


X**  + 


1+p 


X 

X.  +  ~ 

1  °E 


1+p 


(x2~x**) 


(H  -2) 
m 


H  1+P  -  21+P 


AS  = 


(VX1>  ZRP 


(2.15  o’)P 
E 


(i+p) 


X**  +  ~ 

o' 

E 


1+p 


1+p 


+  (x2  -  x*»)  HaP  -  (x2  -  Xi)  2P  j 


nere 


**  ;  x*  >  x. 


X**  = 


x2  ;  X  <  X*  <  X 


(2 


2  ' 
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I 


OBSCURATION  MODEL  ROUTINE 


2.4 


2.4.1  Obscuration  Model 

The  Obscuration  Model  of  AMSORB  is  a  modified  version  of  a  model 
developed  by  Gomez  and  Duncan  (1978)  for  use  in  the  ASL  Smoke  Obscuration 
Model  (AS1.S0M) .  The  ASLSOM  routines  were  based  on  an  obscuration  model 
contained  in  the  SOMT  obscuration  program  (Johnson,  elt  al .  ,  1972)  developed 
for  the  JTCG/ME.  The  Obscuration  Model  calculates  target,  background  bright¬ 
ness  and  cloud  brightness,  light  transmission  through  an  obscuring  cloud, 
scattering  of  light  by  aerosols,  and  target-background  contrasts.  When 
the  model  routine  is  used  in  conjunction  with  dispersion-model  routines 
described  in  Section  2.3.2  and  2.3.3,  the  primary  output  is  the  probability 
of  detecting  a  target.  The  probability  of  detection  depends  on  the  amount 
and  type  of  obscurant  present,  general  lighting  conditions,  relative  con¬ 
trast  of  the  target  and  background,  and  the  amount  and  type  of  natural 
obscurants  which  may  affect  detection. 

The  ASLSOM  visibility  model  contains  options  for  calculating  the 
probability  of  detection  for  visible,  infrared,  laser  and  electro-optical 
(EO)  sensors.  These  options  and  the  option  to  consider  the  effects  of 
adverse  weather  (rain,  fog,  etc.)  contained  in  ASLSOM  have  been  retained 
in  the  Obscuration  Model  used  in  AMSORB,  although  only  the  option  for 
calculating  the  probability  of  detection  in  the  visible  region  of  the 
spectrum  has  been  verified.  The  Obscuration  Model  and  the  modifications 
made  in  the  ALSOM  smoke  visibility-model  routine  for  application  in  AMSORB 
are  described  in  Section  2.4.2  below. 

The  probability  of  detecting  a  target  is  calculated  from  the 
relationship  (Johnson,  £t  a_l.  ,  1972) 


PD  = 


0.5 


jerf  [l.462  ((C/Cfi)  -  l)  +  1 J 


(2-137) 
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where 


C  =  contrast 

2M 

2  -  M 


(2-138 ) 


M  =  modulation  contrast 

2(W 

=  (2-139) 


B  =  target  brightness  as  seen  by  the  observer 


=  background  brightness  as  seen  by  the  observer 


C  =  threshold  constrast  for  a  30  percent  probability  of 
detection  (visible  spectrum) 


Values  of  determined  by  Blackwell  (1946)  for  various  target  illuminations 
and  sizes  are  used  in  the  AMSORB  program  (see  Section  4.4).  Equation  (2-137) 
shows  that  the  constrast  C  must  be  greater  than  the  threshold  constrast  C 

B 

to  have  a  probability  of  detection  of  greater  than  50  percent. 

The  target  brightness  seen  by  the  observer  is 


B 


t 


+  B 


(2-140) 


where 

B^ 


brightness  at  the  target  location  of  a  target  with 
reflectivity  R 
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transmission  factor  from  the  target  to  observer 


B  =  cloud  brightness  with  respect  to  the  target  due  to  light 
scattered  into  the  line-of-sight  between  the  target  and 
observer 


The  subset  of  equations  defining  B  and  t  is 

o 


B  t.  cos  (180 
°i  1 


(2-141) 


where  the  summation  extends  over  only  those  light  sources  where  is  greater 
than  or  equal  to  90°,  since  light  sources  behind  the  target  do  not  contribute 
to  its  brightness,  and  where 


B  =  brightness  of  the  itn  light  source  (sky,  sun,  terrain) 

°i 

=  scattering  angle  for  the  i^  light  source  (see  Figure  4-5 
in  Section  4.4) 


T 


i 


fraction  of  light  reaching  the  target  from  the  i*^  light 
source 


Vyi*zi 


-  a„{X} 


/ 


Vyt' 


X^U.y.z.c}  d£ 


(2-142) 


(2-143) 
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“tu} 


attenuation  coefficient  for  light  of  wavelength  \  and 
the  obscurant  from  the  2,^  source 


concentration  at  time  t  at  the  point  (x,y,z)  given  by 
Equation  (2-41)  for  the  obscurant  from  the  source 


In  Equation  (2-142),  for  example,  £=1  might  be  a  smoke  source,  1=2  might  be 
a  dust  source  and  so  forth  for  the  total  of  P  different  sources.  The  quantity 
d.d  in  Equation  (2-142)  is  oriented  along  the  line-of-sight  between  the  tar¬ 
get  at  (xt,yt,zt)  and  the  light  source  at  (x^.y^.z^).  Since  each  light 
source  or  sky  segment  is  an  infinite  distance  from  the  target,  the  position 
^xi’^i,zi^  sPecifiecl  by  multiplying  the  direction  vectors  to  the  light 
source  by  a  large  number  (Johnson,  £t  al . ,  1972).  In  Equation  (2-140)  the 
fraction  of  light  transmitted  from  the  target  to  the  observer  t  is  given 
by: 


Tt  =  jj  |  exp  (-  ctjUJCL^) 
l 


(2-144) 


where 


C[.£  =  line-of-sight  integrated  concentration  between  the  target 

and  observer  for  the  obscurant  from  Che  source 


x^U.y.z.t)  di 


(2-145) 


where  the  observer  is  located  at  the  coordinate  (x^.y^.z^  and  is  ori¬ 
ented  along  the  line-of-sight  between  the  target  and  observer. 


5b 


The  cloud  brightness  is  defined  by  the  expression 


j**l  1=1  I  i=l 


*  Eh.  E  cnE(8<.i  w  « 


(2-146) 


where 


fraction  of  light  from  the  j  incremental  cloud  volume 
along  the  line-of-sight  reach! :0  the  observer 


;ryr2j 


I  1  exp  ~  J  Xv;£{x,y,Z,t^ 


x  ,y  ,z 
o  Jo  o 


(2-147) 


C  -  concentration  due  to  the  obscurant  source  in  the 

cloud  element,  assumed  constant  within  the  element  and  eoual 
to  Xv;lUryrVt> 

®  fraction  of  light  reaching  the  j^1  cloud  element  from  the  I** 
light  source 


exp  - 


^(A)  J 


(2-148) 


W'j 


fi£  *  fraction  per  unit  concentration  (due  to  the  Vn  obscurant)  of 

the  light  intensify  which  is  scattered  through  an  angle  $  into 
the  observer's  line  of  sight  * 
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F^.X) 


{2-149) 


k  =  propagation  constant 

=  2tt/X 


and  the  scattering  function  (F^l^.A})  is  a  complex  function  of  the  type  of 
particle,  the  particle-si2e  distribution,  the  refractive  index  of  the  aero¬ 
sols,  the  wavelength  of  light  and  the  scattering  angle.  The  computer  pro¬ 
gram  uses  tabular  values  of  this  function  which  have  been  precalculated 
for  several  types  of  smokes  using  a  Hie  scattering  program  (see  Section 
4.4).  No  Hie  fractions  for  obscurants  other  than  UP  or  HC  smoke  have  as 
yet  been  calculated  and  thus  only  these  values  are  currently  used  in  the 
program  regardless  of  the  obscurant  specified.  When  more  than  one  type 
of  obscurant  is  present,  the  summation  over  N  light  sources  and  M  cloud 
elements  indicated  in  liquation  (2-145)  is  performed  for  each  type  of 
obscurant  and  summed  over  the  P  sources  to  obtain  the  total  cloud  brightness. 

The  brightness  of  a  background  of  reflectivity  is  determined 
in  the  same  manner  as  the  target  brightness  in  Equations  (2-140)  through 
(2-148)  by  replacing  the  end  point  of  integration  (x  ,y  ,^t)  with  the  end 
point  (*b.yb.V  a,u*  Rt  If  the  target  and  background  are  coincident, 

the  only  difference  in  brightness  between  them  is  due  to  the  difference  in 
their  reflectivities. 


Detailed  explanations  of  the  theoretical  background  of  the 
obscuration  model  described  above  are  given  by  Johnson,  et.  al^.  (  1972)  and 
by  Domex  and  Duncan  (1978). 
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2.4.2  Revisions  of  the  ASLSOM  Smoke  Visibility  Computer 
Code 

Revisions  in  the  ASLSOM  sm-ke  visibility  computer  code  were  neces¬ 
sary  to  make  the  code  compatible  with  the  concentration  models  described 
in  Sections  2.3.2  and  2.3.3,  with  AMSORB,  to  improve  the  computational 
efficiency  of  the  combined  programs  and  to  allow  consideration  ot  multiple 
obscurants.  The  significant  revisions  made  in  the  code  were: 

•  Modification  of  the  input  routines 

•  Automation  of  the  threshold  contrast  calculation 

•  Addition  of  a  test  to  determine  if  any  obscurant  is  re¬ 
quired  to  obscure  the  tat  get 

•  Replacement  of  the  routine  determining  whether  and  when 
the  obscurant  cloud  intersects  the  line-of-sight 

•  Addition  of  a  routine  to  divide  the  line-of-sight  into 
a  specified  number  of  cloud  elements 

•  Addition  of  coding  to  allow  simultaneous  calculations  for 
different  types  of  obscurants 

•  Modification  of  the  sky-partitioning  routine 

•  Modification  of  the  test  for  large  obscurant  concentrations 

«  Modification  of  the  test  for  small  obscurant  concentrations 

•  Replacement  of  the  output  routine 
Each  of  these  changes  is  discussed  below. 


Modification s  of  tl ic  Input  Routine 


The  input  routines  of  the  ASLSOM  visibility  program  were  modified 
to  automate  the  probability-of-detection  calculation.  In  most  cases, 
the  data  required  to  perform  the  calculations  are  included  in  the  program 
in  the  form  of  "look-up"  tables,  are  calculated  under  various  assumptions, 
or  are  supplied  to  the  obscuration  model  from  other  portions  of  the  program. 
Details  concerning  the  tables,  and  input  calculation  procedures  are  given 
in  Section  A. 4  below.  For  calculating  the  probability  of  detecting  tar¬ 
gets  in  the  visual  spectrum  (visible  scenarios),  the  only  input  required 
is  the  position  of  the  target  and  observer  and  a  code  specifying  target 
and  background  characteristics.  The  target  is  defined  to  be  the  object 
being  viewed  by  the  observer.  For  infrared,  laser  and  other  wavelengths 
outside  the  visible  spectrum,  the  current  version  of  the  program  requires 
that  all  inputs  needed  by  the  original  code  be  supplied. 

AMSOKB  uses  the  tabular  values  of  threshold  contrast  described 
in  Section  4.4  to  calculate  the  probability  of  detection,  where  the 
ASLSOM  and  S0M1  codes  required  the  user  to  input  appropriate  values.  The 
use  of  "look-up"  tables  Incorporated  within  the  program  structure  has  the 
advantage  in  that  the  threshold  contrast  can  be  determined  after  the  target 
and  background  brightness  calculations  have  been  performed.  This  allows 
the  program  to  choose  a  value  based  on  whether  the  target  actually  appears 
brighter  or  darker  than  the  background,  rather  than  basing  the  value  only 
on  a  comparison  of  the  two  reflectivities.  A  target  that  is  brighter 
than  its  background  under  conditions  when  no  obscurant  is  present  may 
appear  darker  than  its  background  due  to  obscurant  above  the  target  attenu¬ 
ating  the  light  which  illuminates  it.  For  a  target  and  background  which 
are  coincident  ,  this  is  not  a  problem  ^nd  a  comparison  of  reflectivities 
can  be  used. 


bO 


Addition  of  a  Test  To  Determine  If  Any  Obscurant 
is  Required  to  Obscure  the  Target 


In  some  cases  the  target  and  background  reflectivities  may  be 
nearly  identical  so  that,  under  the  assumptions  made  in  the  Obscuration 
Model  that  contrast  is  the  only  factor  affecting  detection,  no  obscurant 
is  required  to  obscure  the  target.  The  test  incorporated  in  the  modified 
code  compares  the  reflectivities  of  the  target  R^  and  background  R^  and 
defines  the  contrast  as 


C 


(2-150) 


The  value  of  C  is  compared  to  the  threshold  contrast  from  Blackwell's 
(1946)  data  for  a  given  illumination  and  target  size,  using  the  tabular 
values  described  in  Section  4.4.  If  C  is  less  than  C^,  the  probability 
of  detection  is  calculated  from  Kquation  (2—137)  with  the  value  of  C  given 
by  Kquation  (2-150).  Kor  a  probability  of  detection  less  than  or  equal  to 
5  percent,  the  calculations  are  stopped  and  a  solution  is  printed  stating 
that  no  obscurant  is  required  to  obscure  the  target.  If  the  probability 
of  detection  is  grea  t*-  than  5  percent,  the  complete  calculations  using 
Equations  (2-118)  thro,  '2-149)  arc  performed  to  find  the  actual  value  of 
the  contrast  C. 

Replacement  of  the  Routine  Determining 
Whether  and  When  the  Obscurant  Cloud 
Intersects  the  Linc-ef-Slght 


1  he  modified  obscuration  routine  in  AKSORB  docs  not  require  the 
user  to  specify  times  uheu  the  probability-of-detcetion  or  line-of-sight 
concentrations  are  to  be  calculated  as  in  the  ASLSOM  code.  Instead,  the 
program  automatically  determines  whether  the  obscurant  cloud  Intersects 
the  linc-of -sight  and  the  time  interval  required  for  the  obscurant 
cloud  to  pass  through  the  iiue-yf-sight .  The  cloud  arrival  time  at  the 


ftl 


where 

x  .  =  minimum  alongwind  separation  distance  between  the  source 

min  ,  ,  .  .  .  , 

and  any  point  on  the  line-of-sight 

x  =  maximum  alongwind  separation  distance  between  the  source 
'  A  and  any  point  on  the  line-of-sight 

tg  =  total  source  emission  time  (burn  time  for  quasi-continuous 
‘  smoke  sources)  in  minutes 

=  see  Equation  (2-50) 

u  =  see  Equation  (2-57) 

The  maximum  time  interval  from  the  smallest  value  of  t  greater  than 
zero,  calculated  for  all  sources,  to  the  largest  value  of  tp  for  all 
sources,  is  subdivided  into  the  number  of  times  the  user  specifies  that 
probability-of-detection  calculations  are  to  be  made.  The  probability 
of  detection  is  then  calculated  for  this  number  of  times  during  the  time 
interval  required  for  the  obscurant  from  all  sources  to  pass  through  the 
line-of-sight. 
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The  obscurant  for  any  specific  source  is  not  considered  in  the  calcula¬ 
tion  for  a  specific  time  t^  within  the  maximum  time  interval  if  the  inter¬ 
val  between  and  t  for  that  munition  does  not  encompass  t^.  It  is 
also  not  considered  if  there  is  no  point  on  the  line-of-sight  between 
the  target  and  observer  which  is  within  i5  0^  (see  Equation  (2-62)  for  oz) 
of  the  effective  height  of  the  Cloud  H{ t ^ }  for  that  source.  It  should  be 
noted  that  all  obscurant  sources  are  considered  in  the  calculation  of 
integrated  concentrations  between  the  iCl*  light  source  and  jC*1  cloud 
element  for  all  values  of  t^  within  the  maximum  time  interval.  If  no 
obscurant  from  any  source  passes  through  the  line-of-sight  between  the 
observer  and  target,  no  probability-of-detection  calculation  is  made  and 
a  message  stating  this  condition  is  printed. 

Addition  of  a  Routine  JJsed_t_o  Div ide 
the  Llne-of-Slght  into 
Cloud  Elements 

The  obscuration  code  in  AMSORB  contains  an  automated  procedure 
for  dividing  the  line-of-sight  between  the  target  and  observer  into 
elements  for  Che  cloud  brightness  calculations.  The  code  uses  portions 
of  the  transport  and  dispersion  models  to  define  the  lateral  and  vertical 
boundaries  of  the  obscurant  cloud  at  any  distance  from  the  source  as  5 
standard  deviations  from  the  cloud  centerline.  The  effective  line  source 
within  these  cloud  boundaries  is  defined  and  divided  Into  20  equally- 
spaced  cloud  elements  for  each  calculation  time.  If  the  line  length  of 
each  element  d-i  is  less  than  i  meter  at  this  point,  the  number  of  elements 
is  deci eased  until  the  length  di  is  1  meter  and  this  number  of  elements 
is  used  in  the  calculations  of  cloud  brightness.  If,  on  the  other  hand, 
dt  is  greater  than  l  meter,  the  trapezoidal  rule  of  integration  is  used 
in  combination  with  the  dispersion  model  to  calculate  the  integrated 
concentration  along  the  effective  cloud  length.  The  number  of  cloud 
elements  Is  then  doubled  to  40  and  a  similar  procedure  followed.  If  d£  is 
less  chan  1  meter,  the  number  of  elements  is  reduced  until  d£  equals  1 
meter  and  this  number  of  elements  is  used  in  the  cloud  brightness  calcu¬ 
lations.  If  df.  exceeds  1  meter  for  40  elements,  the  integrated  concentration 


is  calculated  and  compared  with  the  integrated  concentration  calculated 
using  20  elements.  If  the  change  in  integrated  concentration  is  less  than 
1.5x10  ^  percent,  the  number  of  elements  is  decreased  until  the  change  is 
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1.5x10  percent.  The  ei.tire  procedure  of  doubling  the  elements,  checking 

the  length  of  di.,  and  comparing  the  percentage  change  with  the  criteria  of 
-3 

1.5x10  percent  is  continued  until  either  a  1-meter  spacing  or  a  spacing 

-3 

yielding  less  than  a  1.5x10  percent  change  in  the  integrated  concentration 
is  achieved.  Tor  the  trial  calculations  we  have  performed,  the  procedure 
has  led  to  changes  in  the  cloud  brightness  of  less  than  0.1  percent  and 
even  smaller  changes  in  the  calculated  value  of  contrast  C  used  in  the  nro- 
bability-of-detection  calculation. 

Add i_t  ion_  of  Coding  Permitting  Simultaneous 
Calculations  for  Different  Types 
of  Obscurant  Sources 

The  Obscuration  Model  code  in  AMSORB  has  the  capability  of 
performing  probabil ity-of-detection  calculations  when  sources  containing 
different  types  of  obscurants  are  released  at  the  same  or  different  times. 

Lu  the  program,  the  line-of-sight  integrated  concentrations  are  calculated 
separately  for  each  source,  which  allows  the  source  characteristics  and 
the  time  the  source  is  released  to  be  different  for  eacli  type  of  source. 

The  line-of-sight  concentrations  and  elemental  concent . ations  along  the 
line-of-sight  are  summed  tor  each  source  of  the  same  type  and  the  total 
attenuation  calculated  from  the  expression 


t 


(2-153) 


where  there  are  t*  different  types  of  obscurants  and  the  attenuation  coef¬ 
ficient.'.  u  are  known  for  each  type  of  obscurant.  The  cloud  brightness 
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of  the  plume  containing  these  multiple  obscurants  is  calculated  according 
to  Equation  (2-146)  and  takes  into  account  the  difference  in  scattering 
and  attenuation  properties  of  each  type  of  obscurant. 

Modification  of  the  Sky- 
Partitioning  Routine 

The  user  must  input  the  number  and  location  of  light  sources 
to  the  ASLSOM  and  SOMI  visibility  codes.  The  number  of  light  sources 
and  their  location  are  obtained  by  using  a  separate  program.  For  use 
in  the  automated  AMSORB  program,  we  developed  data  tables  in  the  mass 
storage  file  containing  the  location  of  sky  segments  based  on  their  posi¬ 
tion  with  respect  to  a  standard  reference  frame.  For  each  data  case, 
this  standard  reference  frame  is  translated  to  a  new  reference  frame 
that  is  oriented  with  respect  to  the  line-of -sight  between  the  target 
and  observer.  These  data  tables  are  based  on  the  results  of  extensive 
target  brightness  calculations  made  with  the  SOMI  routine  using  hundreds 
of  sky  partitions.  The  results  obtained  using  the  simplified  procedure 
outlined  below  are  generally  more  accurate  than  those  obtained  using  the 
SOMI  routine. 

The  procedure  used  to  define  the  light  source  locations  described 
in  Section  4.4  is  based  on  the  division  of  the  sky  into  either  6  partitions 
of  equal  solid  angle  when  the  sun  is  shining  on  the  face  of  the  target  or 
16  partitions  when  the  sun  is  behind  the  target.  In  both  cases,  another 
partition  of  equal  solid  angle  is  formed  whose  centroid  has  Che  actual 
coordinates  of  the  sun's  position  and  has  the  sun's  brightness.  The  sun's 
partition  may  overlap  some  of  the  ocher  sky  partitions,  but  this  causes 
no  problem  in  the  calculations  since  each  partition  is  handled  independently. 
The  sun  is  thus  placed  in  its  own  partition,  resulting  in  either  a  total 
of  7  or  l?  partitions  being  used  in  the  calculations.  The  sun's  position 
in  the  sky  is  automatically  obtained  from  a  routine  contained  in  the 
mesoacale  model  routine,  which  calculates  the  poo it  ion  as  a  function  of 
date,  time  of  day,  and  the  longitude  and  latitude.  The  eloud  cover  is 
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also  passed  to  the  obscuration  model  from  the  mesoscale  model  routine. 

If  the  sky  is  overcast,  the  sun  is  not  considered  in  the  calculation  and 
16  partitions  or  segments  are  used  in  the  calculations.  The  program  will 
not  currently  handle  nighttime  scenarios  and  some  changes  will  be  required 
so  the  user  can  input  the  brightness  and  position  of  the  moon  or  other 
light  sources  such  as  flares. 

Modification  of  the  Test  for 
Large  Obscurant  Concentrations 

The  ASLSOM  smoke  visibility  model  contains  a  check  to  determine 
if  the  existing  concentration  of  obscurants  is  sufficient  to  obscure  the 
target  without  consideration  of  the  cloud  brightness,  thus  eliminating  the 
necessity  of  performing  the  time-consuming  cloud  brightness  calculations. 

We  have  changed  the  location  of  this  test  in  AMSORB  so  that  the  line- 
integrated  concentration  between  the  observer  and  target  is  checked  prior 
to  any  target  and  background  brightness  calculations.  Computer  runs  using 
data  from  smoke  trials  also  indicated  that  the  critical  level  of  concentra¬ 
tion  required  to  obscure  the  target  without  considering  cloud  brightness 
needed  to  be  increased.  Cioud  brightness  calculations  are  not  considered 
necessary  in  AMSORB  when  the  product  of  the  line-integrated  concentration 
and  the  coefficient  of  attenuation  is  equal  to  or  greater  than  11,  rather 
than  the  value  of  7  used  in  ASLSOM. 

Modification  of  the  Test  for 
Small  Obscurant  Concentrations 


The  ASLSOM  model  also  contains  a  check  to  determine  if  the 
amount  of  obscurant  between  the  target  and  observer  is  so  small  that  the 
probability  of  detection  is  at  least  95  percent.  The  ASLSOM  program 
first  calculates  the  target  brightness,  reduced  by  the  staoke  cioud,  and 
then  calculates  the  attenuation  oi  the  target  brightness  dee  to  the  smoke 
between  the  target  and  observer  without  making  the  cloud  brightness 
calculation.  If  the  target  brightness  is  reduced  less  than  5  percent 
by  the  smoke  betvceu  the  target  ar.d  observer,  the  ASLSOM  code  concludes 
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that  the  probability  of  detection  is  greater  than  95  percent  and  does 
not  perform  the  cloud  brightness  calculation.  In  cases  where  the  re¬ 
flectivities  of  the  target  and  background  are  similar,  even  a  small  amount 
of  obscurant  can  cause  a  reduction  in  the  probability  of  detection.  For 
this  reason,  the  Obscuration  Model  Routine  in  AMSORB,  in  addition  to  the  above 
tests  on  transmission  and  brightness,  computes  the  probability  of  detection 
given  by  Equation  (2-137)  with  the  cloud  brightness  set  to  zero  before 
the  conclusion  is  reached  that  the  probability  of  detection  is  greater  than 
95  percent  and  cloud  brightness  need  not  be  calculated. 


Replacement  of  the  Output  Routine 


The  output  routine  of  ASLSOM  was  changed  to  reflect  the  types 
of  calculations  made  in  AMSORB,  and  the  MS3  Routine  includes  options  for 
selecting  graphical  output  showing  time  profiles  of  line-of-sight  integrated 
concentration,  probability  of  detection  and  light  transmittance. 


2.4.3  Limitations  of  the  MS3  Routine 


The  AMSORB  routines  have  certain  limitations.  As  mentioned 
above,  the  current  version  of  the  program  does  not  contain  sufficient  data 
in  the  mass  storage  data  files  to  perform  completely  automated  calculations 
for  nighttime  scenarios,  and  for  making  probability-of-detection  calcula¬ 
tions  for  lasers  and  infrared  scenarios.  Thus,  when  such  calculations 
are  desired,  the  AMSORB  program  requires  the  user  to  specify  the  same  in¬ 
puts  as  those  required  by  the  ASLSOM  code.  Aiso,  because  of  the  lack  of 
data,  we  have  been  unable  to  check  completely  the  operation  of  these  por¬ 
tions  of  the  code.  The  adverse  weather  module  of  ASLSOM  has  been  included 
in  AMSORB,  but  also  has  not  been  checked. 


The  mass  storage  data  files  contain  source  data  tor  various 
smoke  munitions,  munition  produced  dust,  dust  produced  by  moving  vehicles, 
muzzle  blast  smoke  and  dust,  burning  vegetation,  burning  buildings,  and 
burning  vehicles.  Ue  have  Included  data  for  only  those  smoke  munitions 
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tested  recently  at  Dugway  Proving  Ground.  The  mass  storage  files  contain 
Mie  scattering  fractions  for  only  WP  smoke  in  the  visible  spectrum  and 
HC  smoke  at  the  10.6-micrometer  wavelength  in  the  infrared.  Because  of 
this  limitation,  Mie  scattering  fractions  for  WP  smoke  are  currently  used 
for  all  obscurants  in  the  visible  spectrum  unless  other  values  are  input 
by  the  user.  Attenuation  coefficients  for  all  obscurants  at  all  wavelengths 
of  interest  are  also  not  available,  although  values  exist  for  visible  light 
for  all  of  the  obscurants  mentioned  above. 

Most,  if  not  all,  of  these  limitations  are  caused  by  the  lack  of 
suitable  data  for  use  in  the  program  and  will  be  removed  when  such  data 
become  available  and  are  placed  in  the  mass  storage  data  file.  Details  of 
the  data  obscuration  calculations  now  included  in  the  mass  storage  data  file 
are  given  in  Section  4.4. 
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SECTION  3 

BATTLEFIELD  ENVIRONMENT  SOURCE  CHARACTERISTICS 


The  source  characteristics  of  obscurant  sources  are  required 
for  use  in  the  MS3  and  Battlefield  Environment  (BF.C)  Routines  of  AMSORB. 
This  section  contains  a  description  of  the  source  characteristics  and 
models  used  in  developing  source  characteristics  (SORDAT  routine)  for 
smoke  munition  sources,  dust  clouds  produced  from  munition  bursts  and 
from  moving  vehicles,  smoke  and  dust  clouds  from  muzzle  blasts,  and 
smoke  from  burning  vegetation,  buildings  and  vehicles.  In  some  cases, 
source  characteristics  models  developed  for  the  E-0  SAEL  program  library 
(Duncan,  e_t  al_,  1979)  are  included  in  the  AMSORB  program  and  are  used  to 
develop  source  inputs  in  the  SORDAT  routine  of  AMSORB. 


3.1  SMOKE  MUNITION  SOURCE  CHARACTERISTICS 

Table  3-1  gives  the  source  inputs  required  by  the  dispersion 
models  to  treat  smoke  munition  sources.  The  user  need  only  supply  the 
source  coordinates,  the  detonation  time,  the  direction  of  travel  with 
respect  to  grid  north  of  the  projectile  prior  to  detonation,  and  the 
type  of  smoke  muni r ion  being  used  when  the  munition  is  selected  from  one 
of  the  fifteen  tnuntcions  or  subsuuni c ions  for  which  data  are  available  in 
the  AMSORB  mass  storage  data  file.  Source  input  data  for  these  smoke 
sources  are  shown  in  Table  3-2.  The  parameters  tor  all  sources,  except 
the  bulk  white  phosphorus  (WE)  sources  were  experimentally  determined 
during  wind-tunnel  trials  at  Dugway  Proving  Ground  (Peterson,  1978). 

The  parameters  for  the  bulk  UP  munitions  (last  four  munitions  in  Table 
3-2)  approximate  values  obtained  by  using  CL  profiles  measured  at  short 
distances  from  static-firings  of  the  munitions  in  field  trials  (UPC, 
1978b)  to  develop  the  emission  characteristics.  The  orientation  of  the 
submunition  patterns  of  the  155  mm  hexac loroethanc  (HC)  smoke  projectile, 
and  the  !!>!>  an  XM82‘»  WP  smoke  projectile  are  dependent  on  the  flight 
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TABLE  3-1 

SMOKE  SOURCE  MODEL  INPUTS 


Parameter 

Symbol 

Parameter  Description 

Program  Default 
Value 

x,  y,  z 

Source  coordinates  (m) 

M 

o 

Total  mass  before  ignition  (g) 

MYF 

Munition  yield  fraction 

YE 

Yield  factor  accounting  for  the  effects 
of  moisture  in  the  air 

t 

o 

Time  of  detonation  (s) 

CB 

Total  burn  time  (min) 

ICAT 

Type  of  smoke 

1 

A,  B,  C,  D 

Coefficients  of  the  source  -elease  rate 
expression 

0 

xo 

Standard  deviation  of  the  alongwind  con¬ 
centration  distribution  at  the  source  (m) 

I 

i 

Q 

c 

Standard  deviation  of  the  crosswind  con¬ 
centration  distribution  at  the  source  (m) 

i 

i 

0 

zo 

Standard  deviation  of  the  vertical  con¬ 
centration  distribution  at  the  source  (ra) 

i 

uf 

c 

Effective  rate  of  heat  release  from  quasi- 
coat  inuous  sources  (cal  s~  ) 

! 

0 

UIR 

Direction  ot  the  projectile  flight  with 
respect  to  grid  north  (degrees) 

;  o 

TYPSKK 

Code  specifying  the  type  of  smoke 

i  °  i 

o. — - — — - _ — - ..  ... — 
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CQ 

C6 

o 

cn 

5 


lu 

O 

H 

ui 

o 

gj 

o 

u 
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< 

a 


UJ 
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H 
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W 
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C  •  cvoK  Oj0  for  these  munitions 

Tt;<  s  munition  is  composed  of  92  submuni tions  which  have  been  modeled  by  17  submunitions  as  shown  in  Table  3- 3(b). 
The  3y0  chosen  for  each  submunit ion  is  presented  in  Table  3-3(’o).  All  submunitions  have  Oxoc4.33  and  cj^-0.12 
A  best  estimate  vrss  chosen  for  the  yield  fraction  since  the  calculated  data  fit  for  this  munition  has  not 
been  comp  l e t  ed . 


direction  of  the  projectile  before  impact.  The  character istic  submuni¬ 
tion  impact  pattern  for  the  155  mm  HC  projectile,  which  contains  3-MI 
and  1-M2  canisters,  is  shown  in  Table  3—  3(a)  and  the  characteristic 
pattern  and  other  source  characteristics  for  the  155  tnra  XM825  WP  pro¬ 
jectile  are  given  in  Table  3—  3(b)  (see  Carter,  Dumbauld  and  Rafferty, 

1079  for  the  basis  of  these  patterns).  In  the  AMSORB  program,  the  user 
specifics  the  ground  impact  centroid  location  of  the  submunition  pattern 
and  the  direction  of  the  projectile  flight  (measured  clockwise  from 
north).  The  program  automatically  calculates  the  coordinates  of  each 
submuni  cion  based  on  the  angle  <J>  in  Table  3-3  that  a  line  joining  the 
submunition  position  and  the  centroid  of  the  pattern  makes  with  the 
direction  of  the  projectile  flight,  using  the  distance  between  the 
impact  centroid  location  and  the  submunition  also  given  in  Table  3-3. 

The  HC  munition  produces  elemental  2inc  which,  in  contact  with 
moisture  In  the  air,  produces  the  obscuring  agent,  hydrated  zinc  chloride. 
Similarly,  the  phosphorus  emitted  by  WP  and  UP  munitions  produces  phos¬ 
phoric  acid.  In  the  program,  the  yield  factor  YF  accounting  for  the 
production  of  these  materials  is  linearly  interpolated  from  the  tabular 
values  of  YE  as  a  function  of  relative  humidify  given  in  Table  3-4. 

The  empirical  values  of  the  yield  factors  for  wr  and  RP  were  obtained 
from  a  Dugway  Proving  Ground  Report  (DPC,  1978a)  and  those  for  HC  muni¬ 
tions  from  the  report  by  Petersen  (19/8).  The  relative  humidity  mea¬ 
sured  at  the  surface  weather  station  nearest  the  source  location  is  used 
in  the  calculation  of  V F . 

3.2  EXPLOSIVE  MUNITION  SOURCE  CHARACTERISTICS 

The  explosion  of  an  artillery  projectile  or  detonation  of  a 
mine  produces  dust  clouds  from  the  dust  and  debris  lofted  iuro  the 
atmosphere  by  the  force  Of  the  explosion.  The  amount  of  dust  produced 
depends  ott  the  size  aud  type  of  munition,  the  height  of  the  detonation, 
the  thickness  of  vegetative  cover  and  the  type  uf  soil.  A  completely 
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TABLE  3-4 


YIELD  FACTORS  FOR  WP,  HP  AND  HC  SMOKE  SOURCES 
USED  IN  THE  PROGRAM 


Relative 

WP  and 

HC 

Sources 

Humid i ty 
(percent) 

RP 

Sources 

0 

3.42 

2.09 

5 

3.42 

2.09 

10 

3.57 

2.50 

20 

3.87 

2.98 

30 

4.17 

3.34 

40 

4.46 

3.57 

50 

4.76 

3.90 

60 

5.06 

4.29 

70 

5.36 

4.86 

80 

5.66 

5.63 

100 

5.66 

L 

5.63 

satisfactory  theoretical  description  of  all  aspects  of  an  exploding 
munition  has  not  been  developed.  However,  empirical  relationships 
describing  the  initial  dust  cloud  produced  in  an  explosion  have  been 
developed  principally  by  Kbersole,  £t.  al.  (1979)  and  Zirkind  (1979)  from 
data  measured  during  field  trials  at  Uugway,  Utah  (t)PG,  1978c),  Fort 
Sill,  Oklahoma  (UPC,  1978d) ,  White  Sands,  New  Mexico  (l.indberg,  1979) 
and  Crafenwohr,  Germany  (l.undien,  1979).  The  inputs  required  to  develop 
the  source  parameters  for  the  dispersion  model  are  described  in  Table  3-5. 
The  only  quantities  the  user  needs  to  input,  however,  are  the  source 
coordinate,  the  time  of  detonation,  the  munition  type,  the  soil  type  and 
the  sod  depth. 

The  amount  of  dust  available  for  formation  of  a  dust  cloud  is 
produced  from  the  dust  and  debris  propelled  into  the  air  by  the  e. ‘.plo¬ 
sion.  The  mass  of  dust  lofted  into  the  atmosphere  and  remaining  air¬ 
borne  is  given  by 


M,, 


HP 


s  c 


0-1) 


whcic  o  l-1*  the  soil  density  and  V  is  the  crater  volume.  The  term  n  in 
Y  s  c 

Kquat ton  (3-1)  is  ati  efficiency  factor  for  soil  suspension  in  the  atmo¬ 
sphere  aud  is  considered  by  Ebersole,  td.  ai_.  (1979)  to  be  about  0.10  for 
common  soils.  The  majority  of  the  experimental  evidence  indicates  the 
dust  produced  in  the  detonation  is  partitioned  into  a  buoyant  Just  cloud 
and  a  nonbuoyant  dust  skirt.  The  amount  of  mass  M. ,  found  in  the 
buoyant  cioud  and  the  amount  of  mass  X,  in  the  nonbuoyant  dust  skirt  are 
given  by 


M 


(3-2) 


where  the  values  of  were  estimated  by  Zirkind  (1979)  to  be 


U.90  i  *  l;  buoyant  cloud 


0.10 


l 


2;  nonbuoyant  dust  skirt 


(3-3) 


TABLE  3-5 

SOURCE  INPUTS  FOR  THE  EXPLOSIVE  MUNITION 
SOURCE  MODEL 


Parameter 

Symbol 

Parameter  Definition 

Program  Default 
Value 

->  y 

Scarce  coordinates  (m) 

t 

0 

Time  of  detonation  (s) 

ITYPE 

Type  of  explosive  munition 

5 

W' 

Equivalent  charge  weight  of  TNT  (kg) 

6.98 

f 

Fraction  of  energy  available  for  crater 

formation 

0.6 

dD 

Detonation  depth  above  or  below  the 

surface  (m) 

0 

V  V 

Coefficients  for  determining  the  average 

0.503,  -0.954 

V  V 

scale  1  crater  radius 

0.450,  1.19 

V  V 

Coefficients  for  determining  the  average 

0.251,  -1.17, 

b2’  b3’  \ 

scaled  crater  depth 

-3 

U.494,  4.72,  3.34 

Ps' 

Soil  density  (g  cm  ) 

1.5 

a 

s 

Sod  or  vegetation  cover  depth  vm) 

0 

ISOLTP 

Soil  type  describing  the  soil  texture  and 
consis  tency 

4 

Studies  by  Ebersole,  et_  n_l.  (1979)  showed  that  vegetative  cover 
or  sod  depth  dg  reduces  the  crater  volume.  They  found  that  the  crater 
volume  for  an  explosive  munition  detonation  could  be  expressed  as 


2irr  ^d 
c  c 


1  -  -  — 
3  d 


where  r  is  the  crater  radius  and  d  is  the  crater  depth.  The  crater 
c  c  r 

radius  and  crater  depth  are  given  in  terms  of  a  scaled  crater  radius 
and  a  scaled  crater  depth  d^  according  to 


d  «  d  W 
c  c 


The  effective  energy  W  available  for  crater  formation  is  found  by 

t 

reducing  the  available  energy  W  from  the  explosive  charge  in  equiva¬ 
lent  kilograms  of  TNT  by  some  fraction  f  according  to 

I 

W  »  fW 


A  table  of  f  values  determined  by  Ebersole,  et  al .  are  given  in  Table 

3-6  as  a  function  of  type  of  charge  and  detonation  depth  ithove  or  below 

the  surface  of  the  soil.  A  negative  depth  indicates  the  detonation 

» 

occurred  beneath  the  surface  of  the  soil.  The  available  energy  W  in 
equivalent  kilograms  of  TNT  of  the  exploding  munition  required  in 
Equation  (3-7)  is  given  in  Table  3-7  for  various  munitions.  The  density 
of  the  soil  is  used  in  Equation  (31)  to  determine  the  mass  of  airborne 
dust.  The  soil  density  depends  on  soil  type  and  a  table  of  various  soil 
types  and  densities  is  given  in  Table  3-8.  The  scaled  crater  dimensions 
in  Equations  (3-5)  and  (3-6)  also  depend  on  soil  type.  The  empirical 
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FRACTION  OF  ENERGY  AVAILABLE  FOR  CRATER  FORMATION 


Detonation 

type 


Av til  able  for 
Crater  Formation,  f 


Detonation 
Depth  d.. 

On) 

1 

Description 

0 

Live  fire  tilted  at  30  degrees 

0 

Horizontal  live  fire 

-0.3 

Live  fire  tilted  at  30  degrees 

-0.6 

Live  fire  tilted  at  30  degrees 

0 

Static  bare  charge 

-0.3 

Static  bare  charge 

0.3 

Live  fire 

- 

User  entered  bare  charge 

user  entered  artillery 
projectile 

TABLE  3-8 
SOIL  DENSITIES 


[ - 

1 

Soil  Category 

Soil  Density 
(g  era  ) 

Description 

1 

2.4 

Rock 

2 

2.0 

Clay  or  dry  cohesive  soils 

3 

0.7 

Dry  sandy  soils 

4 

1.5 

Loam,  sandy  soil  and  frozen  ground 

5 

2.0 

Soft  rock 

6 

1.8 

Wet  sand,  wet  loam  or  ice 

7 

1.8 

Wet  cohesive  soils  (not  saturated) 
and  snow 

10 

- 

User  entered  soil  type 

relationship  developed  by  Ebersole,  £t  al_.  (1979)  from  the  field  data 
for  the  scaled  crater  radius  is 


r 

c 


-a  + 
o 


2 

A  a„ A  +  3_  A 
c  2  c  3  c 


and  for  the  scaled  crater  depth  is 


d 

c 


-  b  +  b.  A  + 

o  1  c 


2  3  4 

b_  A  +  b,A  +  b. A 

2  c  3  c  4  c 


where 


A  - 
c 


(3-8) 


(3-9) 


(3-10) 


The  term  d^  in  Equation  (3-10)  is  the  detonation  depth  in  meters  above 
or  below  the  surface  and  can  be  obtained  from  Table  3-6  based  on  the 
detonation  type.  The  polynomial  coefficients  used  in  Equations  (3-8) 
and  (3-9)  are  given  in  Table  3-9  for  the  soil  categories  in  Table  3-8. 

The  dispersion  models  require  the  initial  dimensions  of  the 
buoyant  and  nonbuoyant  dust  clouds  produced  by  the  explosion  and  also 
require  the  energy  available  for  buoyant  plume  rise.  The  initial  radius 
of  the  buoyant  dust  cloud  at  aerodynamic  equilibrium  (approximately  20- 
30  milliseconds  after  detonation)  is  used  to  determine  initial  source 
dimensions.  This  buoyant  cloud  radius  in  meters  is  given  by  (after 
Ebersole,  ec  al . ,  1979) 


R 

o 


(3-11) 


where  is  the  air  density  in  grows  per  cubic  meter  and  pQ  is 

1225  g  ra-~*.  The  height  of  the  centroid  of  the  buoyant  cloud  depends 
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TABLE  3-9 


COEFFICIENTS  USED  IN  DETERMINING  AVERAGE  SCALED  CRATER  DIMENSIONS 

(after  Ebersole,  £t  al,  1979) 

(a)  Coefficients  for  Average  Scaled  Crater  Radius 


'  •  Coefficients 

Soil  - 

Category 

a 

0 

ai 

a2 

a3 

I,  2 

0.271 

-0.684 

0.390 

0.886 

3 

0.386 

-0.849 

0.367 

0.993 

4 

0.303 

-0.954 

0.450 

1.19 

5.  6 

0.629 

-1 .08 

0.264 

1.12 

7 

0.806 

-1.28 

-0.178 

0.852 

(b)  Coefficients  for  Average  Scaled  Crater  Depth 


Coefficients 

Soil 

Category 

bo 

bl 

b2 

b3 

b4 

1 

o.n: 

-0.477 

0.270 

1.84 

1.05 

2 

0.134 

0.571 

0.343 

2.24 

1.31 

3.  5 

0. 189 

-0.840 

0.447 

3.30 

2.10 

4 

0.251 

-1.17 

0.494 

4.72 

3.34 

6 

0.331 

-1.49 

0.579 

4.92 

3.13 

7 

0.449 

-1.82 

0.322 

4.11 

2.02 
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on  the  type  of  charge  used  and  is  higher  for  bare  charges  than  for 
artillery  projectiles.  Experimental  observations  by  Ebersole,  et  al^ 
(1979)  of  dust  clouds  produced  by  explosive  munitions  indicated  that  the 
initial  height  H  of  the  cloud  centroid  is 


I5Rq;  bare  charges  detonated  at  the  surface 
Rq;  artillery  projectiles. 


(3-12) 


where  t  e^ual  to  zero  is  the  detonation  time.  For  bare  charges  detonated 
beneath  the  surface,  the  height  is  given  by  Rq.  The  height  of  the  cen¬ 
troid  of  the  nonbuoyant  dust  cloud  is  assumed  to  be  zero.  The  initial 
standard  deviations  of  the  two  dust  clouds  are  given  in  terms  of  the 
initial  buoyant  cloud  radius  as 


o  =o  -a  =  R  / 2 . 1 5 ;  buoyant  cloud 
xo  yo  zo  o 


a 

xo 


zo 


-a  ■  1.155R 

yo  o 

=  .91  R 

o 


> 


nonbuoyant  dust  skirt  ^ 


(3-13) 


since  the  dust  skirt  is  assumed  to  be  a  cylinder  of  height  and  radius 
equal  to  2Ro-  Finally,  Che  energy  in  calories  available  for  pluae  rise 
is  given  as 


Qx  •  kU  (3-14) 


where  k  « 


l.lxlO6 


calories  per  kilogram  of  TNT. 
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3.3 


MOVING  VEHICLE  SOURCE  CHARACTERISTICS 


A  moving  vehicle  can  generate  dust  through  the  interaction 
of  the  wheels  or  track  with  the  ground.  The  dust  produced  is  generally 
small  particles  less  than  100  vim  in  diameter.  The  inputs  required  to 
develop  the  source  parameters  are  listed  in  Table  3-10.  The  only  quan¬ 
tities  the  user  needs  to  input,  however,  are  the  vehicle  location  prior 
to  moving  (for  example,  the  grid  boundary  for  a  vehicle  entering  the 
grid),  the  initial  start  time  of  vdil  le  movement,  the  total  time  of 
travel,  a  code  specifying  the  vehicle  type,  a  code  specifying  the  soil 
type,  the  direction  of  travel  of  the  vehicle  with  respect  to  grid  north 
and  die  speed  of  travel. 


Some  work  has  been  done  by  various  agencies  to  determine  the 
amount  of  dust  produced  by  a  moving  vehicle.  However,  Just  generation 
depends  on  a  wide  variety  of  factors  which  have  not  been  accurately 
modeled.  The  source  model  in  AMSORb  is  based  on  the  urrently  approved 
EPA  model  developed  by  Midwest  Research  Institute,  which  gives  the 
yearly  average  emission  rate  E  in  grams  per  vehicle  meter  as  (Mann  and 
Cowherd,  1977) 


E 


0.2283 


I 


(3-15) 


where  3^  is  the  silt  content  of  the  surface  soil  in  t  rcent  (percentage 

of  the  weight  in  particles  whose  diameters  are  less  clw  75  pin) ,  V 

_  [  c 
is  tiic  vehicle  spec!  (ms  ),  N'  is  the  number  of  wheels  on  the 

vehicle  and  is  the  number  of  days  during  the  year  when  precipita¬ 

tion  exceeds  0.025  cm  (0.01  in).  Since  the  user  of  AMSORB  is  not  inter¬ 
ested  in  average  emission  rates,  but  requires  current  emission  rates,  we 
suggest  Up  be  set  equal  to  0  for  dry  soil  and  to  365  for  wet  soil. 

The  silt  contents  as  a  function  of  soil  type  are  listed  in 
Table  3-11  and  were  assigned  based  on  work  reported  by  Cowherd  and 
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TABLE  3-10 

SOURCE  INPUTS  FOR  THE  MOVING  VEHICLE  SOURCE  MODEL 


Parameter 
Symbo 1 

Parameter  Description 

Program  Default 
Value 

x,  y 

Starting  coordinates  of  the  moving  source  (m) 

'o 

Start  time  of  vehicle  (s) 

t' 

Vehicle  travel  time  (s) 

0 

I  TYPE 

Code  specifying  the  type  of  vehicle 

1 

N' 

The  number  of  wheels  on  the  vehicle  or  a  best 

estimate  for  track  vehicles 

4 

ISOLTP 

Code  specifying  soil  type 

4 

S 

Silt  content  (percent  of  the  weight  of  part- 

C 

icles  in  the  surface  soil  which  have  diameters 

less  Chau  75  pm) 

25.7 

0 

Direction  ot  travel  with  respect  to  grid  north 

V 

(degrees ) 

0 

V 

Vehicle  velocity  (m  s  S 

13.41 

i 

-1 

V 

s 

Settling  velocity  of  the  soil  particles  (m  s  ) 

0 

1  w 

Number  of  days  during  the  year  with  .025  cm 

P 

(.01  in)  or  more  precipitation 

0 

CJ 

Standard  deviation  of  the  spherical  distribu- 

i 

t ion  ot  the  cloud  at  the  source  (m) 

0.93 

TABLE  3-11 


SILT  CONTENT  OF  SURFACE  SOILS 


Soil 

Category 

s 

c 

Silt  Concent 
(?) 

Descript  ion 

1 

0 

Uock 

2 

95 

Clay  or  dry  cohesive  soils 

3 

18 

Dry  sandy  soils 

4 

25.7 

Loam,  sandy  soil  and  frozen  ground 

5 

12.0 

Soft  rock 

6 

o 

Wet  sand,  wet  loam  or  ice 

7 

0 

Wet  cohesive  soils  (not  saturated) 
and  snow 

10 

L 

User  entered  soil  type 
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Hendriks  (1977)  and  Miller  (1979).  The  emission  iate  depends  on  the 
number  of  wheels  per  vehicle.  A  table  of  the  vehicle  types  and  their 
number  of  wheels  is  given  in  Table  3-12.  The  number  of  wheels  repre¬ 
senting  tracked  vehicles  were  arbitrarily  assigned.  The  dispersion 
model  for  moving  vehicles  is  described  in  Section  2.1.  S 


3.4  MUZZLE  BLAST  SOURCE  CHARACTERISTICS 


The  muzzle  blast  of  a  cannon  or  mortar  can  generate  both  dust 
and  combustion  product  particles  which  act  as  obseurw.ts  on  a  battle¬ 
field.  Particulates  are  generated  from  the  combustion  products  of  the 
propellant  ir.  t  he  gun  after  the  gases  have  emerged  from  the  gun  tube. 

This  propellant  cloud  is  largely  water  vapor  formed  by  the  oxidation  of 
hydrogen  in  the  cloud,  but  also  contains  smaller  amounts  of  carbon 
particles  as  well  as  a  few  other  gases.  A  dust  cloud  is  also  generated 
by  the  interaction  of  the  muzzle  blast  shock  wave  with  the  ground.  The 
ground  ia  compressed  and,  as  the  rarefaction  portion  of  the  shock  wave 
passes,  the  surface  layer  of  dust  is  propelled  into  the  air.  The 
propellant  cloud  is  formed  almost  instantly  whereas  the  dust  cloud 
requires  about  6  seconds  to  form  and  reach  the  same  height  as  the  muzzle. 
The  propellant  and  dust  clouds  combine  very  quick1  y  and  therefore,  we 
treat  them  as  one  cloud.  Only  the  particulate  yield  of  carbon  from  the 
propellant  combustion  will  be  considered  in  this  model.  For  the  visible 
wavelengths  this  a  fairly  good  approximation,  but  there  is  mounting 
evidence  (Stuebing,  c£  al . )  chat  transmission  in  the  infrared  at  10.6  ;ia 
is  affected  to  some  extent  by  the  present  in  the  cloud.  This  section 
describes  the  source  model  used  to  develop  parameters  describing  the 
dust  and  smoke  cloud  formed  by  the  muzzle  blast  of  a  cannon.  The  inputs 
needed  for  the  source  model  are  descrioed  in  Table  3-13.  However, 
only  the  location  of  the  gun  and  height  above  ground  of  the  muzzle,  the 
time  of  firing,  a  weapon  code  specifying  the  type  of  gun  being  fired  and 
a  code  specifying  the  soil  type  around  the  gun  are  required  unless  the 


TABLE  3-12 

VEHICLE  TYPES  MODELED  IN  THE  VEHICLE  MOVEMENT  PRODUCED 
DUST  SOURCE  MODEL 


Vehicle  Description 


Estimated  values  tor  track  vehicles. 


Number  of  Wheels 
N' 


1/4  ton  Jeep 

5/4  ton  M8B0 

2  1/2 

ton  truck 

Dump 

truck 

5  ton 

tractor 

5  ton 

tractor  with  trailer 

APC 

Tank 

User 

entered  vehicle  type 

TABLE  3-13 


SOURCE  INPUTS  FOR  THE  MUZZLE  BI.AST 
SOURCE  MODEL 


Parameter 
Symbo  L 

Parameter  Definition 

Program  Default 
Value 

x,  y 

Source  coordinates  (m) 

h 

Height  of  muzzle  above  ground  (m) 

2 

t 

0 

Time  of  firing  (s) 

I  TYPE 

Code  specifying  the  weapon  type 

5 

I) 

t 

Muzzle  blast  flush  diameter  (m) 

a 

xo 

Standard  deviation  of  the  aior.gwind  con- 

cent  rat  ion  distribution  at  the  source  (m) 

8.07 

0 

Standard  deviation  of  the  crosswind  con- 

yo 

cent  rat  ion  distribution  at  Lite  source  (m) 

2.03 

0 

Standard  deviation  of  the  vertical  con- 

20 

ceotration  distribution  at  the  source  (m) 

2.03 

M 

Mass  of  combustion  products  formed  in  the 

combustion  oi  the  propellant  charge  in  the 

guti  (g) 

M 

si 

Mass  of  soil  propelled  into  the  air  (g) 

tSOLTP 

Code  specifying  the  type  of  soil  at  the 

weapon  location 

4 

►>u 

Soil  density  (g  cm  *) 

1.5 

SS 


us:-  wishes  to  use  a  type  of  weapon  not  considered  in  the  input  selec¬ 
tion  list. 


The  amount  of  smoke  particulates  and  dust  in  the  cloud  can  be 
determined  from  the  expression 


M  »  M  +  M  (3-17) 

o  sc 


where  M  is  the  mass  of  soil  propelled  into  the  air  and  is  the  mass 
of  carbon  particulates  formed  from  the  combustion  products  of  the  pro¬ 
pellant  in  the  gun.  The  mass  of  dust  propelled  into  the  air  was  calcu¬ 
lated  by  using  the  assumption  made  by  Gibson  (1978)  that  the  area  of 
soil  affected  by  the  muzzle  blast  is  determined  by  the  area  of  ground 
intercepted  by  a  spherical  shock,  wave  emanating  from  the  muzzle  and 
striking  the  ground  at  an  incidence  angle  of  60  degrees.  Gibson  also 
assumes  that  only  10T  ot  the  first  millimeter  of  soli  in  this  area 
becomes  airborne.  Thus,  the  mass  of  soil  injected  into  air  by  the  shock 
wave  Is  given  by 


M 

s 


G. 0001 


(  h 

it an (60°) 


■> 


(3-18) 


where  ^  is  the  soil  density  (g  m  ')  i :.u  h  ; the  height  of  the 
muzzle  bore  hole  abo<'_  the  ground  (e>.  rhe  s  -  i.  1  Irn-sity  is  given 
in  Tabic  J-5  as  a  i'unetioft  of  -  LI  type. 

The  particulate  mass  M  in  Equation  (3-1?)  is  determined  bv 
assuming  that  the  mass  or  carbon  in  the  propellant  is  converted  to 
particles  and  the  rest  of  she  conut  it  uents  are  emitted  as  vapoft  that  dc 
not  coalesce  on  any  particulates  after  emission,  the  grams  of  ear boa 
formed  in  combustion  of  TNT  (according  to  Thompson,  1979)  is  228  grass 
of  carbon  per  kilogram  of  TNT .  The  particulate  mass  for  the  different 
types  of  weapons  was  obtained  by  determining  the  amount  of  propellant 
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used  as  an  average  load  in  the  gun  from  the  U.  S.  Army  Field  Mannual 
FM9-13.  The  mass  of  propellants  containing  nitrocellulose  or  nitro¬ 
glycerin  was  converted  to  equivalent  mass  of  TNT  using  a  scaling  factor 
determined  by  comparing  the  energy  available  from  the  different  types  of 
propellant  powder  to  that  of  TNT.  A  table  of  mass  of  carbon  particu¬ 
lates  by  weapon  type  is  given  in  Table  3-14.  The  size  of  the  initial 
cloud  can  be  determined  from  the  dimension  of  the  flash  diameter  D^, 
according  to  Gibson  (1978),  by  the  following  relationships: 


zo 

u 


xo 


yo 

3a. 


yo 


(3-19) 


If  the  dimensions  of  the  flash  in  length  and  width  ere  known,  then 
the  standard  deviation  of  the  concentration  at  the  source  can  be  deter¬ 
mined  from 


l 


a 


zo 


a 


yo 


o 

xo 


1 


and  these  values  can  be  input  to  the  program  instead  of 


(3~2C) 


3.5  BURNING  BRUSH/ VEGETATION  SOURCE  CHARACTERISTICS 

Brush  or  vegetation  can  be  ignited  by  explosions  or  incendiary 
munitions  and  the  subsequent  smoke  produced  in  the  combustion  can  sub- 
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TABLE  3-14 

ESTIMATED  PARTICULATE  MASS  PRODUCED  FROM  COMBUSTION 
OF  PROPELLANT  IN  A  WEAPON 

Description 


81  mm  mortar  14.8 
4.2  in  mortar  47 .9 
120  mm  tankgun  1391 
10 f*  mm  tank  HEP  298.7 

1468 
3080 


Weapon  Type 


Mass  of  Carbon 
Particulate,  Mc  (g) 


155  mm  artillery 
8  in  artillery 
User  entered  type 


stantially  degrade  visibility.  The  particulates  emitted  from  burning 
vegetation  depend  at  a  minimum  on  the  amount  and  kind  of  vegetation 
present,  the  moisture  content,  local  meteorological  conditions  and  the 
availability  of  oxygen.  Table  3-15  gives  the  inputs  used  in  the  AMSORB 
burning  vegetation  source  model.  However,  only  the  location  of  the 
centroid  of  the  burning  brush,  the  time  combustion  starts,  a  code 
specifying  the  type  of  regional  vegetation  involved,  and  the  area 
involved  in  the  burning  need  to  be  specified  for  the  model  since  the 
program  will  automatically  assign  values  to  ether  required  parameters  from 
data  tables  contained  in  the  program. 

The  amount  of  particulate  matter  emitted  by  the  fire  over  its 
entire  burn  time  is  given  approximately  by  the  following  relationship 
used  by  the  Environmental  Protection  Agency  (Vatavuk  and  Yamate,  1975): 

M  =  P  •  %c  •  A  (3-22) 

of 

where  A  is  the  area  of  the  fire,  is  the  fuel  loading  of  vegetation 
and  P  is  the  particulate  yield.  The  particulate  yield  for  most  fires 
varies  between  it)  and  90  g/kg  of  fuel  consumed  according  to  Gomez, 

££  al.  (1979)  and  Eccleston,  ec_  al_ .  (1974).  Tims,  an  average  value  for 
brush  and  gra  >  fires  of 


P  «*  18  g/kg  of  fuel 


(3-23) 


is  used  in  the  computer  program.  The  average  fuel  loadings  by  regions 

of  the  country  are  presented  in  Table  3-16  and  are  the  same  as  those 

»  » 

used  by  Vatavuk  and  Yamate.  The  effective  heat  release  rate  (cal  s  ) 
emitted  by  the  fire  is  given  as 

%  “  \  •  lf/60  (3-24) 
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TABLE  3-15 


SOURCE  INPUTS  FOR  THE  BURNING  BRUSH/VEGETATION 
SOURCE  MODEL 


Parameter 

Symbol 

Parameter  Description 

Program 
Default  Value 

x,  y 

Source  coordinates  of  the  centroid  of  the 
burning  vegetation  (ra) 

t 

0 

Time  combustion  starts  (s) 

WYF 

Yield  fraction  indicating  the  amount  of 
combustion  mass  which  remains  airborne  as 
obscuring  particles 

1.0 

CB 

Total  burn  time  (min) 

A,B,C,D 

Coefficients  of  the  source  release  rate 
expression 

1.0, 0,0,0 

0 

xo 

Standard  deviation  of  the  alongwind  concen¬ 
tration  distribution  at  the  source  (m) 

2.89 

0 

Standard  deviation  of  the  crosswind  concen- 

2.89 

yo 

tration  distribution  at  the  source  (m) 

0 

Standard  deviation  of  the  vertical  concen- 

2.89 

tration  distribution  at  the  source  (m) 

1TYPRC 

Code  specifying  ;he  type  of  region  and 
Indicating  the  type  of  vegetation  present 

3 

Estimated  fuel  loading  (kg  m  ^ ) 

6.7 

rf 

Average  fuel  loading  corresponding  to  the 
burn  rate  R  (kg  m"“) 
a2  -1 

Burn  rate  (a  min  ) 

0.8 

140 

P 

Pollutant  yield  from  the  fire/grams  of 
particulate  emitted  per  kilogram  of  fuel 
consumed 

18 

hk 

Heat  released  by  the  ec^busion  of  the 
particular  fuel  (cal  kg”) 

3.9xl06 

Area  of  the  burning  vegetation  for  total 
burn  (a^) 

100 

93 


f 


TABLE  3-16 


ESTIMATED  AVERAGE  FUEL  LOADINGS 
(after  Vatavuk  and  Yaraate,  19/3) 


Area  or  Region 

Estimated  Average 

Fuel  Loading,  % 

(kg  n  >  f 

Rocky  Mountain  Croup 

8.3 

Region  1 :  Northern 

13.5 

Region  2:  Rocky  Mountain 

6.7 

Region  3:  Southwestern 

2.2 

Region  4:  Intermouatain 

4.0 

Pacific  Group 

4.3 

Region  5:  California 

4.0 

Region  6:  Pacific  Northwest 

13.5 

Region  10:  Alaska 

3.6 

Coast  a 

13.5 

Interior 

2.5 

Southern  Group 

2.0 

Region  8:  Southern 

2.0 

Eastern  Croup 

2.5 

North  Central  Group 

2.5 

Region  9:  Conifers 

2.2 

Hardwoods 

2.7 

User  entered  type 


-1  2  -1 

where  H_  is  the  heat  emitted  (cal  kg  ),  R  is  the  burn  rate  (m  min  ), 

h  d 

is  the  average  fuel  loading  of  the  burning  vegetation  and  the  factor 
of  60  i9  included  to  convert  to  seconds.  If  the  heat  released  is  not 
input  by  the  program  user,  the  program  uses  a  standard  value  for  pine 
wood  and  dry  grass  given  by  Oomcz,  et_  a_l.  (1979)  ns 

Hj,  =  3.9  x  10^  cal/kg  of  fuel  consumed  (3-25) 


The  burn  rate  R_  in  Equation  (3-24)  is  determined  from 

D 


R 


B 


(3-26) 


where  ie  is  an  average  fuel  loading  for  the  vegetation  and  R_  is  the 

l  D 

measured  average  burn  rate  for  an  average  fuel  loading  1  .  The  ABSORB 

2-1  _2  * 

program  uses  values  of  140m  min  and  0.8  kg  m  respectively  for 
Ry  and  based  on  measurements  reported  by  Daniels,  ca  ,  (1975). 

The  burn  time  in  minutes  can  be  determined  from  the  expression 

* 

tB  •  \/Rb  (3-27) 

where  A.^.  is  the  total  area  of  vegetation  involved  in  the  combustion. 

When  the  standard  deviations  of  the  alonguind,  erosswind  and  vertical 
concentration  distributions  at  the  source  are  not  input  by  the  user,  the 
program  determines  those  values  according  to  the  following  relationship 

“  v  "  '  ,  *  /TK?W  0-28) 

assuming  that  a  uniform  distribution  of  vegetation  exists  over  the  entire 
burn  area. 
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3.6 


BURNING  BUILDINC  SOURCE  CHARACTERISTICS 


The  AMSORB  program  uses  an  approach  to  define  source  character¬ 
istics  for  burning  buildings  similar  to  the  approach  used  for  burning 
vegetation  described  in  Section  3.3  above.  Table  3-17  describes  the 
inputs  used  in  the  burning  building  source  model.  However,  only  the 
location  of  the  centroid  of  the  burning  building,  the  time  combustion 
starts,  the  building  height,  the  area  of  the  building  and  the  mass  of 
combustibles  involved  need  be  input  since  the  program  will  automatically 
assign  values  to  the  other  parameters  based  on  the  data  tables  present 
in  the  program. 

The  amount  of  particulate  matter  emitted  by  the  fire  over  its 
entire  burn  time  is  calculated  using  Equation  (3-22).  In  this  application 
for  burning  buildings,  A^  is  the  area  of  the  building,  Is  the  fuel 

loading  in  the  building  and  l>  is  the  particulate  yield.  The  effective 
heat  release  rate  in  calories  per  second  emitted  by  the  fire  is 

found  using  Equation  (3-24).  If  the  heat  released  is  not  input,  the  pro¬ 
gram  uses  the  value  tor  pine  wood  and  dry  grass  (see  Equation  (3-25)  ). 
The  burn  rate  used  to  determine  the  heat  release  rate  as  well  as  the 
burn  time  is  given  by  Equation  (3-26).  The  burn  time  for  the  building 
is  given  by  Equation  (3-27).  The  standard  deviations  of  the  alongwind, 
crosswind  and  vertical  concentration  distributions  at  the  source  can  be 
input.  Ii  they  are  not,  the  program  defines  the  alongwind  and  crosswind 
standard  deviation  from  the  expression 


where  is  the  building  area.  The  vertical  standard  deviation  at  the 

source  for  a  building  of  height  is  given  by 


o 


so 


% 


-29) 


(3-30) 


TABLE  3-17 


SOURCE  INPUTS  FOR  THE  BURNING  BUILDING  SOURCE  MODEL 


Parameter 

Symbol 


Parameter  Description 


Program  Default 
Va  Uie 


x.  y 


Source  coordinates  of  the  centroid  of 
the  building  (ra) 


HYP 


Yield  fraction  indicating  the 

amount  of  combustion  mass  which  remains 

airborne  as  obscuring  particles 


1.0 


h 


B 


Building  height  (m) 


2.0 


M 

o 


Mass  of  combustible  material  in  the 
building  (kg) 


A,  B,  C,  D 


Coefficients  of  the  source  release  rate 
expression 


1.0,  0,  0,  0 


o 

xo 


Standard  deviation  of  the  alongwind  con¬ 
centration  distribution  at  the  source  (m) 


0.631 


o 

yo 


Standard  deviation  of  the  crosswind  con¬ 
centration  distribution  •  the  source  (m) 


0.631 


o 
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Standard  deviation  of  the  vertical  con¬ 
centration  distr ibuation  at  the  source  (m) 


1.15 


A-. 


Estimated  fuel  loading  in  the  building 
(kg  m"') 

Average  fuel  loading  for  the  burn  race 
Rg(kg  m-2) 

Burn  rate  for  the  average  fuel  loading 
specified  by  I£.  (a*  min“i) 

Pollutant  yield  from  the  fire  (grams  of 
particulate  emitted  per  kilogram  of 
fuel  consumed 

Heat  released  by  the  combustion  of  the 
particular  fuel  (cal  kg"*) 

•> 

Area  of  the  burning  building  (m*) 


20.0 


o.a 


140 

18.0 


3.9  x  106 


3.75 


3.7 


SOURCE  CHARACTERISTICS  FOR  BURNING  VEHICLE 


Table  3-18  lists  the  general  source  inputs  required  for  a  burn¬ 
ing  vehicle  as  well  as  Che  default  parameter  values  entered  in  the  con- 
puter  program.  The  default  values  are  based  on  data  for  a  burning 
automobile  obtained  by  Gomes,  et  al_,  (1979).  The  inputs  in  Table  3-18 
are  those  required  by  the  quas i-continuous  source  model  because  it  is 
assumed  that  the  total  vehicle  burn  time  is  on  the  order  o;  10  minutes 
or  longer.  if  the  program  user  lias  source  input  parameter  values  for  a 
burning  vehicle  that  are  more  appropriate  than  the  default  values,  they 
can  ho  entered  in  Che  program  and  used  in  the  diffusion  model  calculations. 

The  default  values  listed  in  Table  3-18  reflect  the  details  of 
the  data  obtained  by  Gomez,  e_t  al.  (1979)  for  a  burning  automobile.  For 
example,  the  fraction  MYF  has  been  set  equal  to  1.0  because  only  the 
total  airborne  mass  of  particulates  was  measured.  Also,  the  default 
values  assigned  to  the  coefficients  of  the  expressions  for  the  source 
release  rate  (A- 1.0;  B,  C,  D=0)  reflect  the  assumption  of  a  constant 
release  rate  for  the  total  burn  time  of  4b  minutes.  The  default  value 
assigned  to  the  standard  deviations  of  the  concentration  distributions 
of  the  source  were  all  set  equal  to  0.866  meters  under  the  assumption 
that  the  horizontal  and  vertical  cloud  dimensions  were  3  meters. 


J.  8  USER- SPECIF  LED  SOURCE  MODELS  AND  INPUT'S 

Proyision  has  been  made  in  the  SORDAT  routine  of  AMSORB  for 
user-specified  source  models  which  can  be  applied  to  both  instantaneous 
and  quasi-eoncinuous  sources.  Table  3-19  lists  the  requisite  source 
parameters  which  must  be  input  by  the  user. 
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TABLE  1-19 

INPUTS  FOR  USER  SPECIFIED  SOURCE  MODELS 


Parameter 

Symbol 


Parameter  Description 


x,  y,  z 


Source  coordinates  (m) 


Total  mass  before  ignition  (g) 


MYF 


t 

o 


t 


B 


A,  a,  C,  D 


0 

xo 


0 

yo 


c 


TYPSMK 


Source  yield  fraction 
Initial  source  release  time  (s) 

Total  burn  time  (min) 

Coefficients  of  cbe  source  release  race 
express  ion 

Standard  deviation  of  the  alongwind  concentra¬ 
tion  distribution  at  Che  source  (in) 

Standard  deviation  of  the  crosswiud  concentra¬ 
tion  distribution  at  the  source  (in) 

Standard  deviation  of  the  vertical  concentra¬ 
tion  distribution  at  the  source  (m) 

Total  heal  released  from  instantaneous  sources 
(cal ) 

Effective  rate  of  heat  release  trom  quasi- 
coni  iuuous  sources  (cai  s-^) 

Code  identifying  type  of  obscurant  (1,  3) 


100 


Obscuration  calculations  require  attenuation  coefficients  and 
Mie  scattering  coefficients  which  depend  on  the  type  of  obscurant.  Th„ 
user  must  therefore  enter  a  parameter  value  for  the  type  of  obscurant, 
denoted  by  the  code  identification  TYPSMK,  in  the  program.  A  value  of 
1  for  TYPSMK  refers  to  HC  smoke  and  a  value  of  2  refers  to  WP  smoke. 

The  program  data  base  contains  values  of  attenuation  coefficients,  Mie 
scattering  fractions  and  the  yield  fractions  YF  for  each  of  these  ob¬ 
scurants.  A  value  of  3  entered  for  TYPSMK  specifies  a  user-spe<* <  f { 
obscurant  and  the  user  must  input  appropriate  attenuation  coefficient- 
and  Mie  scattering  fractions.  Because  no  relative  humidity  yield  tactions 
can  be  entered  in  the  program  for  user-specified  obscurants,  thr  program 
assumes  that  the  yield  factor  YF,  accounting  for  the  effects  of  moisture, 
is  unity. 
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SECTION  4 

AMSORB  MODEL  ROUTINES 


The  computer  program  AMSORB  consists  of  a  number  of  separate 
routines,  each  of  which  contributes  directly  or  indirectly  to  a  battle¬ 
field  environment  solution.  This  section  describes  these  routines  and 
defines  the  basic  model  parameters  contained  in  the  routines.  Section 
4.1  describes  the  Mixing-Layer  Analysis  Routine  (MIXLYR)  that  utilizes 
rawinsonde,  USAF-CWC  wind  and  temperature  data,  and  surface  meteorological 
data  contained  in  mass  storage  data  base  files  to  calculate  initializa¬ 
tion  parameters  required  by  the  Mesoscale  Wind-Field  Model.  The  Battle¬ 
field  Source  Characteristics  Routine  (SORDAT) ,  described  in  Section  4.2, 
processes  user-provided  battlefield  source  input  data  to  generate  a  data 
base  file  containing  source  locations  and  source  characteristics.  The 
input  parameters  required  by  the  various  dispersion  models  are  discussed 
in  Section  4.3.  The  dispersion  routines  use  surface  meteorological  data 
and  mixing-layer  parameters  provided  by  the  Mesoscale  Wind-Field  Model 
as  direct  inputs,  and  source  data  from  user-provided  inputs.  Section 
4.4  discusses  the  Obscuration  Model  and  required  model  input  parameters. 
The  Battlefield  Environment  Routine  (BEC)  described  in  Secton  4.5  uses 
the  battlefield  source  location  and  characteristics  data  base  file  gener¬ 
ated  by  SORDAT  and  meteorological  parameters  generated  by  the  Mesoscale 
Wind-Field  Model. 


4.1  MIXING-LAYER  ANALYSIS  ROUTINE 


The  Mixing- Layer  Analysis  Routine  employs  a  relatively  simple 
approach  in  utilizing  the  various  types  of  data  available  from  the  AMSORB 
mass  storage  data  base  to  calculate  initialization  parameters  for  the 
Mesoscale  Wind-Field  Model.  The  decision  to  use  a  simple  approach  was 
made  so  that  the  effects  of  logic  decisions  on  the  derived  parameters 


could  be  easily  traced  and  analyzed.  The  Mixing-Layer  Analysis  Routine 
is  thus  designed  to  perform  a  modified  "single-station  analysis"  rather 
than  a  more  complex  analysis  of  the  pressure,  wind  and  temperature 
fields  over  the  entire  region  of  interest.  It  should  also  be  noted  that 
some  of  the  logic  decisions  made  in  accepting,  rejecting  and  analyzing 
the  data  are  somewhat  arbitrary.  Where  possible,  we  have  attempted  to 
study  the  effects  of  the  various  decision  processes.  However,  until  the 
routine  is  used  to  predict  dosage  and  concentration  fields  under  a  wide 
variety  of  meteorological  conditions,  the  full  implications  of  many  of 
the  decision  processes  cannot  be  determined. 

Key  features  of  the  Mixing-Layer  Analysis  Routine  (MIXLYR)  are 
shown  in  the  schematic  logic  diagram  of  Figure  4-1.  The  key  features 
are: 


•  Rawinsonde  data  from  the  closest  station  with  a  radius  of 
less  than  100  kilometers  from  the  release  point  and  less 
than  three-hours  old  at  the  time  of  request  are  used 
without  modification  to  obtain  estimates  of  the  input 
parameters  to  the  Mesoscale  Wind-Field  Model 

•  Rawinsonde  data  measured  at  the  closest  station  less  tnan 
100  kilometers  from  the  release  point,  greater  than 
three-hours  old  and  less  than  twelve-hours  old  are  sub¬ 
ject  to  modification  using  OWC-predieted  and  surface 
winds  and  temperatures 

•  Climatological  estimates  of  mixing- layer  depths  and 
surface  wind  speeds  are  used  as  inputs  if  all  rawinsonde 
data  are  more  than  twelve-  hours  old  or  if  there  are  no 
stations  within  a  100  kilometer  radius  of  the  release 
point 
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MODEL 


•  Climatological  estimates  are  also  used  if  the  rawinsonde 
data  are  more  than  six-hours  old  and  the  rawinsonde 
surface  temperature  and  850-millibar  temperatures  have 
not  been  modified  because  GWC  and  surface  data  are  either 
unavailable  or  inadequate 

•  GWC  data  are  used  to  modify  the  rawinsonde  temperatures 
and  winds  at  the  700-  and  850-millibar  levels;  rawinsonde 
temperatures  and  winds  at  the  700-millibar  level  are 
replaced  by  the  GWC-predicted  temperature  for  that  level 
if  the  GWC-predicted  temperature  for  the  hour  of  interest 
is  within  5  degrees  Celsius;  rawinsonde  temperatures  and 
winds  at  850  millibars  are  replaced  if  the  temperature 
from  the  reporting  surface  station  nearest  the  time  and 
point  of  release  is  within  5  degrees  Celsius  of  tne  CWC- 
predicted  temperature  at  the  850  millibar  level  for  the 
hours  of  interest 

•  Surface  temperature  and  wind  data  from  the  closest 
station  less  than  50  kilometers  from  the  release  point, 
less  Chan  three-hours  old,  and  with  a  surface  pressure 
from  the  rawinsonde  observation,  are  used  to  modify  the 
surface  rawinsonde  observation  when  it  is  more  than 
three-hours  old 

•  Surface  wind  speeds  from  the  closest  station  to  the  point 
of  release  are  used  as  default  values  for  the  mean-layer 
wind  speed  if  the  rawinsonde  data  are  not  analysed 

Climatological  estimates  of  the  mixing- layer  depth  shown  in  Table  4-1 
are  used  as  default  values  where  the  rawinsonde  date  are  not  available 
or  are  not  analysed.  These  estimates  of  the  median  mixing  depths  as  a 
function  of  month  and  liae-of-day  are  based  on  the  work  of  Norton  and 


TABLE  4- 


Hoidale  (1975),  who  analyzed  8,236  rawinsonde  observations  made  at  WSKR 
between  1961  and  1972  using  the  technique  described  by  Holzworth  (1967), 

As  shown  In  Figure  4-1,  the  rawinsonde  data  are  analyzed  to 
determine  the  mixing-layer  depth,  mean-layer  wind  direction  and  speed  and 
reduced  gravity  factor.  In  developing  an  automated  procedure  for  deter¬ 
mining  the  mixing-layer  depth,  we  evaluated  four  methods  for  establishing 
that  a  surface-based  or  elevated  inversion  was  indicated  by  the  rawinsonde 
observations.  The  simple  criteria  associated  with  the  four  methods  for 
determining  the  base  of  an  inversion  were: 

Method 


1 


2 


3 


4 


Criteria 


to 

to 


>  .01 


f  >  0 

Az  - 


AO 

w 

-y-  >  0 

Az  - 


AT  . 

v  -4 

>  -5x10 

A  2  -- 


where 


M 

Az 


lapse  rate  of  potential  temperature 


AT 

Az 


lapse  rate  ot  temperature 
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lapse  rate  of  wet-bulb  potential  temperature 

lapse  rate  of  virtual  temperature 

The  MIXLYR  routine,  which  incorporates  these  criteria,  was  used  to 
analyze  twice-daily  rawinsonde  observations  from  Salt  Lake  City,  Utah, 
for  the  first  five  days  of  each  month  during  the  period  June  1972  through 
May  1973  and  for  all  the  month  of  January  1973.  The  heights  of  the 
bases  of  stable  layers  obtained  from  the  routine  were  then  qualitatively 
compared  with  those  selected  by  meteorologists  from  the  plots  of  the 
data.  All  the  criteria  produced  bases  of  stable  layers  in  agreement 
with  those  selected  by  the  meteorologists  for  cases  in  which  clearly- 
defined  stable  layers  were  present.  In  our  view,  however,  the  criterion 
based  on  the  lapse  rate  of  virtual  temperature  (Method  4)  performed  more 
consistently  than  the  other  criteria  when  stable  layers  were  not  as 
clearly  defined  and  for  all  seasons  of  the  year.  For  this  reason,  the 
procedure  finally  used  in  the  routine  is  based  on  this  criterion. 

The  routine  accepts  the  raw  or  modified  rawinsonde  data  and 
begins  the  selection  process  by  checking 


A0 


AT 


AT 


v,  n 


AT  . .  -  T 
v,n+l  v,n 


AZ 


Az 


n+1 


(4-1) 


in  the  first  height  interval  above  the  surface.  If  the  criterion  for 

-4 

the  base  of  a  stable  layer  is  met  (e.g.,  if  AT^  ^  /  Az^  :>  -5x10  ')  and  the 
height  interval  Az  is  greater  than  or  equal  to  100  meters,  the  routine 
defines  a  surface-based  inversion  and  sets,  by  default,  the  surface  mixing 
deptn  to  30  meters.  If  Az^  is  less  than  100  meters  but  ATy  j/Az^-Sxl 0'4 , 
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the  routine  then  checks  AT  „/Az„.  If  AT  »/Az0  >  -5x10 

v,2  2  v,2  2  — 

(Az^  +  kz^)  2.  1^0  meters,  the  routine  checks  the  quantity 


AT' 

Az 


T  0  -  T 
v,n+2  v,n 

z  -  z 
n+2  n 


v>  3  v,l 
2  3  '  Z1 


(4-2) 


If  (At/Az)  _>  -5x10  the  routine  again  defines  a  surface-based  inversion. 
If  no  surface  based  inversion  is  found,  the  routine  increments  n  and 
proceeds  through  the  same  operations  until  the  base  of  an  elevated 


inversion  is  found  or  until  z 


exceeds  3000  meters  AGL.  If  z 


u  ,  i  vawvvuo  o  (wu  •  c,  ,  , 

n+1  n+1 

exceeds  3000  meters  AGL,  the  surface  mixing  depth  is  defaulted  to  3000 

meters. 


After  the  depth  of  the  surface  mixing  layer  H  is  established, 

tn 

the  routine  analyzes  the  rawinsonde  data  between  the  surface  and  H  to 


m 


obtain  estimates  of  the  mean-layer  wind  components  and  the  reduced 
gravity  factor  using  the  expressions: 


“h 


m 


m 


u  (z  -  z  ) 
n  n+1  n 


n=l 


£  <zn+i  ‘  '•„> 


n-1 


(4-3) 


H 

Vn  ^Zn+1  zn^ 

n=l 

H 

m 

L  <,  -  «  ) 

n+ 1  n 

n*»  1 


(4-4) 
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i  - 


(4-5) 


Where  TI,max  is  Che  maximum  temperature  in  the  stable  layer  above  H  . 
In  the  routine,  g*  is  restricted  so  that  if  it  exceeds  0.3  it  is  set™ 
to  0.3  and  if  less  Limn  0.1  is  set  to  0.1. 


4.2  BATTLEFIELD  SOURCE  CHARACTERISTICS  ROUTINE 

The  Battlefield  Source  Characteristics  Routine  (SORDAT)  is 
designed  to  develop  the  source  characteristics  for  a  large  number  of 
sources  required  to  describe  a  battlefield  environment  problem.  From  a 
relatively  few  inputs  supplied  by  the  user,  SORDAT  generates  all  neces¬ 
sary  source  characteristics  required  for  dispersion  modeling  purposes. 
Additionally,  SORDAT  saves  all  source  data  in  the  battlefield  source 
location  and  characteristics  data  base  for  subsequent  use  by  the  Battle¬ 
field  Environment  Routine  (BEC) .  Routine  SORDAT  generates  source  char¬ 
acteristics  tor  any  of  the  various  types  of  sources  discussed  in  Section 
3.  As  shown  in  figure  4-2,  the  routine  reads  user-provided  data,  which 
includes  t.ne  source  location  on  the  battlefield,  the  time  of  day  the 

source  emission  began,  the  type  of  source  or  source  category,  and  other 
data  that  may  be  required  for  a  particular  source  category.  Depending 
on  the  source  category,  the  routine  then  computes  or  retrieves  the 
source  characteristics  required  for  dispersion  modeling  from  data 
tables.  These  source  characteristics  include  the  initial  cloud  dimen¬ 
sions,  the  total  source  emission  time,  the  mass  rate  of  emission,  and 
the  cloud  release  height  or  parameters  required  to  calculate  the  cloud 
height.  Some  source  categories  require  additional  characteristics;  for 
example,  the  vehicle  movement  source  category  requires  the  direction, 
speed  and  time  of  travel  of  the  vehicle  movement.  SORDAT  then  stores 
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FIGURE  4-2.  Battlefield  Source  Characteristics  Routine  Flow  Diagram. 
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FIGURE  4-2.  (Cotit  imiej) 


20 


the  source  data  into  the  data  base  and  loops  back  to  road  another 
source  data  card  set.  After  all  user-provided  source  data  have 
been  processed,  SORDAT  optionally  prints  tables  of  the  source  charac¬ 
teristics  data  based  on  source  category. 


4.3  DISPERSION  MODEL  ROUTINES  AND  METEOROLOGICAL  PARAMETERS  REQUIRED 

BY  THE  DISPERSION  MODELS 

The  meteorological  parameters  required  by  the  dispersion 
models  are  given  in  Table  4-2.  Of  the  meteorological  parameters  shown 
in  the  table,  the  user  need  only  supply  the  program  with  values  of  r 
and  A  if  the  dispersion  models  in  Section  2.3.2  art*  being  used  and  if 
these  parameter  values  respectively  differ  from  the  default  values  of  1 
and  0.  When  the  MSI  Routine  of  AMSORB  is  used,  the  subscript  n 
equals  unity  and  only  the  values  of  the  parameters  in  Table  4-2  at  the 
source  location  are  used  in  the  calculations. 

The  turbulence  parameters  u'  and  o'.  ,  their  power-law 

»  il  C.  *  ll 

coefficients  for  height  dependency  tn  and  q,  and  the  power-law  coef¬ 
ficient  p  are  not  routinely  measured  at  meteorological  stations. 

n 

A  means  of  assigning  appropriate  values  to  these  meteorological  param¬ 
eters  is  therefore  required,  Swanson  and  Cramer  (1963)  made  a  com 
prehens ivs  study  of  the  height  dependence  of  0^  as  a  function  of  wind 
speed  and  tirae-of-day  using  measurements  from  a  62-»cter  tower  at  WSMR. 

We  have  used  the  results  of  this  study  in  conjunction  with  a  stability 
elassif Icatioti  system,  similar  to  the  Pasquill  definition  of  stability 
categories  (Turner,  1964)  used  by  the  Environmental  Protection  Agency, 
to  select  and  assign  values  to  these  parameters.  Figure  4-3  is  a  sche¬ 
matic  diagram  of  the  logic  sequence  used  by  the  routine  iu  defining  the 
parameters  attd  the  use  of  the  parameters  in  the  dispersion  models. 

Key  features  of  the  procedure  shown  in  Figure  4-3  are: 


TABLE  4-2 


METKOROI.OCTCAL  INPUTS  FOR  Til" 
DISPERSION  MODELS 


Parameter 

Definition 

NR 

Net  radiation  index 

H 

m  ,n 

Depth  of  the  surface  mixing  layer  at  the  n1*1 

point  in  the  trajectory  (m) 

li 

Mean  layer  wind  speed  in  the  surface  mixing  layer 

n 

at  the  n1^1  point  on  the  trajectory  (m  s~l) 

°A  n 

Standard  deviation  of  the  wind  azimuth  angle  at 

A,n 

the  nt!'  point  on  the  trajectory  (radians) 

m 

I’ower-law  coefficient  used  to  describe  the  height 
dependency  of  o  in  the  surface  mixing  layer 

A « tl 

o' 

Standard  deviation  of  the  wind  elevation  angle  at 

F « n 

the  n1*1  point  on  the  trajectory  (radians) 

q 

power-law  coefficient  used  to  describe  the  height 

dependency  of  0  ,  in  the  surface  mixing  laver 

l. ,  n 

i> 

Power-law  coefficient  used  to  describe  the  height 

U 

dependency  of  wind  speed  in  the  surface  mixing 
laver 

rt 

Lateral  cloud  expansion  coefficient 

H 

Vertical  cloud  expansion  coefficient 

r 

Fraction  of  material  reflected  at  the  surface 

A 

Washout  coefficient  or  fraction  of  material 
removed  by  precipitation  scavenging  per  unit 
t  ime 

RH 

Relative  (tumidity  (percent)  used  to  compute  the 
yield  factor  (YP)  for  smoke  munitions 
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•  Solar  altitudes  for  the  release  point  and  time  of  request 
are  calculated  using  an  algorithm  suggested  by  Woolf 
(1968) 

•  Til*.'  routine  then  seeks  the  most  recent  cloud  and  ceiling 
data  from  the  reporting  station  nearest  to  the  release 
point 

•  If  there  are  no  stations  reporting  cloud  and  ceiling  data 
within  a  50-kilometer  radius  from  the  release  point  or 

the  data  are  more  tiian  three-hours  old,  the  routine  selects 
a  default  net  radiation  index  from  Table  4-3 

•  When  cloud  data  are  available,  the  routine  calculates  an 
initial  net  radiation  index  based  on  solar  altitude  (also 
see  Table  •<-:)  and  adjusts  the  index  based  on  the  time- 
ot-day,  tc-cal  cloud  cover,  and  ceiling  height 

•  The  net  radiation  index  and  the  wind  speeds  from  the 
Mesuseale  Wind-Field  Model  for  the  release  point  and 
points  along  the  cloud  trajectory  are  used  to  select 
turbulence  input  parameters  and  wind  power-law  coef¬ 
ficients  tor  the  dispersion  models 

«  Values  of  the  surface  mixing  depth  H  for  each  i»E  * 

f9  i  n 

point  along  the  cloud  trajectory  arc  obtained  from  the 
Mesoscale  Wind-Field  Model  solution  using  the  inter¬ 
polation  procedure  given  by  Equation  (2-6) 

Table  4-4(a)  shows  values  of  the  ten-minute  standard  devia¬ 
tions  of  the  wind  asimuth  angle  in  degrees  min)  and  associated 

A 

values  of  t lie  power-law  coefficient  a  for  describing  o  variation 

A 

with  -height  above  the  surface  based  on  the  studv  of  Swanson  and  Cramer. 
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TABU-  4-3 

NF.T  RADI  AT  t  DM  TJJMCFS 


Solar  Altitude 

Daytime 

In  Degrees 

Index 

V  >  GO 

4 

35  <  i]  %  GO 

3 

15  <  r)  <  35 

2 

0  <  TJ  15 

1 

Nighttime 

N/A 

Default 

Index 


TABLE  4-4 


Ti'RBLLTLCE  PARAMETERS  FOR 
DISPERSION  MODELS 


(a)  Ten-Minute  Standard  Deviations  of  the  Wind  Azimuth  Angle  in  Degrees 


(a.)  and  Corresponding  Values  of  the  Power-Law  Exponent  m 
A 


Wind  Speed , 
at  5  Meters) 
(m/sec) 


Net  Radiation  Index 


m 

°A 

m 

°A 

m 

°A 

m 

-.109 

26 

-.109 

21 

-.133 

15 

-.112 

-.  109 

22 

-.  128 

16 

-.123 

11 

-.  112 

122 

16 

-.  123 

13 

-.  103 

9 

-.120 

-.  105  | 

11 

-.117 

11 

-.117 

8 

-.148 

-.100 

9 

-.117 

8 

-.  110 

8 

-.148 

aA  m  aA  m 


189,8.  6 


(b)  Standard  Deviation  of  the  Wind  Elevation  Angle  In  Degrees  (a  )  and 
Corresponding  Values  of  the  Power -Law  Exponent  q 


Wind  Speed 
at  5  Meters 
(m/sec) 


Net  Radiation  Index 


u  <  1  8.7 


E 


10  8.  7 


0  7.0  .07 


1  <  u  <  3  8.7  .  10  7.4  .  13  id.  3  .07 


3  <  u  <  5  6.  3  .  15  5.3  .15  4.3  .  13  3.  0  0  3.  0  0  2.6  -.20  2.6 


<u<  7  4.7  .20  3.7  .20  3.  7  .15  2.  7  0  2.7  0  2.2  -.15  2.2  -.  15 


3.0  .25  3.0  .20  2.7  .20  2.7  0  2.7  0  2.2  -.  15  2.2 


— 

-2 

°E 

q 

2.9 

-.30 

2.9 

-.25 

2.6 

-.20 

2.2 

-.15 

2.2 

-.  15 

nr* 


Values  of  the  standard  deviation  of  the  wind  elevation  angle  in  degrees 
and  the  power-law  coefficient  q  describing  the  variation  of  with 

height  above  the  surface  are  shown  in  Figure  4-4 (b).  These  values  of 
or  were  assigned  under  the  assumption  that  (after  Cramer,  et  al . , 

1972)  turbulence  is  isotropic  for  short  averaging  times  and  that  the 
one-fifth  power  law  can  be  used  as  a  scale  factor  to  adjust  0^  for 
averaging  times  less  than  10  minutes.  Thus, 


°E  -  OA{T=i0  min}  (M> 


1/5 


(4-6) 


oa(t=10  min}  /2 . 99 


(4-7) 


The  values  of  the  power-law  coefficient  q  were  assigned  under  the  assump¬ 
tions  that: 

•  Under  stable  conditions  (NR  =  -1,  -2) 


°E,Z  -  °E,ZrUr 


"P. 


(4-8) 


Under  near-neutral  conditions  (NR  ®  0,1) 


°E,z  —  °E,z 


(4-9) 


«  Under  unstable  conditions  (NR  a  2,3,4) 


E,z  E,zr  VZj^ 


■»ZP  V«D  ' 


0.3  -p 


(4-10) 
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The  values  of  the  wind  power-law  coefficient  as  a  function 

of  net  radiation  index  and  mean  layer  wind  speed  used  in  the  routine  are 
shown  in  Table  4-5.  These  values  are  similar  to  those  suggested  for  use 
at  WSMR  by  Swanson  and  Hoidale  (1962)  and  for  Dugway  Proving  Ground  by 
Cramer,  et^  al .  (1972). 

The  plume-rise  models  described  in  Section  2.3.1  and  2.3.2 
above  require  estimates  of  the  ambient  air  temperature  and  density  and, 
under  stable  conditions,  the  lapse-rate  of  potential  temperature.  If 
surface  meteorological  data  are  available,  the  routine  automatically 
selects  the  ambient  air  temperatures  and  densities  according  to  the 
procedures  shown  in  Figure  4-3  and  the  lapse-rate  of  potential  temper¬ 
ature  is  calculated  from  the  rawinsonde  sounding  established  using  the 
procedure  described  in  Figure  4-1.  If  no  surface  data  are  available, 
the  default  values  of  surface  temperature  shown  in  Table  4-6  and  a 
default  air  density  of  1204.2  grams  per  cubic  meter  are  used.  The 
temperatures  in  the  table  are  based  on  the  monthly  normal  and  maximum 
temperatures  at  El  Paso  extracted  from  the  report  "Climatic  Atlas  of  the 
United  States"  published  by  the  U.  S.  Department  of  Commerce  in  June 
1968.  The  values  shown  for  a  net  radiation  index  of  -1  and  -2  corre¬ 
spond  to  the  normal  minimum  .  temperature  for  each  month.  Similarly,  the 
values  shown  for  net  radiation  indices  of  0  and  1  and  for  indices  of  2 
and  3  respectively,  correspond  to  the  normal  average  and  normal  maximum 
temperatures  for  each  month.  Finally,  the  surface  temperatures  for  a 
net  radiation  index  of  -4  were  obtained  by  taking  an  average  of  the 
recorded  maximum  temperature  and  the  normal  maximum  temperature  at  El 
Paso  for  each  month. 

If  no  rawinsonde  data  are  available,  the  routine  selects 
default  values  of  the  lapse  rate  of  potential  temperature  from  Table  4-7 
based  on  the  net  radiation  index  and  the  calculated  wind  speed  at  the 
effective  source  height.  A  default  value  of  1  is  always  used  for  the 
vertical  cloud  expansion  coefficient  3*  A  default  value  of  1  for  the 
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TABLE  4-5 


WIND  POWER-LAW  EXPONENT  p  FOR 
THE  DISPERSION  MODELS 


Mean  Layer 
Wind  Speed 
u  (m/sec) 

Net  Radiation  Index 

4 

3 

2 

1 

0 

-1 

-2 

u  <  1 

.2 

.20 

.20 

.20 

.20 

.3 

1  <  u  <  3 

.17 

.20 

.20 

.20 

.20 

.25 

3  <  u  <  5 

.15 

.15 

.  17 

.17 

.20 

.20 

.20 

5  <  u  <  7 

.  10 

.10 

.15 

.15 

.  15 

.15 

.15 

7  <  u 

.05 

.10 

.  10 

.10 

.  15 

.15 

.15 

Net  Radiation  Index 


Month 

4 

3,2 

1,0 

-1,-2 

Jan 

l~ 

292 

286 

279 

272 

Feb 

296 

290 

282 

275 

Mar 

300 

294 

286 

277 

Apr 

303 

299 

290 

282 

May 

308 

304 

295 

287 

Jun 

312 

308 

300 

292 

Jul 

312 

308 

301 

294 

Aug 

310 

307 

300 

293 

Sep 

308 

304 

297 

289 

Oct 

303 

299 

291 

283 

Nov 

297 

292 

284 

275 

Dec 

293 

287 

280 

272 

TABLE  4-7 

DEFAULT  VALUES  OF  THE  LAPSE  RATE  OF  POTENTIAL 
TEMPERATURE  $  (deg  K/m) 


Wind  Speed 
u{h} 

(m  sec"1) 

Net  Radiation  Index 

1 

0 

— 

-1 

-2 

0  -  1.4 

0.005 

0.015 

0.  030 

0.040 

1.5  -  2.9 

0.003 

0.010 

0.020 

0.030 

3.0  -  4.  9 


0.001 


0.005 


0.015 


0.020 


lateral  cloud  expansion  coefficient  a  is  used  for  instantaneous  sources 
and  a  value  of  0.9  is  used  for  quasi-continuous  and  continuous  sources. 

The  dispersion  models  described  in  Section  2.3.2  require  a 
yield  factor  (YF).  The  factor  YF  is  calculated  as  function  of  relative 
Inanity  (RH)  by  the  dispersion  models.  If  nv_  relative  huraidy  value  is 
available  from  the  Data  Base,  the  routine  selects  a  default  value  sho\m 
in  Table  4-8  based  on  the  time  of  day  and  season.  The  values  used  in 
Table  4-8  are  for  El  Paso  and  are  obtained  from  the  "Climatic  Atlas 
of  the  United  States"  published  by  the  U.  S.  Department  of  Commerce  in 
June  1968. 


4.4  OBSCURATION  MODEL  ROUTINE  AND  MODEL  PARAMETERS 

The  function  of  the  Obscuration  Model  Routine  is  to  calculate 
the  probability  of  detecting  a  target.  Figure  4-4  is  a  program  logic 
flow  diagram  of  the  Obscuration  Model.  Before  a  probability  of  detection 
value  can  be  computed,  several  sets  of  calculations  must  be  performed, 
as  shown  in  Figure  4-4.  First,  the  line- of-sight  (LOS)  integrated  con¬ 
centration  from  the  observer  to  the  target  is  calculated  for  all  contrib¬ 
uting  obscurant  sources.  The  Obscuration  Model  Routine  references  the 
Dispersion  Model  Routines  in  order  to  compute  integrated  concentration 
values.  Prior  to  performing  the  more  time  consuming  calculations  such 
as  cloud  brightness,  the  routine  checks  to  determine  if  sufficient  ob¬ 
scurants  are  present  to  reduce  the  transmittance  through  the  cloud  to  a 
level  at  which  targets  displaying  the  largest  possible  contrast  are 
always  obscured.  If  it  is  determined  that  this  condition  exists,  the 
routine  indicates  the  target  is  obscured  and  returns  control  to  the 
calling  program.  Otherwise,  the  routine  continues  processing  and  com¬ 
putes  the  brightness  of  the  target  as  viewed  from  the  observer.  The 
Dispersion  Model  Routines  are  again  referenced  for  computing  the  observer- 
to-background  LOS  integrated  concentration  and  the  brightness  of  the 
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Season  2200  -  0300 

Winter  36 

Spring  31 


Summer 


48 


TABLE  4-8 


UES  OF  RELATIVE  HUMIDITY  (PERCENT)  FOR 
MPUTING  THE  YIELD  FACTOR  (YF) 


K  i 


\ 


i 

'fc  i 


yy 


1 

i 


FIGURE  4-4. 


Obscuration  Model 


OBSCURANT 
DISPERSION 
MODEL  ROUTINE 


SET 

YES 

PROBABILITY  OF 

DETECTION  >  57. 

FLAG 

*« 


OBSCURANT 
DISPERSION 
MODEL  ROUTINE 


OBSCURANT 
DISPERSION 
MODEL  ROUTINE 


,ne  Fiou  Diagram. 
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(Continued) 


background  as  viewed  from  the  observer.  These  calculations  and  the 
calculations  of  the  target  brightness  are  then  used  to  determine  if  the 
LOS  integrated  concentration  is  sufficient  to  obscure  the  target  from 
the  observer.  If  not,  a  flag  is  set  indicating  that  the  observer  has  a 
95  percent  or  greater  probability  of  detecting  the  target  and  program 
control  returns  to  the  calling  program.  Should  an  obscuration  be  present, 
the  program  continues  processing  and  calculates  the  brightness  of  the 
obscurant  clouds  along  the  LOS,  and  finally  calculates  the  probability 
of  detection  as  shown  in  the  flow  diagram.  The  Obscuration  Model  Routine 
then  returns  to  the  calling  program. 

The  parameters  required  by  the  Obscuration  Model  described  in 
Section  2.4  are  given  in  Table  4-9.  Of  the  parameters  listed  in  the 
table,  onlv  the  coordinates  of  the  observer  and  target  need  be  input  by 
the  user  for  visible  scenarios  unless  the  user  has  specified  a  source  of 
obscurants  not  included  in  the  data  base.  In  addition,  a  code  for  the 
target  background  description  must  be  input.  As  noted  earlier,  the 
inputs  required  by  the  Obscuration  Model  are  the  same  as  those  required 
by  the  ASLSOM  code  for  infrared,  laser,  and  nighttime  scenarios,  or  when 
adverse  weather  is  to  be  considered  in  the  calculations.  Details  of  the 
inputs  required  by  cne  Obscuration  Model  Routine  for  calculations  under 
these  conditions  are  contained  in  Volume  II,  Section  A. 4  of  the  program 
users'  manual. 

Target  and  background  codes  or  target  and  background  reflec¬ 
tivities  for  the  wavelength  of  interest  as  well  as  the  target  size  must 
be  entered,  as  mentioned  above.  To  assist  the  program  user,  examples  of 
codes  currently  used  in  the  model  to  define  the  type  of  target  and  back¬ 
ground,  the  target  size,  and  the  target  and  background  reflectivities 
for  radiation  in  the  visible  spectrum  are  presented  in  Table  4-10. 

As  an  example  of  the  use  of  the  code  to  specify  a  scenario,  assume  the 
program  user  wishes  to  calculate  the  probability  of  detecting  a  vehicle 
with  ordinary  camouflage  paint  against  a  background  of  vegetation. 
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TABLE  4-9 


INPUTS  FOR  THE  OBSCURATION  MODEL 


Parameter 

Definition 

V  V  20 

Coordinates  of  the  observer  (tn) 

V  yt-  2t 

Coordinates  of  the  target  (m) 

d 

Average  diameter  of  the  target  (m) 

Rt 

Target  reflectivity  or  a  code  specifying 
target  type 

Background  reflectivity  or  a  code  specifying 
background  type 

B  , 
oi 

t  h 

Brightness  of  the  i  light  source  (candles 
ra*"  in  the  visible  region  of  the  spectrum 
and  W  m-2  steradian  “1  in  the  infrared) 

a  U) 

Attenuation  coefficient  for  wavelength 

X  (m2  rag*1) 

CB 

Threshold  contrast 

♦l 

£  h 

Scattering  angle  for  the  i  light 

source  (deg) 

F 

k2 

Mie  scattering  function  divided  by  the 
square  of  the  propagation  constant 

X 

Wavelength  of  light  or  a  code  specifying 

the  wavelength  or  region  of  interest  (urn) 
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TABLE  4-10 


OBSCURATION  ROUTINE  CODE  FOR  TARGET  AND  BACKGROUND 
REFLECTIVITIES  AND  TARGET  SIZES 


(a)  Target  Information 


Code 

Reflectivity 
in  the 

Visible  Region 

Average 

Target 

Diameter 

(m) 

Target  Type 

Reference 

for 

Reflectivity 

Values 

2 

0.08 

5 

Army  vehicle 
(ordinary  camouflage 
paint) 

Dolce  (1974) 

3 

0.30 

10 

Red  brick  building 

Handbook  of  Chemistry 

4  Physics  (1 9 5 7 , p2 74 5 ) 

4 

0.40 

10 

Wood  (pine)  buildi 

ng 

Handbook  of  Chemistry 
&  Physics  (1957  ,p2745) 

5 

0.55 

5 

Steel  target 

American  Institute  of 
Physics  Handbook  (1957) 

6 

0.69 

5 

Polished  aluminum 
target 

American  Institute  of 
Physics  Handbook  (1957) 

7 

0.72 

l 

White  painted 
sign 

Handbook  of  Chemistry 

4  Phys ics  (1957 ,p2745) 

8 

0.72 

10 

White  painted 
i  building 

Handbook  of  Chemistry 

4  Physics  (1957,p2745) 

(b) 

Jackground  Information 

Code 

Reflectivity 
in  the 

Visible  Region 

Background  Type 

Reference  for 
Reflectivity  Values 

2 

0.08 

Soil  (damp) 

Finklestein  (1964) 

1 

0.15 

Vegetation 

Dolce  (1974) 

4 

0.40 

Granite 

(estimated  value) 

5 

0.93 

Snot; 

American  Institute  of 
Physics  Handbook  (1957) 

6 

l  .0 

Sky 

Dolce  (19/4,  p56) 

Inspection  of  Table  4-10  shows  a  code  of  2  for  the  target  and  3  for  the 
background  provide  the  requisite  input  variables  describing  such  a 
scenario.  The  actual  target  and  background  reflectivities  and  target 
size,  if  known,  car  be  entered  by  the  program  user. 

Most  of  the  remaining  parameters  required  by  the  Obscuration 
Model  are  automatically  obtained  from  data  tables  in  the  mass  storage 
fij.c  of  AMSORB.  Attenuation  coefficients  for  various  types  of  obscur¬ 
ants  and  for  five  wavelengths  or  wavelength  ranges  are  given  In  Table  4-11, 
as  well  as  the  type  of  experimental  conditions  (atmospheric  or  smoke 
chamber)  under  which  the  data  were  obtained  and  the  reference;  citing 
the  experiments.  The  routines  automatically  select  the  proper  attenu¬ 
ation  coefficient  based  on  the  obscurant  type  specified  by  the  user  (see 
Section  2.4). 

The  Obscuration  Model  partitions  the  sky,  as  noted  in  Section 
2.4,  into  either  6  or  16  partitions  depending  on  whether  the  sun  is 
behind,  or  in  front  of,  the  target.  For  overcast  skies  (total  sky  cover 
greater  thau  or  equal  to  0.85),  the  routine  uses  16  sky  partitions. 

Table  4-12  shows  the  azimuth  and  zenith  angles  used  to  represent  the 
centroids  of  the  sky  partitions  for  the  brightness  calculations.  The 
azimuth  angles  shown  in  Table  4-12  are  based  on  a  reference  coordinate 
system  in  which  0  degrees  corresponds  to  grid  north  and  90  degrees  to 
grid  east.  The  routine  automatically  translates  these  azimuth  angles  to 
a  reference  frame  in  which  a  ray  with  an  azimuth  of  0  degrees  is  oriented 
along  the  observer  to  target  line-of-sight.  The  sun's  position  and  the 
total  sky  cover  are  passed  to  the  Obscuration  Model  through  the  data 
base  aastageaent  system  from  the  Mixing-Layer  Analysis  Routine  (see 
Section  4,1).  The  sun's  position  is  also  used  to  determine  the  position 
of  the  sun's  partition  when  direct  radiation  is  used  in  the  calculations 
for  clear  to  partly-eloudy  days.  Scattering  angles  for  the  sun  and  sky 
partitions  are  calculated  under  the  assumption  that  rhe  ligh*  sources 
are  at  very  large  distances  from  the  target  and  observer.  As  shown  in 
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TABLE  4-11 


ATTENUATION  COEFFICIENTS  FOR  VARIOUS  TYPES 
OF  OBSCURANTS 


Type  of 
Obscurant 

Wavelength 

Range 

(pm) 

Coefficient  of 
Attenuation 
(m2  mg-1) 

Test 

Conditions 

Reference 

HC  Smoke 

.4  -  .7 

.0033 

Atmospheric 

DPG  (1978b) 

WP  Smoke 

.4  -  .7 

.0021 

Atmospheric 

DPG  (1978b) 

Fog  Oil  Smoke 

.4  -  .7 

.0029 

Smoke  chamber 

Allen  and  Simonson  (1970) 

FS  Smoke 

.4  -  .7 

.0047 

Smoke  chamber 

Allen  and  Simonson  (1970) 

Explosive  Munition 
Dust 

.4  -  .7 

.00032 

Atmospheric 

DPG  (1978c) 

Vehicle  Movement 
Dust 

.4  -  .7 

.00024 

Atmospheric 

DPG  (1978d) 

Muzzle  Blast 

Smoke  and  Dust 

.4  -  .7 

.00024 

© 

DPG  (1978d) 

Burning  Brush 

.4  -  .7 

.0042 

Atmospheric 

Turner,  e£  al. ,  (1979), 
and  Eccleston,  et  ai., 
(1974) 

Burning  Building 

.4  -  .7 

.0076 

Atmospheric 

Turner,  et  al,(1979) 

Burning  Vehicle 

.4  -  .7 

.0076 

© 

- 

HC  Smoke 

1.06 

,00098 

Atmospheric 

DPG  (1978b) 

WP  Smoke 

1.06 

.00059 

Atmospheric 

DPG  (1978b) 

Explosive  Munition 
Dust 

1.06 

.00035 

Atmospheric 

DPG  (1978c) 

Vehicle  Movement 
Dust 

1.06 

.00019 

Atmospheric 

DPG  (1978d) 

Muzzle  Blast 

Smoke  and  Dust 

1.06 

.00019 

© 

DPG  (1978d) 

HC  Smoke 

3.4 

.00011 

Atmospheric 

DPG  (1978b) 

WP  Smoke 

J.4 

.00023 

Atmospheric 

DPG  (1978b) 

Explosive  Munition 
Dust 

1  Vehicle  Movement 

3.4 

.00027 

Atmospheric 

DPC  (1978c) 

3.4 


Vehicle  Movement 
Duse 


. 00016 


Atmospheric 


DPG  (1978d) 


TABLE  4-11  (Continued) 


Type  of 
Obscurant 

Wavelength 

Range 

(pm) 

Coefficient  of 
Attenuation 
(m"  mg-1) 

Tes  C. 

Conditions 

Reference 

Muzzle  Blast 

Smoke  and  bust 

3.4 

.00016 

© 

DPG 

(1978d) 

HC  Smoke 

9.75 

.000044 

Atmospheric 

DPG 

(1978b) 

WP  Smoke 

9.75 

.00027 

Atmospheric 

DPG 

(1978b) 

Explosive  Munition 
Dust 

9.75 

.00021 

Atmospheric 

DPG 

(1978c) 

Vehicle  Movement 
Dust 

9.75 

. 00013 

Atmospheric 

DPG 

(1978d) 

Muzzle  Blast 

Smoke  and  Dust 

9.75 

.00013 

© 

DPG 

(197  8d ) 

HC  Smoke 

10.6 

.00015 

Atmospheric 

DPG 

(1978b) 

WP  Smoke 

10.6 

.00048 

Atmospheric 

DPG 

(1978b) 

Fog  Oil  Smoke 

10.6 

.00025 

Smoke  Chamber 

Allen 

and 

Simonson 

(1970) 

FS  Smoke 

10.6 

.00052 

Sroke  Chamber 

Allen 

and 

Simonson 

(1970) 

Assumed  to  be  identical  to  the  attenuation  coefficient  measured  for 
dust  clouds  produced  by  surface  dust  generated  by  a  moving  vehicle. 

The  attenuation  coefficient  for  a  burning  vehicle  is  assumed  to  be 
approximately  equal  to  the  coefficient  for  a  burning  building. 


© 

© 
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TABLE  -4-12 

ANGLES  REPRESENTING  THE  CENTROIDS  OF 
THE  SKY  PARTITIONS 


£ 


\3 


(a)  6  Sky  Partitions  (used  when  the  sun  is  in  front  of  the  target) 

Azimuth 

Zenith 

4>  (deg) 

8  (deg) 

60.0 

14.4775 

180.0 

14.4775 

300.0 

14.4775 

60.0 

48.5904 

180.0 

48,5904 

300.0 

48.5904 

(b)  16  Sky  Partitions  (used  when 

the  sun  is  behind  the  target  or 

obscured  by  clouds) 

Azimuth 

• 

Zenith 

ip  (deg)- 

0  (deg) 

22.5 

14.4775 

67.5 

14.4775 

112.5 

14.4775 

157.5 

14.4775 

202.5 

14.4775 

247.5 

14.4775 

292.5 

14.4775 

337.5 

14.4775 

22.5 

48.5904 

67.5 

48.5904 

112.5 

48.5904 

157.5 

48.5904 

202.5 

48.5904 

247.5 

48.5904 

292.5 

48.5904 

337.5 

48.5904 

Figure  4-5,  the  light  rays  arriving  from  each  source  can  thus  be  assumed 

th 

parallel,  and  the  scattering  angle  <jn  for  the  i  light  source  is 
determined  from  the  cosine  law  as  follows: 

cos  <j>^  =  sin  8*  ♦  sin  0^  +  cos  0'  •  cos  0^  •  cos  ^  (4-11) 

where 

<jj^  =  azimuth  angle  with  respect  to  the  observer-target  line-of- 
sight  for  the  i^  light  source 

til 

0^  =  elevation  angle  of  the  i  light  source 

0*  =  elevation  of  the  line-of-sight  with  respect  to  the  plane 

of  the  horizon 

Table  4-13  gives  two  sets  of  the  Mie  scattering  function  F{<jn,  X} 

divided  by  the  square  of  the  propagation  constant  k  for  scattering  angle 

increments  of  2  degrees  currently  contained  in  the  data  base.  The 

values  for  the  dilute  WP  in  the  visible  spectrum  shown  in  Table  4-13  were 

obtained  from  the  SOMI  program  (Johnson,  et_  al ,  1972)  and  the  values  for 

HC  at  the  10.6-micrometer  wavelength  were  obtained  from  the  ASLSOM  code 

(Gomez,  1978).  The  obscuration  routines  use  linear  interpolation  to 

2 

determine  the  exact  value  of  F{<j>^,X}/k"  for  the  scattering  angle  calcu¬ 
lated  from  Equation  (4-11) .  In  the  absence  of  values  for  HC  smoke  in 
the  visible,  the  program  uses  the  values  shown  in  Table  4-13  for  dilute 
WP  smoke.  Should  Mie  fractions  for  other  obscurants  and/or  other  wave¬ 
lengths  be  required,  they  must  be  entered  by  the  user  (see  Section  A. 4 
of  Volume  II) . 

The  Obscurant  Dispersion  Model  routine  is  called  when  concen¬ 
trations  or  line-of-sight  integrated  concentrations  are  to  be  calculated 
for  use  in  the  transmission  and  brightness  calculations.  At  this  point, 
the  obscuration  routine  also  performs  the  tests  described  in  Section  2.4 
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TABLE  4-13 

VALUES  OF  F{4)it A>/k2  FOR  DILUTE  WHITE  PHOSPHORUS 
AND  HEXACHLOROETHANE  SMOKES 


Scattering 

Angle 

Scattering 

Fraction 

Scattering 
.  Angle 

Scattering 

Fraction 

Scattering 

Angle 

Scattering 

Fraction 

(deg) 

(deg) 

(deg) 

(a)  Dilute 

White  Phosphorous  (WP)  in  the  visible 

0.0 

1.17889x10"* 

62.0 

1.32288x10? 

124.0 

4.19813x10"^ 

2.0 

1.15003x10 

64.0 

1. 17992x10  ? 

126.0 

4.42084x10  ^ 

4.0 

1.06873x10  Z. 

66.0 

1.05331x10  c 

128.0 

4.65831x10 

6.0 

9.49033x10"^ 

68.0 

9.41838x10  , 

130.0 

4.91241x10“^ 

8.0 

8.09441x10 

70.0 

8.43939x10 

132.0 

5.18116x10 

10.0 

6.67796x10 

72.0 

7.57852x10  j? 

134.0 

5.45844x10"^ 

12.0 

5.37472x10  * 

74.0 

6.82063x10  ^ 

136.0 

5.73895x10 

14.0 

4.25876x10  ^ 

76.0 

6.15570x10  , 

138.0 

6.02246x10  1 

16.0 

3.35111x10  ^ 

78.0 

5.57700x10  ^ 

140.0 

6.31207x10"^ 

18.0 

2.63791x10  * 

80.0 

5.07699x10 

142.0 

6.60882x10“^ 

20.0 

2.08918x10  ^ 

82.0 

4.64580x10 

144.0 

6.91061x10  ° 

22.0 

1.67184x10  * 

84.0 

4.27384x10 

146.0 

7.21815x10  : 

24.0 

26.0 

1.35602x10  * 
1.11700x10  , 
9.35254x10"? 

86.0 

88.0 

3.95518x10 

3.68754x10 

148.0 

150.0 

7.53955x10  j? 
7.88504x10 

28.0 

90.0 

3.46875x10  ^ 

152.0 

8.25509x10  Z 

30.0 

7.95836x10  , 

92.0 

3.29384x10 

154.0 

8.63627x10  ^ 

32.0 

6. 87646xl0_? 

94.0 

3.15606x10 

156.0 

9.01061x10  * 

34.0 

6 . 02473x10? 

96.0 

3.05042x10 

158.0 

9.36426x10 

36.0 

5.34125x10  7 

98.0 

2.97533x10 

160.0 

9.67548x10"^ 

38.0 

4.77872x10  7 

100.0 

2.93072x10 

162.0 

9.88716x10 

40.0 

4.30176xl0_7 

102.0 

2.91520x10 

164.0 

9.90212x10  j! 

42.0 

44.0 

3.88561x10  7 
3.51418x10^7 

104.0 

106.0 

2.92579x10  : 
2.96015x10  J 

166.0 

168.0 

9.63625x10  * 
9.10971x10“^ 

46.0 

3. 17742x10? 

108.0 

3.01796x10 

170.0 

8.49746x10  ^ 

48.0 

2.86893xl0_? 

110.0 

3.09933x10  j? 

172.0 

8.06542x10  ;? 

50.0 

2 . 58465xlO_? 

112.0 

3.20259x10  , 

174.0 

8.00654x10  ^ 

52.0 

2.32238x10  7 

114.0 

3.32457x10 

176.0 

8.29279x10  ^ 

54.0 

2 . 08134xl0_? 

116.0 

3.46332x10 

178.0 

8.67358x10 

56.0 

1 . 86144x10? 

118.0 

3.61965x10 

180.0 

8.84540x10 

58.0 

1 . 66238x10? 

120.0 

3.79510x10  ^ 
3.98892x10  3 

60.0 

1.48329x10 

122.0 
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TABLE  4-13  (Continued) 

VALUES  OF  F { '4» £ ,  A  }  / k 2  FOR  DILUTE  WHITE  PHOSPHORUS 
AND  HEXACHLOROETHANK  SMOKES 


Scattering 

Angle 

(deg) 


Scattering 

Fraction 


Scatter  ing  | 
.Angle  i 
(deg)  [ 


Scatter  ing 
Fraction 


Scat  t  or ing 

Angl  e 
(deg) 


Scattering 

Fraction 


(b)  llexachloroethane  (110)  in  the  Infrared  (10.6  )tm) 


1.45197x10 

1.44700x10' 

1.43221x10' 

1.40794x10' 

1. 37476x10' 

1.33343x10' 

1.28485x10' 

1.23007x10" 

1.17022x10' 

1 . 1 0648x10 

1.04004x10' 

9.72039x10' 

9.03570x10' 

8.35627x10' 

7.69090x10' 

7.04710x10" 

6.43100x10' 

5.84734x10' 

5.29951x10' 

4.78963x10' 

4.  31870x1  O' 

3.88673x10' 

3.49290x10' 

3. 13576x10' 

2.81337x10' 

2. 52346x10“ 

2.26356x10' 

2.0311 3x10' 

1.82363x10' 

1.638  3x10' 

1.47382x10' 


62.0  1.32705x10  5 

64.0  1.19637x10.3 

66.0  1. 08002x1 0_6 

68.0  9.76416x10_£ 

70.0  8.841 54x1 0^ 

72.0  8 . 01 987x1 0l6 

74.0  7.28812xl0_f) 

76.0  6. 63652x1 0_^ 

78.0  6.05637x10^3 

80. 0  5.53996x10_£ 

82.0  5. 08040x1 0_£ 

84.0  4.671  52x1 0_c1 

86.0  4.30779x10_6 

88.0  3. 98427x 10_6 

90.0  3 . 696  52x1 0_  ^ 

92.’o  3.44060x10 _ 

94,0  3.21297x10^ 

96.0  3. 01051x10 

98.0  2. 83046x1  O'’ 

100.0  2 . 67037x1  O*1 

102.0  2.52811x10 ' 

104.0  2.40180x10  ’ 

106.0  2 . 28981x10”° 

108.0  2.19069x10  (' 

110.0  2.1031 8x1 0~h 

112.0  2 . 026 15x1 0-6 

114.0  1 . 95864x1 0“6 

116.0  l .  89975x10“*’ 

118.0  1 . 84872x1 0“6 

120.0  1 . 80484x1 O-^ 

122.0  1 . 76752x1 0-6 


124.0 

126.0 

128.0 

130.0 

132.0 

134.0 

136.0 

138.0 

140.0 

142.0 

144.0 

146.0 

148.0 

150.0 

152.0 

154.0 

156.0 

158.0 

160.0 

162.0 

164.0 

166.0 

168.0 

170.0 

172.0 

174.0 

176.0 

178.0 

180.0 


1 .73621x1015 
1.71 047x1 0_6 
1 .68992xl0_6 
1.67422x10  , 
1.66315x10  , 
1.65651x10“^ 
1 .6541  7xloI6 
1 .65602x10_& 
1 . 66200x10_6 
1 .67204x10_6 
1 .68605x10_6 
1 . 70393xl0_6 
1  .  7255 1 xlO_6 
1 .  75057x10_6 
1.77878x10_6 
1 .80974x10_6 
1.84294x10_6 
1.87779x10_6 
1 . 9 1 359x10_6 
1. 94956x10. 6 
1.98490x10_6 
2 . 01 874x1 0_6 
2 . 05022x1 0_8 
2.07852xl0_6 
2. 1 028  5x1 0_6 
2. 12255xl0„6 
2 . 1  3704x 10_6 
2 . 1 4  59 1 x 1 0_6 
2.14889x10 


to  determine  if  any  obscurant  is  required  to  obscure  the  target  or  if  there 
is  too  little  obscurant  present  to  affect  the  probability  of  detection. 

The  amount  of  direct  and  diffuse  illumination  of  the  scenario  is  an  important 
parameter  in  the  calculation  of  the  target  and  background  brightness,  cloud 
brightness  and  in  determining  the  threshold  contrast.  The  present  version 
of  the  obscuration  routines  provides  for  three  levels  of  sky  luminance:  a 
clear  to  partly  cloudy  day  (total  cloud  cover  less  than  0.85)  with  the  sun 
near  its  meridian,  a  mostly  overcast  day  (total  sky  cover  equal  to  or 
greater  than  0.85)  and  a  clear  to  partly-cloudy  day  at  sunrise  or  sunset. 
Luminance  levels  for  these  sky  conditions  used  in  the  routines  are  shown  in 
Table  4—1 4  (Handbook  of  Chemistry  and  Physics,  1957).  The  direct  sun 
luminance  of  6x10^  candles  per  square  meter  for  sunrise  and  sunset  under 
clear  to  partly  cloudy  sky  conditions  is  considered  to  occur,  for  calcula¬ 
tion  purposes,  between  one  half-hour  prior  to  and  one  half-hour  after 
sunset  and  sunrise,  respectively.  The  diffuse  sky  luminance  of  8000 
candles  per  square  meter  is  an  average  brightness  for  the  sky  according  to 
the  Handbook  of  Chemistry  and  Physics,  (1957).  The  brightness  for  each  sky 
partition  B^  is  obtained  by  dividing  the  indirect  radiation  value  given  in 
Table  4-14  by  6  or  16,  depending  on  the  number  of  sky  partitions  used  in 

the  calculation.  For  clear  to  partly-cloudy  skies,  when  direct  radiation 

9 

is  present,  the  sun's  partition  is  assigned  the  brightness  of  1.6x10 
candles  per  square  meter.  For  infrared  scenarios,  the  sky  brightness  is 
not  input.  Instead,  as  described  in  Section  A. 4  of  Volume  II,  radiance 
values  for  the  target,  sky  and  terrain  in  microwatts  per  square  meter  per 
steradian  must  be  input  by  the  user. 

The  threshold  contrasts  for  a  50-percent  probability  of  detection 
used  in  the  routines,  and  shown  in  Table  4-15  as  a  function  of  illumination 
level  and  target  size,  are  compiled  from  data  presented  by  Blackwell 
(1946).  It  should  be  noted  that  these  threshold  contrasts  are  for  viewing 
in  the  visible  spectrum  by  the  human  eye.  Two  sets  of  threshold  contrast 
values  are  given,  the  first  for  use  when  the  brightness  of  the  target  has 
been  calculated  to  be  more  than  the  background  brightness  and  the  second 
set  when  the  target  brightness  is  less  than  the  background  brightness.  In 
the  test  to  determine  if  smoke  is  required  to  obscure  the  target,  the 
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TAPl'.E  4-14 

LUMINANCE  LEVELS  FOR  VARIOUS  SKY  CONDITIONS 


1  I 


1 


Direct 

Radiat  ton^ 
(caiulles/n  ) 

Indirect 

Radiat ion^ 
(candles/m ') 

Skv 

Condition 

l.6xl09 

8000 

Clear  to  Partly  Cloudy 
(TSC*  <  0.85) 

0 

8000 

Mostly  Overcast 
(TSC  >  0.85) 

6.0xl06 

_ _ _ 

8000 

Clear  to  Partly  Cloudy 
Sunrise  &  Sunset 
(TSC  <  0.85) 

if 

TSC  is  the  total  sky  cover  or  fraction  of  the  sky  covered  by  clouds. 


»l 


*1 


i 

i 

t 
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THRESHOLD  CONTRAST  VALUES 


Target  Size 

Total  Available 

2 

Luminance  (candles/m  ) 

(minutes  of  arc) 

>104 

103 

9 

10“ 

10 

1 

(a)  Threshold 

Contrasts 

for  Targets 

Brighter  than  the  Background 

150 

0.0027 

0.0027 

0.0027 

0.0030 

0.0038 

100 

0.0027 

0.0027 

0.0028 

0.0032 

0.0044 

50 

0.0028 

0.0028 

0.0030 

0.0034 

0.0050 

20 

0.0031 

0.0033 

0.0036 

0.0042 

0.0073 

l  5 

0.0034 

0.0038 

0.0044 

0.0051 

0.0095 

10 

0.0042 

0.0050 

0.0064 

0.0073 

0.0150 

8 

0.0052 

0.0062 

0.00?' 

0.0096 

0.020 

6 

0.0068 

0.0083 

0.012 

0.014 

0.030 

4 

0.011 

0.0135 

0.021 

0.025 

0.055 

2 

0.025 

0.036 

0.062 

0.082 

0.190 

1 

0.072 

0.130 

0.230 

0.340 

0.800 

(b)  Threshold 

Contrasts 

for  Targets 

Darker  than 

the  Background 

150 

0.0070 

0.0070 

0.0072 

0.0073 

0.0092 

100 

0.0074 

0.0074 

0.0077 

0.0079 

0.0102 

50 

0.0037 

0.0087 

0.0094 

0.0100 

0.0138 

20 

0.0130 

0.0130 

0.014 

0.0160 

0.025 

15 

0.0155 

0.0155 

0.017 

0.0195 

0.031 

10 

0.021 

0.021 

0.023 

0.027 

0.045 

8 

0.026 

0.026 

0.029 

0.034 

0.059 

6 

0.035 

0.034 

0.040 

0.047 

0.085 

4 

0.056 

0.054 

0.065 

0.082 

0.150 

2 

0.  150 

0.140 

0.175 

0.320 

0.410 

1 

0.470 

0.450 

0.520 

0.660 

1.00 

target  and  background  reflectivities  are  compared  to  determine  which  set  to 
use  in  selecting  the  threshold  contrast  for  comparison  with  the  contrast 
calculated  from  Equation  (2-150).  The  angular  size  of  rhe  target  in 
minutes  of  arc  6  is  determined  from  the  expression 

8  =  120  tan'1  (d/2D)  (4-12) 

where  d  is  the  diameter  of  the  target  obtained  either  from  Table  4-13 
when  a  code  is  used  to  specify  the  scenario,  or  from  input  specified  by 
the  user.  The  distance  D  between  the  target  and  observer  is  given  by 


For  target  sizes  intermediate  between  the  tabulated  values  shown  in 
Table  4-15,  the  routines  use  linear  interpolation  to  obtain  the  requi¬ 
site  threshold  contrasts.  Although  Blackwell's  data  show  that  the  thres¬ 
hold  contrast  is  not  linear  with  target  size,  the  errors  introduced  by 
linear  interpolation  between  the  target  sizes  in  Table  4-.\5  do  not 
greatly  affect  the  probability-of-detection  calculation. 

Finally,  as  shown  in  Figure  4-4,  the  cloud  brightness  is  cal¬ 
culated,  if  necessary,  and  the  probability  of  detection  is  calculated 
for  the  first  time  step.  After  the  time  profile  of  probability-of- 
detection  has  been  completed,  the  results  are  stored  in  AMSORB's  data 
files. 

4.5  BATTLEFIELD  ENVIRONMENT  ROUTINE 

The  first  function  of  the  AMSORB  Battlefield  Environment 
Routine  (BEC)  is  to  determine  sources  on  the  battlefield  that  most 


likely  contribute  to  the  observer-target  LOS  concentration,  and  to  the 
obscuration  of  the  target.  BEC  then  references  the  Trajectory/Transport 
Routine,  Obscurant  Dispersion  Model  and  Obscuration  Model  Routine  to 
perform  the  concentration  and  obscuration  calculations  for  these  selected 
battlefield  sources.  Figure  4-6  is  a  schematic  diagram  showing  the  flow 
logic  of  BEC. 

As  shown  in  the  first  two  boxes  of  the  BEC  schematic  diagram 
(Figure  4-6),  the  routine  reads  user-provided  input  data  discussed  in 
Section  A. 4  of  Volume  II  and  reads  the  raesoscale  analysis  meteorological 
data  base  created  from  the  Analysis  Phase  routines  of  AMSORB.  From 
these  input  data,  the  routine  computes  parameters  remaining  constant 
throughout  the  calculations.  Due  to  the  si2e  of  the  battlefield  grid 
and  the  need  to  conserve  calculation  time,  a  closed  boundary  is  defined 
using  the  Trajectory/Transport  Routine  discussed  in  Section  2.2  to 
eliminate  sources  which  can  not  possibly  affect  the  observer-target  LOS 
concentration  or  probability  of  detection  calculation.  As  shown  on  the 
second  page  of  figure  4-6,  BKC  loops  over  all  battlefield  source  loca¬ 
tions  to  determine  sources  within  the  closed  boundary,  reading  the 
source  locations  from  the  battlefield  source  location  and  characteris¬ 
tics  data  base  previously  created  by  the  SOKDAT  Routine  discussed  in 
Section  4.2.  The  BEC  Routine  sets  a  flag  (eligible  source  flag)  for 
those  sources  within  the  boundary,  effectively  eliminating  sources 
outside  the  boundary  from  further  consideration,  and  references  the 
Tcajectory/Transpor t  Routine  to  calculate  an  average  wind  speed,  travel 
distance  and  travel  time  between  the  grid  point  closest  to  each  source 
and  the  LOS. 


As  shown  on  pages  3  and  4  of  Figure  4-6,  BEC  makes  a  second 
loop  over  all  battlefield  sources  and  again  references  the  battlefield 
source  location  and  characteristics  data  base  created  by  SORDAT.  Al¬ 
though  processing  battlefield  source  data  would  be  more  efficient  if  the 
data  were  initially  transferred  from  the  mass  storage  data  base  to  core 
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memory,  BFC  assumes  that  the  initial  number  ot'  battlefield  sources  can 
exceed  BF.C's  core  memory  array  allocation.  The  purpose  of  this  second 
and  successive  loops  over  the  battlefield  source  data  from  a  mass  storage 
data  base  is  to  select  the  sources  most  likely  to  contribute  to  the  LOS 
concentration  and  the  obscuration  of  the  target.  For  "eligible"  sources, 
BEC  determines  if  the  source  has  been  released  prior  to  the  user- Input 
calculation  time.  If  Che  source  has  been  released,  BEC  determines  if 
the  source  is  within  one  mesoscale  grid  spacing  from  the  line-of-sight . 

In  our  diagram,  we  have  called  this  source  a  "close"  source  and  all 
other  more  distant  sources  are  called  "fir"  sources.  If  the  source  is 
"close",  BEC  references  cite  appropriate  obscurant  dispersion  model  and 
computes  its  contribution  to  tiie  observer-target  LOS  integrated  concen¬ 
tration.  The  computed  concentration  is  then  summed  with  the  contribu¬ 
tions  from  other  "close"  sources.  If  BEC  determines  that  the  summed 
contributions  are  sufficient  to  obscure  the  target,  all  processing  is 
terminated  and  the  user  is  notified  that  the  target  is  obscured. 

If  no  obscuration  exists  at  this  point,  BEC  has  recognized  all 
sources  which  have  been  released  prior  to  the  user-specified  calculation 
time  and  has  calculated  LOS  integrated  concentrations  for  "close"  sources. 
The  BEC  logic  shown  on  page  5  of  Figure  4-6  is  designed  to  determine  if 
the  clouds  from  the  "far"  sources  are  in  the  vicinity  of  the  LOS  at  the 
specified  time,  defines  the  "close"  sources  and  the  "far"  sources  whose 
clouds  are  In  the  vicinity  of  the  LOS  at  the  user  specified  time  as 
"active"  sources,  and  assigns  a  priority  to  these  active  sources  based 
on  distance  from  the  LOS  ("far"  sources)  arid  concentration  ("close" 
sources).  l‘p  to  this  point,  BEC  lias  retrieved  and  stored  all  battle¬ 
field  source  data  l rum  and  to  mass  storage  data  bases.  To  increase  the 
speed  of  subsequent  processing,  BEC  stores  all  information  for  the 
active  sources  with  highest  priority  into  the  program's  core  memory 
arrays  as  stiuwn  on  page  b  of  Figure  4-6.  It  should  be  noted  that  the 
number  of  active  sources  could  exceed  the  routine's  core  memory  array 
limits  (currently  equal  to  100).  If  this  occurs,  the  routine  prints  a 
'.'anting  message  and  processes  the  700  highest  priority  sourecs. 
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The  primary  purpose  of  the  BEC  logic  shown  on  page  7  of 
Figure  4-6  is  to  effectively  u3e  the  Trajectory /Transport  Routine  de¬ 
scribed  in  Section  2.2  and  the  Dispersion  Models  described  in  Section 
2.3.1  in  moving  and  dispersing  clouds  from  a  point  of  origin  "far"  from 
the  DOS  to  a  distance  not  greater  than  one  mesoscale  grid  spacing  from 
the  I, OS.  In  effect,  the  models  are  used  to  redefine  a  "new"  source 
close  to  the  I.OS  with  dimensions  equivalent  to  those  predicted  to  occur 
for  cloud  travel  over  a  trajectory  from  the  source  far  from  the  I.OS. 

The  "new"  source  has  a  release  time  equivalent  to  the  original  time  of 
release  plus  the  predicted  travel  time  along  the  trajectory  from  the 
point  of  origin  to  the  location  of  the  "new"  source.  The  appropriate 
cloud-rise  equations  described  in  Section  2.3.1  are  also  used  to  define 
the  cloud  rise  from  buoyant  sources  and  thus  the  effective  height  of  the 
"new"  source.  As  shown  on  page  H  of  Figure  4-6,  the  contribution  of  the 
"new"  source  to  the  l.OS  integrated  concentration  is  then  calculated  and 
added  to  the  sum  of  concentrations  from  the  sources  originally  "close" 
to  die  line  of  sight  an:!  all  previously  processed  "new"  sources  whose 
origin  was  closer  in  distance  to  the  LOS.  If  the  summed  concentration 
clearly  exceeds  the  concentration  required  to  form  an  obscuration,  the 
progroa  ceases  to  process  sources  and  informs  the  user  the  target  is 
obscured.  Otherwise,  sources  are  processed  until  there  are  no  more 
active  sources  to  consider.  At  this  point,  the  program  calculates  the 
probability  »  !'  detecting  the  target,  if  the  obscuration  option  has  been 
selected  by  the  user,  using  the  Obscuration  Model  Routine  described  in 
Section  2.4.  the  Obscuration  Model  Routine  also  eotnpuf'  ?i  the  target, 
background  and  cloud  brightness,  and  the  target-back).- r  und  contrast. 

When  all  model  processing  is  completed,  tLKC  -utputs  the  input 
data  and  the  results  tu  a  printer  device  and  to  a  h-t flctield  results 
data  base  as  shown  on  the  final  page  of  Figure  4-T.  §KC:  outputs  to  the 
printer  device  some  ot  the  user-provided  input  d.ata.  Che  observer-target 
Hne-of-sight  Integrated  concentration  and  the  observer- tar get  trans¬ 
mittance.  For  the  obscuration  option,  the  fltoseura .  lorn  Model  calculation 
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results  are  also  printed.  Finally,  a  brief  message  is  listed  indicating 
what  portions  of  BEC  were  utilized.  The  messages  may  also  indicate 
diagnostics  or  errors  that  occurred  during  program  execution.  Refer  to 
Section  6  for  an  example  print  output  produced  by  3EC.  Upon  completion 
of  all  model  processing,  the  battlefield  results  data  base  contains  all 
user-provided  data  and  meteorological  constants  used  in  the  data  case 
and  all  calculated  results.  The  data  base  also  contains  a  calculation 
result  code  which  indicates  t he  portions  of  UEC  utilized  in  the  proces¬ 
sing  and  is  analogous  to  the  messages  listed  on  the  printer  device. 


HI 


SECTION  5 

STRUCTURE  OK  THE  COMPUTER  PROGRAM: 

ANALYSIS,  DATA  PROCESSOR,  APPLICATION  AND 
GRAPHICS  PHASES  OF  AMS ORB 


Figure  5-1  shows  the  executive  components,  in  the  form  of  a 
block  diagram,  of  the  AMSORB  system.  The  AMSORB  executive  system 
manager,  labeled  XKCAMS  in  Figure  5-1,  is  where  all  processing  begins 
and  ends.  The  executive  systems  manager  on  option  passes  control 
to  the  Supervisors'  labeled  ANALIZ,  DATPRO,  APPLY  and  PLTPRG. 

The  Analysis  Supervisor  ANALIZ  controls  all  processing  in 
the  Analysis  Phase  of  AMSORB  and  passes  control  to  the  Mixing-Layer 
Analysis  and  Mesoscale  Wind-Field  Models  foreman  MIXLYR  or  to  other 
meteorological  analysis  model  routines.  The  operational  elements  of 
MIXLYR  are  described  in  Section  5.1  below. 

The  Data  Processor  Supervisor  DATPRO  controls  all  processing 
in  the  Data  Processor  Phase  of  AMSORB  and  passes  control  to  the  Battle¬ 
field  Source  Characteristics  Routine  foreman  SORUAT.  Section  5.2 
contains  a  disc-.sion  of  the  operational  elements  of  SORDAT.  The  Appli¬ 
cation  Supervisor  APPLY  passes  control  on  option  to  the  Trajectory/Transport 
Dispersion  Model  -pplication  foreman  TRY DIF,  to  'he  MS3  Obscurant  Dis¬ 
persion  and  Obscuration  Models  application  foreman  MS3,  to  the  Battle¬ 
field  Environment  Routine  (BEG)  Obscurant  Dispersion  and  Obscuration 
Models  application  foreman  BEC,  or  to  ocher  application  routines. 

The  operational  elements  of  TRYDtF,  MS3  and  BEC  are  described  in 
Section  5.3.1,  5.1.2  and  5.3.3,  respectively.  At  the  completion  of  any 
of  these  phases,  program  eoncrol  is  returned  to  XKCAMS.  The  executive 
systems  manager  can  address,  on  option,  the  Graphics  Supervisor 
PLTPRG  to  obtain  either  printed  or  plotted  solutions  developed  in  the 
application  routines.  Further  details  of  the  Graphics  Supervisor 
PLTPRG  are  described  in  Section  ‘.4  below. 
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5.1 


ANALYSIS  PHASE  OF  AMSORB 


As  shown  in  Figure  5-1,  the  foreman  M1X1.YR  of  the  Analysis 
Phase  of  AMSORB  is  comprised  of  the  two  operational  elements  MXOBAN  and 
SOLVMX.  The  operational  subroutines  of  MXOBAN  are  shown  in  the  block 
diagram  In  Figure  5-2.  The  operational  element  MXOBAN  controls  the 
analysts  of  all  meteorological  data  and  passes  the  results  of  the 
data  analysis  to  the  subroutine  SOLVMX  at  the  completion  of  execution. 

The  AMSORB  System  user  requests  transport  and  dispersion  information  for 
a  given  time  and  source,  which  is  then  passed  to  MXOBAN  via  XECAMS .  The 
MXOBAN  Mixing- Layer  Analysis  Model  routine  then  proceeds  to  analyze 
rawLnsonde,  CWC  and  surface  meteorological  data  according  to  the  logic 
process  described  in  Section  4.1.  The  subroutine  MXUADT  obtains  ra- 
winsonde  data  from  the  AMSORB  data  base  for  the  desired  solution  time  and 
source  location.  CWC  prediction  data  and  surface  weather  observations 
for  use  in  the  analysis  are  obtained  from  the  data  base  via  the  sub¬ 
routines  MXCWDT  and  MXSFOT.  These  subroutines  also  use  the  CWC  and 
surface  data  to  modify  the  raulnsonde  data  if  appropriate.  The  sub¬ 
routine  MXSFOT  also  calculates  the  net-radiation  index  for  the  surface 
meteorological  station  using  information  obtained  from  the  LONLAf  and 
DAYNIT  subroutines.  The  rawinsonde  data  are  then  used  in  the  subroutine 
LYKDHM  to  calculate  the  depth  or  the  surface  mixing  layer  and  the  mean 
wind  components  in  the  layer  as  described  in  Section  4.1.  If  insuffic¬ 
ient  data  are  available  for  calculating  the  mixing-layer  depth  or  wind 
components,  default  data  are  supplied  via  the  LYRDHM  and  MXSFOT  sub¬ 
routines.  Control  is  then  returned  to  the  operational  element  MIXLYK 
and  processing  continues  by  branching  to  the  subroutine  SOLVMX. 

A  block  diagram  showing  operational  subclcaenta  of  SOLVMX  is 
given  in  Figure  5-i.  Processing  begins  through  examination  of  the  Meso- 
scale  Wind-Field  Mrdel  routine  input  parameters  passed  to  SOLVMX  from 
MIXLYK  to  determine  which  of  two  procedures  are  to  be  u- e J  in  initial- 
icing  t tie  Mesoscale  Wind-Field  Model  routine.  Thus,  the  input  param¬ 
eters  from  MtXLYR  are  eoapared  in  the  subroutine  MX S TXT  with  the  input 
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■Sure  5-2.  Operational  Elements  of  1MX0BAN  Analysis  Phase 


F;<*ure  5-3.  Operational  Elements  of  S0LVMX  Analysis  Phase. 


parameters  used  to  generate  the  last  solution  of  the  mesoscale  routine 
resident  in  the  AMSORB  data  base.  The  mesoscale  routine  is  initialized 
with  the  new  inputs  from  MXOKAN  if  the  mean  layer  wind  speeds  and  ele¬ 
vation  of  the  mixing  depth  above  mean  sea  level  (MSI.)  differ  by  more 
than  twenty-five  percent  and  the  mean  wind  directions  differ  by  more 
than  ten  degrees  (cold  start).  The  initial  layer  depth  and  momentum 
components  for  a  cold  start  are  given  by  the  expressions 
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where 

=  mixing  depth  (MSL)  calculated  for  the  present  hour 

U  =  easterly  component  of  the  mean  layer  wind  for  the 
present  hour 

Miin  =  rainiraum  terrain  elevation  (MSL)  in  the  solution  matrix 

=  northerly  component  of  the  mean  layer  wind  for  the  pre¬ 
sent  hour 

If  the  differences  are  less  than  those  specified  above,  subroutine  MXSTRT 
calculates  initialization  parameters  for  a  "hot"  start  from  the  expres¬ 
sions 
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where 


mixing  uepth  (MSI)  at  each  solution  matrix  grid  point 
for  the  previous  wind-field  solution 


terrain  elevation  (MSL)  at  each  solution  matrix  grid 
point 


H'  =  mixing  depth  (MSL)  used  to  initialize  the  previous 
wind- field  solution 


easterly  component  of  the  mean  layer  wind  at  each  solu¬ 
tion  matrix  grid  point  for  the  previous  wind-field  solu¬ 
tion 


U  =  easterly  component  of  the  mean  layer  wind  used  to  initial¬ 
ize  the  previous  wind-field  solution 
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V 


i.J 


northerly  component  of  the  mean  layer  wind  at  each  solu¬ 
tion  matrix  grid  point  for  the  previous  wind-field  solu¬ 
tion 


V 

o 


northerly  component  of  the  mean  layer  wind  used  to  init¬ 
ialize  the  previous  wind-field  solution 


After  the  initialization  procedure  has  been  completed,  the  input  param¬ 
eters  are  passed  to  the  subroutine  MXSOLV  via  SOLVMX  for  execution  of 
the  Mesoscale  Wind-Field  Model  described  in  Section  2.1. 


Subroutine  MXSOLV  will  allow  the  mesoscale  routine  to  execute 
for  a  maximum  model  time  of  approximately  7200  seconds  for  a  cold  start 
and  1800  seconds  for  a  hot  start  or  until  the  solution  converges,  which¬ 
ever  occurs  first.  The  subroutine  begins  to  check  convergence  criteria 
after  1200  seconds  model  time  has  elapsed  for  a  cold  start  or  600  seconds 
for  a  hot  start.  The  convergence  criteria  are  then  checked  after  every 
application  of  the  nine-point  filter  (after  every  seventh  time-step  If  the 
default  value  is  used).  Subroutine  MXSOLV  begins  the  convergence  check 
by  calculating  the  relative  variance  of  the  change  in  mixing  layer  depth 
V  (AH  )  over  the  solution  matrix  according  to  the  expression 


(5-7) 
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where 


°  layer  depth  calculated  for  the  present  time  step  at  each 
point  in  the  n  by  m  matrix 


n  layer  depth  for  the  previous  time  step  at  each  point  in 
the  n  hy  m  matrix 


The  subroutine  determines  the  slope  of  the  curve  describing  as 

a  function  of  time  and  the  difference  between  V  (AH  )  at  the  present 

K  TU 

and  previous  application  of  the  filter.  The  solution  is  assumed  to 
have  converged  if  the 


•  Slope  of  the  curve  describing  V^(AH  J  as  a  function  of 
time  is  equal  to  or  less  than  zero 

•  Present  value  of  V^(AH  )  is  less  than  the  previous 

value  and  the  difference  between  the  values  is  less 

-5 

than  or  equal  to  1  x  10 

The  Analysis  Phase  saves  the  converged  Mesoscale  Wind-Field  Model  solu¬ 
tion  and  other  parameters  calculated  in  the  Mixing-Layer  Analysis  Model 
routines  by  passing  them  to  the  AMSORB  data  base.  At  this  point  M1XLYR 
is  terminated  and  control  is  returned  to  XECAMS. 


5.2  DATA  PROCESSOR  PHASE  OF  AMSORB 

The  AMSORB  executive  XECAMS  accesses  the  Data  Processor  Phase 
by  calling  Supervisor  DATPRQ.  Supervisor  DATPRo  passes  program  control 
to  t».«  Battlefield  Source  Characteristics  Routine  foreman  SORDAT,  or 
passes  program  control  to  other  data  processor  foreman  routines. 
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Figure  5-4  shows  a  block  diagram  of  Che  Data  Processor  Phase 
of  AMSORB  after  program  control  is  passed  from  Supervisor  DATPRO.  Foreman 
SORDAT,  which  is  also  the  main  routine  of  the  source  characteristics  program 
reads  the  first  data  card  and  calls  one  of  the  following  source  type 
or  category  subroutines:  RDSMK ,  BRNBR ,  BRNBL,  EXPMUN,  MUZZLE,  or  VHMOVE, 
to  read  the  remaining  input  data  cards  for  the  source  category  specified 
in  the  first  data  card.  These  specific  subroutines  provide  default  values 
and/or  compute  tin*  source  input  parameters  particular  to  each  type 
source  category.  The  default  values  and  source  model  used  for  a  particu¬ 
lar  source  category  are  based  on  the  values  and  source  models  discussed 
in  the  corresponding  source  category  subsection  in  Section  3.  These 
routines  write  all  model  input  parameters  to  a  mass  storage  data  base 
that  is  in  turn  referenced  by  Che  Battlefield  Environment  Routine 
(BEC)  discussed  in  Section  5.3.3.  After  the  processing  of  data  for  one 
source  is  completed,  a  return  is  made  by  the  category  subroutine  to 
SORDAT  to  read  the  next  set  of  source  cards  until  an  end-of-file  card 
is  encountered. 


At  the  users  option,  SORDAT  calls  DRVPRT  to  print  source 
characteristics  tables  containing  input  values  and  model  input  parameters. 
DRVPRT  utilizes  the  following  subroutines  in  creating  these  tables: 

PKTSOR,  PRTSMK ,  PUTINS,  PRTQC,  PRTBBR,  PRTBBL,  PRTBVH ,  PRTEXP,  PRTMUZ 
and  PUTVHM.  Program  control  is  then  returned  from  SORDAT  to  the  Data 
Processor  Supervisor  DATPRO  which  in  turn  passes  control  to  XECAMS. 

5.3  APPLICATION  PHASE  OF  AMSORB 

5.3.1  Foreman  TRY DIF 

Inspection  of  Figure  5-1  above  allows  the  foreman  TRYDIF  of 
the  Applications  Phase  of  AMSORB  is  comprised  of  the  two  operational 
elements  TRNPRT  and  DtFFUZ.  These  operational  elements  perform  whe 
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Rout incB  Controlled  by  the  Data  Processor 


FIGURE  5-5.  Routines  controlled  by  the  Application  Phase  TRYDIF  Foreman  of  AMS ORB . 


calculations  outlined  in  the  block  diagram  shown  in  Figure  5-5  after 
receiving  program  control  through  the  foreman  TRYD1F  from  XECAMS  and 
Supervisor  APPLY. 

The  subroutine  TRNPRT  calculates  the  trajectory  of  the  cloud 
centroid  downwind  from  the  source  location  using  the  output  from  the 
Mesoscale  Wind-Field  Model  routine  and  the  procedures  outlined  in  Section 
2.2  above.  Subroutine  TRNPRT  calculates  a  point  on  the  trajectory  at 
intervals  equal  to  one-tenth  the  spacing  of  the  solution  matrix  of  the 
Mesoscale  Wind-Field  Model  routine  (default  value  is  500  meters).  At 
each  of  the  calculated  points  on  the  trajectory,  tne  subroutine  also 
determines  the  appropriate  mixing-layer  depth  and  wind  speed  for  use 
in  the  dispersion  models,  using  the  interpolation  techniques  described 
in  Section  2.2.  After  completion  of  the  calculations,  program  control 
is  passed  to  subroutine  L'lFFUZ  via  TRYD1F. 

The  subroutine  lHFFL'Z  calculates  concentration  and  dosage  with 
optional  cloud  depletion  due  to  decay  and  either  deposition  losses  from 
gru' i t a* i uua 1  settling  or  precipitation  scavenging,  using  the  Instan¬ 
taneous  and  continuous-source  dispersion  models  described  in  Section  2.3 
at.o'e.  At  ttie  user's  option,  the  program  also  calculates  concentration 

and/or  dosage  isopieth  half-widths  at  each  point  along  the  trajectory 
for  Instantaneous  sources  and  concentration  isopieth  half-widths  for 
continuous  sources.  At  the  completion  of  the  dispersion  model  calcu¬ 
lations,  the  program  writes  the  solution  to  the  AMSORB  print  file  and 
returns  control  to  XECAMS. 

5.3.2  Foreman  MS3 

Figure  5-1  shows  that  the  Foreman  MS3  of  the  applications 
phase  of  AMSORB  consists  of  the  two  primary  operational  elements  HECSOM 
and  Pl'OUT.  Figure  5-6  is  a  block  diagram  showing  the  function  of  these 
primary  operational  elements  after  receiving  program  control  through 


FIGURE  ,  Routines  controlled  by  the  Application  Phase  MS  3  foreman  of  AMSORB . 


foreman  MS3  from  XF.CAMS  and  the  Supervisor  APPLY.  Figure  5-6  also 
shows  the  relationship  between  these  primary  operational  elements  and 
subordinate  elements. 


Subroutine  HECSOM  supervises  the  number  of  times  to  be  processed 
in  the  time  profiles,  controls  concentration  calculations  and  contains 
the  MS3  Obscuration  Model  described  in  Section  2.4.1.  For  each  time  in 
the  time  profile,  this  subroutine  calculates  the  target,  background  and 
cloud  brightnesses  which  are  then  used  to  compute  the  time  profiles  of 
probabl 1 i ly-of-detect ion  and  transmittance.  This  subroutine  also 
passes  program  control  to  the  subordinate  routine  ALLCN  when  a  line- 
integrated  concentration  is  required,  or  when  the  number  of  smoke  cloud 
elements  and  the  concentration  at  each  element  are  required.  Program 
control  ii  passed  back  to  foreman  MS 3  after  completion  of  all  calcu¬ 
lations  performed  by  subroutine  HECSOM. 

The  MSI  Dispersion  Model  described  in  Section  2.3.2  is  con¬ 
tained  within  routines  ALLCN,  CH1PT  and  CHUN.  Subroutine  HECSOM  passes 
the  observer  and  target  locations  (as  end  points  of  the  line-of-sight) 
to  ALLCN.  ALLCN  controls  the  integration  of  concentration  along  the 
l iue-of-sight  between  the  observer  and  target  and  also  performs  the 
calculation  of  line  integrated  concentrations  between  light  sources  and 
the  target,  background  atiu  cloud  elements.  ALLCN  also  calculates  the 
number  and  position  of  cloud  elements  required  to  obtain  an  accurate 
description  of  the  cioo-j  for  use  in  the  eloud  brightness  calculations 
performed  by  the  MS 1  Obscuration  Model  resident  in  HECSOM.  When  con¬ 
centrations  at  points  along  the  observer-target  line-of-sight  arc  re¬ 
quired,  subroutine  ALLCN  passes  control  to  CHIPT  or  CHI  IN  where  the  concen¬ 
trations  at  the  points  are  calculated  using  Equation  (2-40)  or  (2-64)  after 
uhich  control  is  returned  to  ALLCN.  Program  control  is  passed  baek  to 
the  calling  routine  HECSOM  after  completion  of  the  calculations  performed 
in  routine  ALLCN. 
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Subroutine  PPOUT  controls  the  print  output  and  writes  all 
calculations  to  a  save  file  for  later  use  by  the  'Graphics  Phase  of 
AMSORB.  This  routine  serves  the  MS3  Obscurant  Dispersion  and  Obscura¬ 
tion  Models  for  print  output  and  for  saving  all  calculations  in  a  data 
file.  Program  control  is  returned  to  foreman  MS3  upon  completion  of  the 
routine . 


5.3.3  Foreman  BKC 

As  shown  in  Figure  5-1,  program  control  is  passed  to  foreman 
BF.C  from  Supervisor  APPLY.  Subroutine  UEC  contains  most  of  the  program 
logic  of  the  Battlefield  Environment  Routine.  This  subroutine  reads  the 
user  input  data,  determines  which  battlefield  sources  may  contribute  to 
the  observer-target  LOS  concentration,  computes  the  probability  of 
detecting  the  target,  and  prints  and  saves  the  calculation  results. 
Subroutines  TRNPRT,  DtFFUZ  and  ALLCN  are  called  by  subroutine  BEC. 
Figure  5-7  shows  the  relationship  of  these  subroutines  and  foreman  BEC. 

Subroutine  TRNPKT  has  the  same  prograte  logic  as  described  in 
Section  5.3.1.  it  is  used  by  subroutine  BEC  to  determine  the  boundary 
of  active  battlefield  sources  located  within  the  battlefield  grid  and  te 
compute  the  trajectory  of  clouds  from  battlefield  sources. 

Subroutine  DtFFUZ  has  the  same  logic  as  described  in  Section 
3.3.1.  it  Is  used  by  subroutine  BEL  to  account  for  dispersion  from 
sources  distant  from  the  iinc-of-slght  (see  Section  4.5). 

Subroutine  ALLCN  lus  the  same  program  logic  as  described  in 
Section  5.3.d.  It  is  called  by  subroutine  SEC  when  an  integrated 
concent  rat  ion  calculation  along  a  line  is  required.  Uhca  subroutine  bEC 
is  computing  the  cloud  brightness  in  the  probability  of  defection  calcu¬ 
lation,  ALLCN  is  also  called  to  compute  the  number,  position  and  concen¬ 
tration  of  the  cloud  elements  required  'to  approximate  the  cloud.  For 
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{  vppli cation  Phase  foreman  BEC  and  its  supporting  subroutines. 


quasi-continuous  sources,  ALLCN  calls  subroutine  CHIPT  to  calculate  the 
concentration  at  a  point  on  the  line-of -sight .  CHIPT  has  the  same  logic 
as  discussed  in  Section  5.3.2.  CHIIN  and  LNSRC  are  similar  in  purpose 
to  subroutine  CHIPT. 

5. A  GRAPHICS  PHASE  OF  AMS ORB 

Figures  5-8  and  5-9  show  block  diagrams  of  the  Graphics  Phase 
of  AMSGRB.  The  program  routine  PL'fPRG  is  the  Supervisor  that  controls 
all  graphics  output. 

For  c!.c  outputs  from  the  Mining -Payer  Analysis  Model  routine 
and  the  Trajectory/Transport  Dispersion  Model  application  routine,  PLTPRG 
retrieves  the  solutions  from  the  data  files  and  prints  and/or  plots  the  solu¬ 
tions  on  option  as  directed  by  the  user  through  XECAMS.,  »t>  shown  in  Figure 
5-8.  Options  include  the  printing  and  plotting  of  Intermediate  and  final 
wind-field  patterns  produced  by  the  mesoscale  model,  isopleths  of  con¬ 
centration  and  dosage,  cloud  trajectories  and  mixing-layer  depth  contours. 

The  plotted  solutions  are  shown  on  a  base  map  of  the  solution  matrix 
grid  with  a  background  of  terrain-height  contours.  The  subroutine 
labeled  iiKl.GOA  in  Figure  5-8  plots  and  labels  the  axes  of  each  plot  and 

the  terrain  height  contours  are  added  to  the  plot  by  subroutines  HEC005 
and  HECGG9.  Wind  vectors  are  calculated  and  plotted  by  subroutine 
Hgt‘0G6.  The  program  can,  in  subroutine  HECG08,  recalculate  (at  the 
user's  option)  additional  cloud  trajectories  for  other  source  locations 
using  the  wind-field  solution  in  the  AMSGRB  data  base.  The  supervisor 
PLTPRC  eaa  also  calculate  and  display  concentration  and  dosage  isopleths 
for  additional  lever*  at  the  user's  option.  This  feature  is  useful  if 
the  initially-requested  levels  do  not  show  details  in  which  the  user  is 
interested  or  if  the  initial  levels  are  either  so  large  or  so  small  that 
they  do  not  occur  in  the  solution  matrix.  It  should  be  noted  that 
PLTPRG  or  its  subroutines  do  not  recalculate  the  wind-field  nor  the 
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FIGURE  5-9.  Operational  Elements  of  PLTPRG  Graphics  Phase  for  the  MS3  Routine  Data  Output 


dispersion  parameters,  but  use  information  from  the  last  solutions  to 
recalculate  additional  trajectories  and  isopleths.  The  subroutine 
HEC007  controls  the  plotting  of  the  relative  variance  of  the  change  in 
layer  depth,  the  relative  variance  of  the  change  in  momentum,  and  the 
total  mass  and  total  momentum  as  a  function  of  model  solution  time.  As 
noted  in  Section  5.1,  the  relative  variance  of  the  change  in  layer  depth 
is  used  as  criteria  for  establishing  that  the  Mesoscale  Wind-Field  Model 
solution  has  converged.  Since  the  Mesoscale  Wind-Field  Model  theoreti¬ 
cally  conserves  mass  and  momentum,  the  plots  of  total  mass  and  momentum 
versus  model  time  are  useful  in  monitoring  the  effects  of  numerical 
approximations  in  calculations  of  the  wind-field  solutions. 

For  the  MS3  Routine  output,  PLTPRC  passes  program  control  to 
subroutine  SMKPLT.  Details  of  subroutine  SMKPLT  arc  shown  in  Figure  5- 
9.  Subroutine  SMKPLT  retrieves  the  output  data  from  data  file  SMKFIL 
and  passes  program  control  to  subroutine  SOMUP  along  with  the  plot 
option  chosen  by  the  user  through  PLTPRC.  Depending  on  the  plot  option 
chosen,  subioutine  SOMUP  passes  program  control  to  subordinate  routines 
when  necessary.  As  shown  in  Figure  5-9,  subroutine  DL1NE  plots  the 
solution  line  which  represents  either  line-of-sight  integrated  concen¬ 
tration,  probability  of  detection  or  transmittance.  Subroutine  PROBAX 
plots  the  probabllity-of-detection  axes,  subroutine  CHI AX  plots  the 
line-of-sight  integrated  concentration/ transmittance  axes  and  subroutine 
TIMEAX  plots  the  time  axes.  Program  control  is  passed  from  subroutine 
SOMUP  back  to  subroutine  SMKPLT,  and  then  from  SMKPLT  to  routine  PLTPRG 
when  all  plot  output  has  been  accomplished. 

At  the  completion  of  the  plotting  and  print  operations  from 
the  Craphics  Phase,  program  control  is  returned  from  routine  PLTPRC  to 
XKCAMS. 


SECTION  6 

EXAMPLE  CALCULATIONS 

The  Analysis  and  Application  Phases  of  AMSORB  arc  used  to  ill¬ 
ustrate  the  capability  of  the  battlefield  environment  components  (DEC) 
of  AMSORU  to  calculate  line-of-sight  integrated  concentrations  and 
probability  of  detection  for  meteorological  conditions  that  occurred 
in  the  vicinity  of  WSMH  during  14  November,  1974.  These  meteorological 
conditions  have  previously  been  used,  for  example,  to  illustrate  the 
capabilities  of  the  Mesoscalo  Wind-Field  Model  (Dumbauld  and  Bjorklund, 
1977)  and  the  MS 3  component  (Dumbauld,  Saterlie  and  Cheney,  1979). 

The  surface  weather  pattern  in  the  western  United  States  on 
Che  12th  of  November  was  characterized  by  a  high-pressure  ridge  extend¬ 
ing  from  the  Pacific  Northwest  southeastward  to  Texas.  The  high-pressure 
ridge  began  to  break  down  on  the  13th  of  November  with  the  rapid  approach 
of  a  low-pressure  system  moving  into  the  central  U.  S.  from  western 
Canada.  The  cold  front  associated  with  this  system  tended  to  remain 
east  of  the  Rocky  Mountains  as  it  pushed  south  and  east  during  the  day. 

By  the  morning  of  the  14th,  Che  low-pressure  cell  had  moved  into  the 
northeast  and  the  cold  front  extended  southwestward  through  the  eastern 
U.  S.,  into  Texas  and  northeast  Mexico,  Chen  abruptly  northward  along 
the  eastern  slopes  of  the  Rockies.  The  cold  front  had  actually  passed 
El  Paso  during  the  night  from  the  "back-door,"  or  moving  from  east  to 
west.  The  passage  was  not  accompanied  by  significant  weather  at  WSMR, 
but  broken  middle  and  high  clouds  were  reported.  The  high-pressure  sys¬ 
tem  behind  the  cold  front  had  moved  over  WSMR  by  the  morning  of  the  ISth 
of  November  and  calm  winds  u'U’e  reported  at  many  surface  stations  in  the 
area. 


The  temperature  profiles  from  rau-usonde  observations  made  at 
,E1  Paso  during  the  period  are  shown  in  Figure  6-1.  Inspection  of  Figure 
6-1  shows  that  the  cold  front  had  not  reached  El  Paso  by  1700  L  ou  the 
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13th  of  November.  By  0300  L  on  the  14th,  the  cold  front  had  passed  El 
Paso  and  the  temperature  inversion  at  about  800  millibars  associated 
with  the  frontal  surface  persisted  in  the  rawinsonde  observations  made 
at  1700  L  on  the  14th  and  0500  L  on  the  15th  of  November.  Surface 
observations  made  at  El  Paso  indicate  that  relatively  shallow  surface- 
based  inversions  formed  during  the  nighttime  and  were  eliminated  cy 
surface  heating  during  the  early  morning  hours.  The  mesoscale  wind- 
field  solution  for  0700  L  on  Che  14th  of  November  is  shown  in  Figure  6-2. 
Southeasterly  winds  dominate  the  flow  pattern  In  the  southern  portions 
of  the  Tularosa  Valley  and  the  Valley  Jornada  Dei  Muerto  west  of  the 
Organ  Mountains.  The  San  Andres,  Organ  and  Sacremento  Mountains 
extend  above  the  depth  of  the  surface  mixing  layer  and  form  an 
effective  barrier  to  the  southeasterly  flow,  diverting  tlte  wind  flow 
in  the  surface  layer  northward  in  Che  northern  portion  of  the  Tularosa 
Valley. 


Section  6.1  below  illustrates  cite  use  of  Che  Battlefield  Source 
Characteristics  Routine  (SORDAT)  to  define  a  bactlefield  scenario. 

Example  calculations  of  obscuration  information  for  this  bactlefield 
scenario  are  described  in  Section  6.2. 

6.1  EXAMPLE  BATTLEFIELD  SCENARIO 

The  Battlefield  Source  Characteristics  Routine  (SORDAT)  de¬ 
scribed  in  Section  4.2  is  used  to  generate  the  battlefield  scenario  in 
the  form  of  a  battlefield  source  location  and  characteristics  data 
base  file.  For  illustration  purposes,  we  have  used  the  smoke  and  dust 
sources  shoun  in  Table  6-1  for  our  calculations.  As  indicated  in 
Table  6-1,  we  have  used  4  smoke  munitions,  3  dust  clouds  from  exploding 
munitions,  burning  brush/ vegetation,  a  burning  building  and  6  muscle 
blasts  from  a  ISS-ata  artillery  piece.  The  last  souree  in  tins  table, 
labeled  with  a  source  identification  number  of  204.  is  used  with  the 
other  sources  in  a  second  problem  described  in  Section  6.2.  The  loca- 
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FI  CURE  6-2.  Mesosenle  uirnl-f  it-ld  taodel  solution  for  0700L, 
14  November  1974. 


TABI.E  6-1 

SOURCES  USED  IN  THE  EXAMPLE  CALCULATIONS 


Source 

Identification 

Number 

Source 

Category 

Source 

Description 

Time  of 
Release 
(hr, min, s) 

Source  coordinates 
x,  y,  z 
(m) 

101 

Smoke 

Munition 

155  mm  HC 

9:37:10 

380000,3636000,0 

102 

Smoke 

Munition 

105  ram  HC 

9:37:15 

374000,3639000,0 

103 

Smoke 

Munition 

155  mm  XM825 

WP  wedge 

i 

8:40:00  | 

376000,3613000,0 

104 

Smoke 

Munition 

155  mm  M110E2 
bulk  WP 

8:30:00  ! 

412000,3582000,0 

201 

Explosive 

Munition 

155  mm 
artillery 

i 

9:39:00 

376000,3639000,0 

202 

Explosive 

Munition 

155  ram 
artillery 

9:34:18 

377000,3643000,0 

203 

Explosive 

Munition 

8  in 

artillery 

9:21:46 

373000,3621000,0 

301 

Burning  Brush 
v/egetation 

User  Entered 
Burning  Brush 

8:53:00 

374000,3621000,0 

401 

Burning 

Building 

Burning 

Building 

9:38:00 

378000,3643000,3 

301 

Muzzle 

Blast 

155  mra 
artillery 

9:14:45 

386000,3632400,2 

502 

Muzzle 

Blast 

155  mia 
artillery 

9:37:19 

386000,3632395,2 

503 

Muzzle 

Blast 

155  mra 
artillery 

9:37:17 

386000,3632390,2 

504 

Muzzle 

Blast 

1  55  mra 
artillery 

9:37:15 

386000,3632385,2 

505 

Muzzle 

Blast 

155  mra 
artillery 

9:37:13 

386000 1 36  32  380 , 2 

S06 

Muzzle 

Blast 

155  mra 
artillery 

9:37:11 

386000,3632375.2 

204 

Explosive 

Munition 

15S  mra 
artillery 

9:39:19 

380300,3636300.0 

* 

coordinates  ai  the  centroid 
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cicm  of  these  sources  with  respect  to  the  line-of-sight  between  the 
observer,  target  and  background  is  shown  by  the  MX’'  marks  in  Figure  6-3. 

Figure  6-4  shows  the  card  input  data  to  SORDAT.  The  user 
instructions  for  completing  this  card  input  are  given  in  Section  A. 3  of 
Volume  H.  The  first  data  card  directs  the  AMSQRB  executive  routine 
XECAMS  to  call  SORDAT.  Because  the  executive  variable  TYPE3  (column  9 
of  the  first  data  card)  is  greater  than  zero,  SORDAT  expects  FORTRAN 
NAMELIST  input  data  using  name  QLST4.  The  QLST4  input  data  are  shown  in 
the  second  data  card.  The  remaining  input  data  cards,  except  the  last 
card,  describe  the  source  types,  locations,  release  times  and  other 
required  information.  For  example,  input  data  cards  3  through  5  in 
Figure  6-4  describe  source  number  101  in  Table  6-1,  the  155  ram  HC  smoke 
munition.  The  number  101  is  arbitrarily  assigned  as  the  source  identi¬ 
fication  number.  Following  the  instructions  for  smoke  munitions  in 
Section  A. 3.1  of  Volume  II,  the  number  in  column  7  identifies  this  card 
as  Che  first  card  for  this  munition.  The  number  1  in  column  11  sig¬ 
nifies  HC  smoke.  The  value  of  0  in  column  14  signifies  a  155  ram  HC 
munition.  Columns  16  through  20  contain  the  number  93710  indicating  the 
munitions  bursts  at  0937  hours  and  10  seconds  in  local  time.  Finally, 
the  UTM  x  and  y  location  of  the  burst  and  the  burst  height  are  given  in 
the  required  columns  of  card  number  1.  The  second  card  for  the  155  m» 

HC  smoke  projectile  shows  the  number  2  in  column  7,  indicating  this  is 
the  second  data  eard  for  this  source.  The  rest  of  the  columns  on  this 
card  are  blank  because  we  will  use  the  program  default  data  for  the  155  mm 
HC  munition.  The  third  card  for  this  souree  is  also  blank  except  for  the 
3  in  colua'>  7  indicating  the  third  data  card  for  this  source  and  the  value 
45  appearing  in  columns  39  and  40.  The  value  45  indicates  the  smoke 
projectile  approached  the  burst  location  with  an  angle  of  45  degrees  from 
grid  north  (towards  the  northeast).  The  remaining  data  cards  are  com¬ 
pleted  in  a  similar  manner  following  the  user  instructions  given  in 
Section  A. 3  of  Volume  11.  The  last  input  data  card  indicates  an  end  of 
data  to  the  computer  program.  The  SORDAT  routine  also  requires  the  mass 
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FIGURE  6-3 


Location  ef  the  smoke  and  dust  sourees  and  the  observer 
target-background  itnv-of-sighc  sed  in  Ehe  exaaple  eai 
eolations. 
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LIME  NO 


FIGURE  6-4.  (Continued) 


storage  data  base  file  SMKDAT  described  in  section  A. 8  of  Volume  II. 

The  data  input  example  in  Section  A. 7  of  Volume  II  ensures  that  SMKDAT 
is  assigned  in  the  runstream. 

The  optional  print  output  produced  by  SORDAT  is  shown  in 
Figures  6-5  and  6-6.  These  output  tables  contain  all  the  source  input 
information  required  by  the  dispersion  and  plume  rise  models.  For 
example,  the  exact  position  of  the  4  canisters  comprising  the  155  mm 
HC  smoke  munition  labeled  with  source  identification  number  101,  the  mass 
of  each  submunition,  the  burn  coefficients,  source  dimensions  (sigmas) 
and  heat  content  of  each  submunition  are  shown  in  Figure  6-5.  Figure  6-6 
contains  additional  information  for  the  same  submunitions.  Two  output 
tables  are  presented  because  not  all  information  could  be  included  in 
a  single  output  table.  Some  information  is  repeated  in  the  output  tables 
for  user  convenience.  Also,  for  user  convenience,  the  output  tables  are 
partitioned  into  different  categories  of  battlefield  sources,  i.e., 

■  moke  munition  sources,  burring  brush/vegetation  sources,  burning 
ocilding  sources,  etc. 

The  appropriate  dispersion  model  and  plume-rise  model  input 
data  are  also  written  to  the  battlefield  source  location  and  character¬ 
istics  data  base  file  BKCSOR  (described  in  Section  A. 8  of  Volume  II)  for 
subsequent  use  in  the  transport  and  dispersion  model  calculations. 


6.2  BATTLEFIELD  ENV1  KONMKNT  ROUT  INK  (RFC)  CALCULATIONS 

Two  solutions  of  BEC  are  described.  The  first  solution  is  for 
all  the  sources  shown  in  Table  6-1  except  the  source  identified  as  204 
and  the  second  solution  includes  204.  Figure  6-7  shows  the  card  input 
data  to  BEC.  The  first  data  card  directs  the  AHSORB  executive  routine 
XECAMS  to  call  I4KC  in  the  Applications  Pliase  of  AMSOUB  and  the  remaining 
cards  are  read  by  BEC.  As  discussed  in  Section  A. 4  of  Volume  U,  BEC 
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FIGURE  6-5.  (Continued) 
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FIGUPJi  6-C.  (Continued) 


FIGURE  6-7.  Example  card  input  data  to  EEC. 


uses  the  FORTRAN  NAMELIST  method  of  card  input  associated  with  NAMELIST 
name  QLST2.  The  third  card  in  the  sequence  specifies  the  time  the  user 
desires  to  obtain  the  solution,  either  in  terms  of  line-of-sight  inte¬ 
grated  concentration  or  obscuration  information.  Figure  6-7  shows  that 
a  solution  is  desired  for  0940  local  time.  The  fourth  card  in  the 

sequence  defines  the  coordinates  of  the  observer  (X  ,  Y  ,  Z  )  and 

0  0  0 

target  (Xfc ,  Y^,  Z^)  and  the  distance  DISBAK  between  the  target  and 
background  along  the  line  of  sight.  The  values  in  card  4  of  the  se¬ 
quence  correspond  to  the  locations  of  the  observer,  target  and  background 
shown  in  Figure  6-3.  Besides  card  input  data,  BEC  requires  mass  storage 
data  base  files  as  input.  In  this  case,  the  mesoscale  wind-field  solu¬ 
tion  shown  in  Figure  6-3  (from  file  MXLYRF) ,  along  with  file  BECSOR, 

SMKUAT  and  WSTRRN  were  used  in  the  calculations.  These  files  are  dis¬ 
cussed  in  Section  A. 8  of  Volume  II  and  the  example  runstream  for  assign¬ 
ing  these  files  is  given  in  Section  A. 7  of  Volume  II. 

Figure  6-8  shows  the  print  output  produced  by  BEC  for  the 
first  data  case  (source  204  omitted).  The  data  base  identl.^ers  are 
output  first  to  show  the  user  the  time  and  date  of  the  computer  run,  the 
solution  time,  the  time  and  date  of  the  mesoscale  wind-field  solution, 
and  the  date  and  time  the  battlefield  source  location  and  characteris¬ 
tics  data  base  file  was  gener  ced  using  SORDAT.  The  user  generated 
input  for  the  observer,  target,  and  background  locations  are  listed 
next,  followed  by  the  calculation  results. 

As  shown  in  Figure  6-8,  BFC  has  calculated  a  line-of-sight 

3 

integrated  concentration  (combined  dust  and  smoke)  of  1.41  x  10  milli¬ 
grams  per  square  meter  and  a  transmittance  in  the  visible  (.4  -.7  pa) 
of  0.108.  For  our  example  problea,  29  of  the  total  of  34  individual 
sources  were  considered  in  the  concentration  and  obscuration  calculations. 
That  is,  6  sources  were  excluded  fro«  the  calculation  either  because  they 
uere  not  released  at  a  time  that  allowed  thee  to  be  in  the  vicinity  of 
the  line-of-sight  at  the  calculation  time  or  because  the  sources  are 
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located  in  regions  where  winds  will  not  carry  the  material  through  the 
line-of-sight  (for  example,  sources  202  and  401  in  Figure  6-3).  Figure 

6-8  shows  that  the  calculated  target  brightness  is  620  candles  per 

-2  -2 
square  meter  (cd  m  ) ,  the  background  brightness  is  345  cd  m  ,  and  the 

observer  to  target  and  observer  to  background  brightness  of  the  obscur- 

_2 

ant  cloud  are  948  cd  m  .  The  calculated  target-background  contrast  is 
0.21  for  our  example  calculation,  which  yields  a  probability  of  detection 
equal  to  or  greater  than  95  percent.  The  output  also  indicates  that  all 
visibility  model  computations  were  performed  and  that  there  were  11 
cloud  elements  along  the  line-of-sight  used  in  the  obscuration  calcula¬ 
tions.  In  addition  to  the  printed  output,  BEC  outputs  results  of  the 
calculations  and  model  inputs  to  the  mass  storage  data  base  file  BECSOL, 
discussed  in  more  detail  in  Section  A. 8  of  Volume  II.  For  this  data  case, 
the  BEC  calculations  required  290  seconds  of  UNIVAC  time. 

Figure  6-S  3hows  the  results  of  calculations  made  for  all  the 
sources  described  in  Table  6-1,  including  the  source  identified  by  the 
number  204.  In  this  case  the  line-of-sight  integrated  concentration  is 
4.7  x  10^  milligram  per  square  meter.  Since  the  sum  of  the  products  of 
the  line-of-sight  integrated  concentration  of  each  pollutant  and  its' 
respective  coefficient  of  attenuation  in  the  visible  (.4  -  .7  ym)  is 
greater  than  11  (see  Section  2.4.2)  for  the  user  specified  calculation 
time  with  the  addition  of  the  source  identified  as  204,  AMSGRB  does  not 
make  the  time  consuming  brightness  calculations.  Instead  the  trans¬ 
mittance  is  automatically  set  to  zero  and  the  probability  of  detection 
is  less  than  5  percent.  For  this  data  case,  the  computation  time  is  a 
little  greater  than  2  seconds,  or  more  chan  a  factor  of  10  less  than  the 
calculation  time  required  when  the  brightness  calculations  had  to  be 
made  when  the  source  identified  by  204  was  omitted  from  the  scenario. 
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FIGURE  6-S .  Example  print  output  produced  by  BEC  for  the  second  data  case 
(source  204  included). 
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